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Abstract 

We review the quantum interference effects in a system of interacting electrons con- 
fined to a quantum dot. The review starts with a description of an isolated quantum 
dot. We discuss the Random Matrix theory (RMT) of the one-electron states in the 
dot, present the universal form of the interaction Hamiltonian compatible with the 
RMT, and derive the leading corrections to the universal interaction Hamiltonian. 
Next, we discuss a theoretical description of a dot connected to leads via point 
contacts. Having established the theoretical framework to describe such an open 
system, we discuss its transport and thermodynamic properties. We review the evo- 
lution of the transport properties with the increase of the contact conductances 
from small values to values ~ /ttTl. In the discussion of transport, the emphasis 
is put on mesoscopic fluctuations and the Kondo effect in the conductance. 

PACS numbers: 73.23.Hk, 73.23.-b, 72.10.Bg 

Key words: Coulomb Blockade; Mesoscopic Fluctuations; Random Matrix 
Theory; Bosonization; Kondo effect 



Preprint submitted to Elsevier Preprint 1 February 2008 



Contents 



1 


Introduction 








2 


The model 





2.1 Non- interacting electrons in an isolated dot: Status of RMT 





2.2 Effect of a weak magnetic field on statistics of one-electron states in 


the dot 





2.3 Interaction between electrons: The universal description 



2.3.1 Universal description for the case of short-range interaction 



2.3.2 RMT for the intra-dot Coulomb interaction 



2.3.3 Final form of the effective Hamiltonian 



2.4 Inclusion of the leads 



3 Strongly blockaded quantum dots 



3.1 Mesoscopic fiuctuations of Coulomb blockade peaks 



3.2 Mesoscopic fiuctuations of Coulomb blockade valley^ 



3.3 Kondo effect in a strongly blockaded dot 



3.4 Overall temperature and gate voltage dependence of the conductance 



4 Weakly blockaded dots 



4.1 Main results for the conductance and differential capacitance of an 
almost open quantum dot 



4.2 Finite-size open dot: Introduction to the bosonization technique and 
the relevant energy scales) 





4.3 The limit of low temperature: The effective exchange Hamiltonian 


and Kondo effect 





4.4 Simplified model for mesoscopic fluctuations! 





4.5 Rigorous theory of mesoscopic fluctuations: formalism foi 


multichannel junctions 





4.5.1 Derivation of the effective action 



4.5.2 Bosonization of the effective Hamiltonian 



4.6 Results for the differential capacitance 



118 



4?7 Conductance through a dot with two almost reflectionless junctions 127 

132 



4.8 Conductance in the strongly asymmetric setup 



Conclusions 



Appendices 



"K Correlation of the wave-functions in ballistic dots 



B Effect of the Interaction in the Cooper channel 
C Scattering states and derivation of Eq. (99) 



D Mesoscopic fluctuations of elastic co-tunneling far from the peaks] 
E Derivation of the Hamiltonian (238) 



F Canonical versus grand canonical ensembles 



G Correlation functions 



H Derivation of Eqs. (340) - (343) 



Derivation of Eqs. (360) - (361) 



References 



138 
141 
141 
143 
144 
147 
149 
151 
153 
156 
160 
162 



3 



1 Introduction 



Conventionally, electric transport in bulk materials is characterized by the 
conductivity a. Then, the conductance G of a finite-size sample of dimen- 
sions Lx X Ly X Lz can be found by combining the conductances of its smaller 
parts, G = aLyLz/ Lx- This description, however, is applicable only at suffi- 
ciently high temperatures, at which the conductivity can be treated as a local 
quantity. It was discovered about two decades ago [1,2] that the quantum 
corrections to the conductivity are non-local on the scale of the (temperature 
dependent) dephasing length L^, which is much larger than the elastic mean 
free path. If the sample is small, or the temperature low, so that exceeds 
the sample size, the concept of conductivity loses its meaning. Moreover, the 
conductance G acquires significant samplc-to-sample fiuctuations and can no 
longer be treated as a self- averaging quantity. For a temperature T below the 
Thouless energy Et — Tl/td-, these sample-to-sample fluctuations are of the 
order of the conductance quantum e^/27r7l, independent of sample size or mean 
free path {tjj is the time it takes for an electron to diffuse through the system). 
This phenomenon is called "Universal Conductance Fluctuations" (UCF), and 
it was studied in great detail both theoretically and experimentally, see [3-7] 
for a review. It is important to emphasize that the lack of self-averaging of 
the conductance is a feature of mesoscopic samples, i.e., samples with linear 
dimensions smaller than but still much larger than the Fermi wave length. 

These ideas that defined the field of mesoscopic physics were initially devel- 
oped in the context of transport through disordered metals. A convenient way 
of creating more controllable mesoscopic samples was developed with the use 
of semiconductor heterostructures, see [8] for a review. As the electron mo- 
tion in these systems is quantized in the direction perpendicular to the plane 
of the hctcrostructurc, a two-dimensional electron gas is formed at the inter- 
face between the layers. With the help of additional electrostatic confinement, 
usually in the form of metal gates deposited on top of the hetrostructure, 
the two-dimensional electron gas (2DEG) can be tailored to form a finite-size 
sample, referred to as a quantum dot. Its size, shape, and connection with 
the rest of the 2DEG can be controlled by means of a sophisticated system of 
electrostatic gates, see Fig. 1. 

Links connecting the quantum dot to the 2DEG (quantum point contacts) are 
characterized by the number A^ch of electron modes (or channels) propagating 
through the contact at the Fermi level, and by the set of transmission proba- 
bilities for these modes. With the increase of the point contact cross-section, 
the number of modes, Ach, is increasing. If the boundaries of the point contact 
are smooth on the length scale set by the electron wavelength, then the contact 
acts as an electron wave guide [10]: all the propagating modes, except the one 
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GaAs AL^Ga^.^As 



Fig. 1. Micrograph of the quantum dot [9]. By changing the voltages on the gates 
conductances of point contacts and the shape of the dot can be changed. 

which opened last, do not experience any appreciable backscattering.|3] For 
the last mode, the transmission coefficient varies from to 1 in the crossover 
from the evanescent to propagating behavior [5] . In the specific case of a gate- 
controlled contact, see Fig. 1, the distance between the gate and the depleted 
edge of the 2DEG is typically large compared to the electron wavelength in 
the electron gas [12]. The curvature radius of the edge of the 2DEG at the 
contact is of the order of this distance, and therefore the smoothness condition 
is satisfied. 

A statistical ensemble of quantum dots can be obtained by slightly varying 
the shape of the quantum dot, the Fermi energy, or the magnetic field, keeping 
the properties of its contacts to the outside world fixed. Statistical properties 
of the conductance for such an ensemble depend on the number of modes, 
A^ch; in the junctions and their transparency. The "usual" UCF theory ade- 
quately describes the conductance through a disordered or chaotic dot in the 
limit of large number of modes, Nch 3> 1. When all channels in the two point 
contacts are transparent, the universal value of the root mean square (r.m.s.) 
conductance fiuctuations ~ e'^/2nh is much smaller than the average 

conductance {G) ~ Ncte"^ f^nh. The distribution of conductances is essentially 
Gaussian, and can be found theoretically by means of the standard diagram- 
matic expansion for disordered systems [13], or with Random Matrix Theory 
(RMT), see Refs. [14] and [15] for reviews. Upon closing the point contacts 
{i.e., decreasing Nch), the average conductance decreases whereas the r.m.s. 
fiuctuations retains its universal value. At N^h — 1, the fiuctuations and the 
average of the conductance are of the same order, and the conductance dis- 
tribution function is no longer Gaussian. Still, if one neglects the interaction 
between the electrons, the statistics of the conductance can be readily obtained 
within RMT [14,16,17]. 

^ We refer the reader to Ref. [11] for a pedagogical presentation of the theory of 
adiabatic electron transport 
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The RMT of transport through quantum dots can be easily justified for non- 
interacting electrons. However, real electrons are charged, and thus interact 
with each other. It is known from the theory of bulk disordered systems, that 
the electron-electron interaction results in corrections to the conductance [18] 
of the order of unit conductance quantum e^/27r^. For a quantum dot con- 
nected to a 2DEG by single-mode junctions, this correction is not small any- 
more in comparison to the average conductance {G). 

We see that at (G) ~ 6^/271% and low temperatures, both mesoscopic fiuctua- 
tions and interaction effects become strong. Quantum interference and interac- 
tion effects can not be separated, and there is no obvious expansion parameter 
to use for building a theory describing the statistics of the conductance in this 
regime. This paper reviews the current understanding of physical phenomena 
in a wide class of mesoscopic systems - quantum dots - under the conditions 
in which effects of electron-electron interaction are strong. 

A first step to approach this problem is to disregard the mesoscopic fiuctua- 
tions altogether. It means that one neglects the random spatial structure of 
the wave functions in the dot and the randomness of the energy spectrum; 
the discrete energy spectrum with spacing A between the one-electron levels 
in a closed dot is replaced by a continuous spectrum with the correspond- 
ing macroscopic density of states. For reahstic systems, this approach can be 
justified sometimes at not-so-low temperatures. The remaining problem of ac- 
counting for the Coulomb interaction is non-perturbative at small conductance 
G ^ e^/27r/i, and therefore still not trivial. It constitutes the essence of the 
so-called Coulomb blockade phenomenon [19,8]. 

The Coulomb blockade manifests itself most profoundly (see, e.g., Ref. [20,8]) 
in oscillations of the dot's conductance GiVg) with the variation of the voltage 
on a gate, which is capacitively coupled to the dot, see Fig. 2. The resulting 
dependence GiVg) exhibits equidistant Coulomb blockade peaks separated by 
deep minima (Coulomb blockade valleys). The peaks occur at the charge de- 
generacy points, i.e., specific values of Vg at which changing the dot charge by 
a single quantum e does not cost any energy. The idea of the Coulomb blockade 
was suggested in the early experimental paper [22] in the late 1960's, though 
the term "Coulomb blockade" was coined only two decades later [23]. A quan- 
titative theory in terms of rate equations describing the transport through 
a blockaded quantum dot or metal grain at G <^ e^/27r/i, was formulated 
in Ref. [24,25], and was generalized to systems with a controllable gate in 
Ref. [23] . This theory is commonly referred to as the "orthodox theory" . The 
main conclusion of the orthodox theory is that the conductance through a 
blockaded grain at low temperatures is exponentially suppressed. The first 
experiments on gated metallic Coulomb blockade systems were reported in 
Ref. [26]. 
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Fig. 2. (a): Schematic view of the dot with the gate electrode; (b): conductance of 
the dot as a function of the gate voltage (from Ref. [21]). 



At larger values of the dot's conductance (corresponding to larger values of the 
conductances Gi and G2 of the point contacts connecting the dot to the out- 
side world), the orthodox theory and, eventually, the phenomenon of Coulomb 
blockade break down because of quantum fluctuations of the charge in the dot. 
The fluctuations destroy the Coulomb blockade at the value G ~ e'^/nh. Quan- 
tum fluctuations affect both transport properties (such as linear conductance 
and I-V characteristics) and thermodynamics {e.g., differential capacitance). 
The I-V characteristic in the Coulomb blockade valleys at relatively small 
conductance was considered in [27] in second order perturbation theory in G; 
later the corresponding non-linear contribution to the current found in [27] 
became known as "inelastic co-tunneling" [28]. Corrections induced by charge 
fluctuations at G ^ e'^/nh were calculated in [29] for the thermodynamic char- 
acteristics, and in [30,31,40,32] for the linear conductance. At G <^ e^/nh the 
problem of quantum charge fluctuations was mapped [33] on the multi-channel 
Kondo problem [34,35], which in some cases can be solved exactly [36,37]. 

Another limit which allows for a comprehensive solution, is the case of a 
point contact with one almost open conducting channel linking the dot to 
the lead (conductance G is close to the conductance quantum) [38-40]. In 
this limit, and still neglecting mesoscopic fluctuations (A = 0), the Coulomb 
blockade almost vanishes: the contribution periodic in the gate voltage Vg to 
the differential capacitance, which is the prime signature of Coulomb blockade, 
is small compared to the average capacitance [39]. Moreover, it was shown 
within this model that the Coulomb blockade disappears completely in the 
case where the channel is fully transparent. 

In all of the aforementioned works the effects of mesoscopic fluctuations, as- 
sociated with the finite size and the discrete energy spectrum of the isolated 
dot (nonzero level spacing A), were ignored. However, these fluctuations dras- 
tically affect the results. In open dots, mesoscopic fluctuations reinstate peri- 
odic oscillations of the differential capacitance and conductance as a function 
of the gate voltage V^, as a remnant the charge discreteness in the dot, even 
in the limit of completely transparent channels in the point contacts. If the 
dot is almost closed, the chaotic nature of the wave functions in the dot gives 
rise to mesoscopic fluctuations of observable quantities. Both in the Coulomb 
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blockade valleys and at the peaks, the conductance turns out to have large 
fluctuations, and is very sensitive to a magnetic fleld. 

Our main goal in this paper is to review the rich variety of effects associ- 
ated with the flnite level spacing A, the chaotic nature of the wave functions, 
and the Coulomb interaction for both closed and open dots. Here is a one- 
paragraph guide to the review. Sections 2, 3, and 4 are devoted, respectively, 
to the physics of an isolated quantum dot, transport through a dot weakly 
connected to leads, and thermodynamic and transport properties of a dot 
connected to lead(s) by one or two almost reflectionless quantum point con- 
tacts. We begin in Section 2 with the formulation of the model describing the 
statistical properties of the system. We clarify the conditions of the equiva- 
lence of the usual diagrammatic technique and Random Matrix representa- 
tion, and establish the validity of the capacitive interaction approximation. 
The central result of this part can be found in Subsection 2.3.3, where we 
present the universal form of the interaction Hamiltonian describing the dot, 
and estimate the accuracy of the universal description. In Section 3 we re- 
view the theory of mesoscopic fluctuations of the conductance in peaks and 
valleys of the Coulomb blockade,^ and the theory of the Kondo effect in quan- 
tum dots. Weak dot-lead coupling allows us to use the familiar technique of 
the perturbation theory in tunneling strength. We believe it makes Section 3 
easy to read. Nevertheless, a very pragmatic reader who wants to grasp the 
main results flrst, without going through the details, may start by reading the 
summary in Subsection 3.4. Section 4 deals with the similar phenomena in 
quantum dots connected to leads by transparent or almost transparent point 
contacts, where the interplay between quantum charge fluctuations and meso- 
scopic physics becomes especially interesting. Methods described and used in 
Section 4 are technically more involved than in Section 3, and include the 
bosonization technique along with the effective action formalism. To make the 
navigation through the material easier, we decided to start this Section with 
a presentation of the main results for the conductance through a dot and for 
its thermodynamics, see Subsection 4.1. Along with providing the general pic- 
ture for the conductance and differential capacitance, this Subsection points 
to formulas valid in various important limiting cases, which are derived later 
in Subsections 4.6-4.7. Subsections 4.2-4.4 review the main physical ideas 
built into the theory, while the details of the rigorous theory are presented in 
Subsection 4.5. Last but not least, we illustrate the application of theoretical 
results by briefly reviewing the existing experimental material, see Subsec- 
tions 3.1-3.3, and 4.8. 



^ There is certain overlap with Ref. [15] 
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2 The model 



In Subsections 2.1 and 2.2 we discuss the electronic properties of an isolated 
quantum dot. Under the assumption that electron motion inside the dot is 
chaotic, we establish a hierarchy of energy scales for the free-electron spectrum: 
Fermi energy, Thouless energy, and level spacing. The correlation between 
electronic eigenstates and eigenvalues can be described by Random Matrix 
Theory (RMT) if the difference between the corresponding energy eigenvalues 
is smaller than the Thouless energy. In subsection 2.3 we proceed with the 
inclusion of effects of the interaction between the electrons. It is shown that 
the two-particle interaction matrix elements also have hierarchical structure. 
The largest matrix elements correspond to charging of the quantum dot, as 
described by the charging energy in the constant interaction model. Finally, 
open quantum dots are discussed in Subsection 2.4, where we introduce the 
Hamiltonian of the junctions connecting the dot with the leads and relate the 
parameters of this Hamiltonian to the characteristics of the quantum point 
contacts. 



2. 1 Non-interacting electrons in an isolated dot: Status of RMT 



Let us start from a picture of non-interacting electrons, postponing the intro- 
duction of electron-electron interactions to Sec. 2.3. The dot is described by 
the Hamiltonian 



Hf= dr 



—V%i)^Vip + U ir) "ip^ 
2m 



where the electron creation (annihilation) operators obey canonical anticom- 
mutation relations 

{^(fi),^(r2)} =0, {^^(fi),^(f2)} = 5(ri-r2). (2) 



The potential U{r) describes the confinement of electrons to the dot, as well as 
the random potential (if any) inside the dot.0 We will neglect the spin-orbit 
interaction; The one-particle Hamiltonian then is diagonal in spin space, and 
for the time being we omit the spin indices. 



We put ?i = 1 in all the intermediate formulae. 
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The eigenfunctions of the Hamihonian (1) arc Slater determinants built from 
single-particle states with wavefunctions defined by the Schroedinger equation 



V2 



(3) 



In this orthonormal basis 

HF = T.^Jii'a, (4) 



where the fermionic operators ipa are defined as ipa = f drip (r) 0* (r) and 
have the usual anticommutation relations, cf. Eq. (2). 

Each particular eigenstate depends sensitively on the details of the random 
potential U{r), which is determined by the shape of the quantum dot. How- 
ever, we will not be interested in the precise value of observables that depend 
on the detailed realization of the potential U{r). Instead, our goal is a statis- 
tical description of the various response functions of the system with respect 
to external parameters such as magnetic field, gate voltage, etc., and of the 
correlations between the response functions at different values of those param- 
eters. Hereto, the statistical properties of a response function are first related 
to the correlation function of the eigenstates of the Hamiltonian (3). Then, we 
can employ the known results for statistics of the eigenvalues and eigenvectors 
in a disordered or a chaotic system [41-44]. 

For a disordered dot, the correlation functions are found by an average of the 
proper quantities over the realizations of the random potential. Such averaging 
can be done by means of the standard diagrammatic technique for the electron 
Green functions 

a'''''fer-:.f.) = i:«§i^, (5) 



where plus and minus signs correspond to the retarded (R) and advanced (A) 
Green functions respectively. 

The spectrum of one-electron energies is fully characterized by the density of 

states 

u{e)^Y.5{e-e^) = ^J df [g\e, f, r) - g'^ie, r, f^] , (6) 



where the last equality follows immediately from the definition (5). Even 
though the density of states is a strongly oscillating function of energy, its 



10 



average is smooth, 



(7) 



where A(£:) is the mean one-electron level spacing. It varies on the charac- 
teristic scale of the order of the Fermi energy Ep, measured, say, from the 
conduction band edge, which is the largest energy scale in the problem. Since 
we are interested in quantities associated with a much smaller energy scale, wc 
can neglect the energy dependence of the mean level spacing. The average in 
Eq. (7), denoted by brackets (...), is performed over the different realizations 
of the random potential for the case of a disordered dot, or over an energy 
strip of width ^ A, but <^ Ep for a clean {i.e., ballistic) system. 

The average density of states does not carry any information about the corre- 
lations between the energies of different eigenstates. Such information is con- 
tained in the correlation functions for the electron energy spectrum. Probably, 
the most important example is the two-point correlation function TZ^'^\uj), 

n^^\uj) = A^iy{e)iy{e + cu)) - 1. (8) 

To calculate TZ'^'^^lu), one substitutes the expression (6) for the density of 
states in terms of the Green functions and G"^ into Eq. (8), 

n^'Koj) = ^Re J dr^dr^iQ'' {e + u;, Q"" {e, r^, ^2)) - 1- (9) 

When the dot contains many electrons — which is the case of interest here 
— the size of the dot and the region of integration in Eq. (9) greatly exceed 
the Fermi wavelength. In that case, averaging of the products Q^Q^ using 
the diagrammatic technique yields two contributions to the integrand in (9), 
proportional to the squares of diffuson and Cooperon propagators in coinciding 
points, respectively, see Fig. 3. The result of such a calculation [41], 

7e(^)M = ^ReE 7-^-2, (10) 

is expressed in terms of the eigenvalues 7„ of the classical diffusion operator, 

-DVVn(f)=7n/n(r), (11) 

supplemented by von Neumann boundary conditions at the boundary of the 
dot. Here D is the electron diffusion coefficient in the dot, and the Dyson 
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(a) (b) 




Fig. 3. Diagrammatic expansion of the correlation of the density of states for the 
disordered system [41]. In the presence of magnetic field (unitary ensemble, (3 = 2) 
the Cooperon diagram (b) is suppressed. In the figure, the symbols P and Q refer 
to the momenta of the electrons; P is the large momentum, \P — kp] <C kp, Q is the 
small momentum \Q\ <C kp; ^ = P'^ /2m — /U, where /i is the chemical potential; in 
two dimensions, the mean free time r is related to the diffusion constant D through 
D = VpT/2, vp being the Fermi velocity. 

symmetry parameter /5 = 1 (2) in the presence (absence) of time reversal 
symmetry. Ensembles of random systems possessing and not possessing this 
symmetry are called Orthogonal and Unitary ensembles, respectively. Being 
derived by means of diagrammatic perturbation theory, Eq. (10) and Eq. (12) 
below are valid for energy differences uj ^ A only. The exact results, valid for 
all u, are also available, see, e.g., Refs. [43-47]; however, the approximation 
(10) is sufficient for the present discussion. An expression similar to Eq. (10) 
is believed to be valid for a chaotic and ballistic quantum dot. The only differ- 
ence with the diffusive case is that instead of the eigenvalues of the diffusion 
operator one has to use eigenvalues 7„ of a more general (Perron-Frobenius) 
operatorQ of the classical relaxational dynamics [49] ;0 in this case the eigen- 
values 7„ can be complex with Re7„ > 0, see also Appendix A. 



The corresponding non-linear cr- model was first suggested in Ref. [48], however, 
this classical (Perron-Frobenius) operator was erroneously identified with the Liou- 
ville operator. 

^ This identification does not properly handle the repetitions in periodic orbits [50] , 
and therefore is applicable only for the systems, where all the periodic orbits are 
unstable. 
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The relation of TZ^'^\uj) to the dynamics of a particle in a closed volume allows 
one to draw an important conclusion regarding the universality of 7l^'^\uj). 
Note that a spatially uniform particle density satisfies the relaxational dy- 
namics or the diffusion equation. Because of particle number conservation, 
such a solution is time independent. Therefore, the lowest eigenvalue 70 in 
Eq. (11) must be zero, independent of the shape of the dot and the strength 
of the disorder, 

= + #^Rc ^ ^ ,2 ■ (12) 



The first term in Eq. (12) is universal. One can immediately see that if hcu is 
much smaller than the Thouless energy Ex = ^Re7i, the remaining terms in 
Eq. (12) are small compared to the first one and can be neglected. 

The universal part of the correlation function (12) can be reproduced within 
a model where the Hamiltonian (1) is replaced by a hermitian matrix with 
random entries 



The coefficients Haj in Eq. (13) form a real {(3 — 1) or complex (/3 = 2) random 
Hermitian matrix Ti of size M x M, belonging to the Gaussian ensemble 

{Ha^Ha'^') — 5— [(^a7'(^a'7 + (2//? — 1) (^aa'(^77'] , = 1, 2. (14) 

TT 



The Hamiltonian (13), with the distribution of matrix elements (14) repro- 
duces the universal part of the spectral statistics of the microscopic Hamilto- 
nian (1) in the limit M ^ oo of large matrix size. The condition A ^ Et-, 
provides the significant region uj <^ Et where the non-universal part of the 
two level correlator 7^2 (f-^) could be neglected, and where the replacement of 
the microscopic Hamiltonian by the matrix model (13), is meaningful. This 
condition may be reformulated as the requirement that the dimensionless con- 
ductance g of the dot be large, where 

_ n^i Et , . 

^=^ = ^»1- (15) 



A large value of g indicates that the dot can be treated as a good conductor. 

We now present explicit expressions for the dimensionless conductance g for 
a disc-shaped dot for two simple models: (i) a dot of radius it! exceeding the 
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electron mean free path /, and (ii) a ballistic dot of radius R with a boundary 
which scatters electrons diffusively, see Ref. [51,52] and references therein. 

For a diffusive dot, one needs to solve Eq. (11) to obtain 71 = Dx^/R^, 
where Xi ~ 2.40 is the first zero [53] of the Bessel function, Jq{xi) = 0. 
Using the electron density of states in two-dimensional electron gas, one finds 
A = 2h^/mR^, and the conductance is 

Xi 27lh Xi 

Here kp is the Fermi wave vector, and Rn is the resistance per square of the 
two-dimensional electron gas. Note that g is independent of the radius of the 
dot in the case of the diffusive electron motion. 

For the model (ii), the eigenvalue Re 71 ^ 0.38vf/R of the Liouville operator 
with the diffusive boundary conditions was found in Ref. [51], and one finds 

g = O.lQkpR. 

2.2 Effect of a weak magnetic field on statistics of one-electron states in the 
dot 

Physical properties of a mesoscopic system are manifestly random. The statis- 
tics of the random behavior can be studied by a measurement of the fiuctu- 
ations caused by the application of an external magnetic field B. The role 
of such a field is primarily to alter the quantum interference pattern. As a 
result, the characteristic fields which significantly affect the properties of a 
mesoscopic system are rather weak in a classical sense: the cyclotron radius of 
an electron remains much larger than the linear size of the dot, and the effect 
of the magnetic fields on classical trajectories can be neglected. 

A magnetic field can be included in the random matrix model of Eq. (13). 
In order to include the magnetic field, we have to take the matrix Hmn in 
Eq. (13) from a crossover ensemble that interpolates between the orthogonal 
{j3 = 1) and unitary (/3 = 2) ensembles [54-56]: 

Ha^ = (1 + X^) (h^^ + iXni^) . (16) 

Here Ti.^^ is the random realization of independent Gaussian real symmetric 
(antisymmetric) M x M matrices, 

{W^fir}^,) = [5^o.'5,,, ± 5^,>5^,,] , (17) 
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where "+" ("— ") sign corresponds to the symmetric (antysimmetric) part 
of the Hamihonian. In Eq. (16), X is a real parameter proportional to the 
magnetic field (the precise relation between X and B is given below), so that 
the time reversal symmetry is preserved, 

Ha-riB) — H-yai—B). (18) 



The prefactor (1 + X^) in Eq. (16) is chosen in such a way that the mean 
level spacing A remains unafi^ected by the magnetic field. 

The correlation function of elements of the Hamiltonian (16) at different values 
of the magnetic field Bi^2 can be conveniently written similarly to Eq. (14) as 



{T-Ca'y{Bi)Ha''y' {B2)) 



(19) 



X 



D C 

The dimensionless quantities N^^ ' , which characterize the effect of the mag- 
netic field on the wave-functions and the spectrum of the closed dot, are related 
to the original parameters X as 



4M, 



AM 



where the normalization is chosen in a way that the size of the matrix M does 
not enter into physical quantities in the physically relevant limit M ^ 1. 

The relationship between the parameter X and the real magnetic field applied 
to the dot is best expressed with the help of and A^^, 



where 5? 3> 1 is the dimensionless conductance of the closed dot, see Eq. 
(15), $1(2) is the magnetic flux through the dot at the magnetic held -61(2), 
$0 = hc/e is the flux quantum, and x is a geometry dependent numerical 
factor of order unity. A derivation of Eq. (20) is sketched at the end of this 
subsection. We also refer the reader to Refs. [57,58]. Here we give the values 
of the constant x ^ot two speciflc cases: (i) % = tt/xi ~ 1.31 for a disk-shaped 
diffusive dot, and (ii) x ~ 2.78 for a disc-shaped baUistic dot with diffusive 
boundary scattering [52]. 

The random matrix description (16) — (19) is valid provided X/f'*" -C g. At 
larger values of the magnetic fleld, a universal description is no longer appli- 
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cable. However, for most physical quantities, the crossover between the or- 

D C 

thogonal and unitary ensembles takes place at A^^^ ' 1, i.e., well within the 
regime where the RMT description is appropriate. 

A small magnetic field significantly affects the correlations of the eigenvectors 
and eigenvalues of the random Hamiltonian (16). For the pure orthogonal and 
unitary ensembles {[3 = 1 and 2, respectively), the eigenvectors and eigenvalues 
are independent of each other. Moreover, in the limit M — > cxd, the eigenvectors 
of different eigenstates are independently distributed Gaussian variables, 

iUmj)) - - (I - 1) (21) 



In other words, an eigenvector can be represented as 



where M is the size of the matrix, u, v are independent real Gaussian variables, 
{unUm) = ^nm, {vnVm) = Snm, {unVm) = 0, and ta = (l) for the orthogonal 
(unitary) ensemble [59,54]. Equation (22) is a consequence of the symmetry of 
the distribution of the random matrices with respect to an arbitrary rotation 
of the basis. At the crossover between those two ensembles, this symmetry no 
longer exists, which leads to a departure of the distribution of il^ain) from a 
(real or complex) Gaussian [60,61], and to correlations of the values of the 
wave functions at different sites [62,63]. It was noticed in Refs. [64,65] that 
the crossover between the orthogonal and unitary ensembles can be described 
by using the decomposition (22) where the parameter ta is no longer fixed at 
the extremal values t = or i = 1, but a fluctuating quantity < i < 1 with 
the distribution function [60-63] 



W{t) 



1-t' 



X 



exp 


- (to 


l + t^Y 4 


2t 





(f){x) = / dyexp —x(^l — y^^ 



(23) 



where Nj^ is given by Eq. (20) with $i = $2 = being the magnetic flux 
through the dot. We see, that the crossover between two ensembles occurs 
at characteristic scale of the magnetic held ~ 1, or using (20), at $ = 
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$o/(2y'X^). Equation (23) completely describes the distribution of a single 
eigenvector. In order to find the statistics of the quantities contributed by 
several levels, it is necessary to find the joint distribution of of several 
levels, W{ti,t2, . . .). The exact form of such a distribution function is not 
known. 

In the opposite limiting case when a physical result is contributed to by a large 
number of levels, it is possible to use the conventional diagrammatic technique 
[13] to calculate different moments of the Green functions 



see also Eq. (5). The diagrammatic calculation shown in Fig. 4 yields 



(25) 

where Nj^' and N/^ are defined in Eq. (19). The irreducible averages (Q^Q^) or 
{Q^Q^) are smaller than Eq. (25) by factor of uj/{MA) and can be disregarded 
in the thermodynamic limit M oo. All the higher moments can be expressed 
in terms of the second moments (25) with the help of the Wick theorem.j^ 

Equation (25) can serve as a starting point to relate the quantities A^^f, Ny^ to 
the real magnetic field, thus providing a derivation of Eq. (20). To establish 
such a relation, one needs to recalculate the average (25) using the microscopic 
Hamiltonian (1). The result is again given by Eq. (10), where now the magnetic 
field is introduced into the diffusion equation of Eq. (11) by the replacement 
of d/dfhy the covariant derivative d / dr + i{e / ch) [r) , and a modification of 
the boundary condition to ensure that the particle flux through the boundary 
remains zero. The vector-potentials A^{r) are related to the magnetic fields 
Bi and B2 by curlA='= = Bi ± B2- The corresponding lowest eigenvalues % 
and 70" are identified with (A/27r)A^/P and (A/27r)A^^ respectively. A similar 
replacement can be also made in the Liouville equation for the ballistic system. 

We close this subsection with a more elaborate discussion of Eq. (25) and its 
physical interpretation. Hereto, we consider the case of equal magnetic fields 
Bi and i?2, ^-e., Nj^ = 0. Then, one sees from Eq. (25) that a weak magnetic 



" It is noteworthy, that such a decomposition is legitimate only at energy scales 
greatly exceeding the level spacing A. It is not valid for the calculation of the 
properties of a single wavefunction. 
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(a) J- J 

n' m' 
— >~ 

R,A 

(b) © = ^ +^ + ... 




Fig. 4. (a) analytic expressions for the lines on the diagrams; (b) diagrams for 
self-energy S. The second term in the self-energy includes an intersection of the 
dashed lines ant it is smaller than the first term by a factor 1/M. (c), (d) diagram- 
matic representation for irreducible averages (25). 

field introduces new energy scale 

h/r^ = AN^/2rr (26) 

If this scale is smaller than the Thouless energy Et, the universal description 
holds, and we find that the system under consideration is at the crossover 
between the Orthogonal and Unitary ensembles. Phenomena associated with 
the energy scale smaller than fi/r^', are described effectively by the unitary 
ensemble, whereas phenomena associated with larger energy scales still can be 
approximately described by the orthogonal ensemble. This conclusion allows 
for a simple semiclassical interpretation of the relation (20) between iV^ and 
the real magnetic field B, as we now explain. 

Consider a classical trajectory of an electron starting from a point r and 
returning to the same point, see Fig. 5. The quantum mechanical amplitude 
for this process contains an oscillating term e^'^"^^^, where Sd is the classical 



_\_ 

S + /0 
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action along this trajectory. In a weak magnetic field the classical trajectory 
does not change, while the classical action acquires an Aharonov - Bohm phase: 

Sci ^ 5ci + — Ai, (27) 
c 



where Ad is the directed area swept by the trajectory, see Fig. 5. If the time it 
takes for the electron to return along the trajectory, t, would be of the order 
of the ergodic time H/Et, the characteristic value of this area would be of 
the order of the geometrical area of the dot \Aci\ — Adot (the meaning of the 
Thouless energy Ex is discussed in previous subsections). If, however, t » 
h/Er, the electron trajectory covers the dot A^^ ~ tEx/h times. Each winding 
through the dot adds a value of the order of Adot to Ad] these contributions, 
however, are of random signs, (an example is shown in Fig. 5). Therefore, the 
total area accumulated scales with -\/A^, 

^ AotV^t - Adot [-j^j ■ (28) 



The orthogonal ensemble is different from the unitary ensemble by the fact 
that in the orthogonal ensemble the quantum mechanical amplitudes corre- 
sponding to a pair of time reversed trajectories have the same phase, while 
these phases are unrelated for the unitary ensemble. For a small magnetic 
field, this means, that the trajectories which acquired an Aharonov-Bohm 
phase smaller than unity still effectively belong to the orthogonal ensemble 
whereas, those which acquire a larger phase are described by the unitary en- 
semble. To obtain the characteristic time scale separating these two regimes, 
we require 



eB\A, 



cl 



h. (29) 



Substituting estimate (28) into Eq. (29), and using $ = -B^ot, one finds 

1 Et r ^ " 



t n V$o 



(30) 



Up to a numerical coefficient, this estimate coincides with the energy scale 
given by Eqs. (26) and (20). The small energy scale physics is governed by 
long trajectories with a typical accumulated fiux larger than the fiux quan- 
tum. For such trajectories the interference between the time reversal paths 
is already destroyed. At larger energy scales, the contributing trajectories are 
shorter, and they do not accumulate enough magnetic field fiux to destroy the 
interference. 
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Fig. 5. Example of the self-returning trajectories. 



2.3 Interaction between electrons: The universal description 

Let us now turn to the discussion of interaction between electrons in a quantum 
dot. In the basis of eigenfunctions 0q, of the free-electron Hamiltonian, the 
two-particle interaction takes the form: 

^int = \Y. ^"/375^1,<xiV'i,a2^7,'^2V'5,ai- (31) 

Hereinafter, we will write explicitly the spin indices a for the fermionic oper- 
ators. The generic matrix element of interaction is: 

'Hap^S = J dfidf2V{fi - f2)0„(fi)0^(f2)0*(f2)05(fi). (32) 

Our goal is to show that the matrix elements of the interaction Hamiltonian 
also have hierarchical structure: Only a few of these elements are large and 
universal, whereas the majority of them are proportional to the inverse di- 
mensionless conductance 1/g, and thus small [67,69,66,68,70]. As the result, 
Hamiltonian (31) will be separated into two pieces: 

U — £7(0) I /y(l/9) /oo\ 
-f^int - -Hint + -Hint ■ l^dj 

The first term here is universal, does not depend on the geometry of the dot, 
and it does not fluctuate from sample to sample for samples differing only by 
realizations of disorder. The second term is small as a power of 1/g, and it 
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fluctuates. This term only weakly affects the low-energy {E ^ Et) properties 
of the system. 

The form of the universal term can be established using the requirement of 
compatibility of this term with the Random Matrix model (13). Since the 
Random Matrix distribution is invariant with respect to an arbitrary rotation 
of the basis, the operator H^^^^ may include only operators invariant under 
such rotations. In the absence of the spin-orbit interaction, there are three 
such operators: 

a,(T a a 



The operator n is the total number of particles, S is the total spin of the dot, 
and operator T corresponds to the interaction in the Cooper channel. (For the 
unitary ensemble the operator T is not allowed.) 

Gauge invariance requires that only the product T^T of the operators T and 
can enter into the Hamiltonian. At the same time, SU (2) symmetry dic- 
tates that the Hamiltonian may depend only on (5^)^, and not on separate 
components of the spin vector. Taking into account that the initial interaction 
Hamiltonian (31) is proportional to tp^, we find for its universal part: 

t = EcU^ + Js (^) ' + Jcf^f. (35) 

The three terms in the Hamiltonian (35) have a different meaning. The first 
two terms represent the dependence of the energy of the system on the total 
number of electrons and total spin, respectively. Because both total charge and 
spin commute with the free-electron Hamiltonian, these two terms do not have 
any dynamics for a closed dot. We will see that the situation will change with 
the opening of contacts to the leads. Finally, the third term corresponds to 
the interaction in the Cooper channel, and it does not commute with the free- 
electron Hamiltonian (1) or (13). This term is renormalized if one considers 
contributions in higher order perturbation theory in the interaction. For an 
attractive interaction, Jc < 0, the renormalization enhances this interaction, 
eventually leading to the superconducting instability. (We will not consider 
the case of an attractive interaction here; for a recent review on the physics of 
small superconducting grains, see Ref. [71]. Effects in larger grains, E^. ^ A, 
are reviewed in Ref. [72].) For the repulsive case, > 0, this term renormalizes 
logarithmically to zero. We will be dealing with the latter case throughout the 
paper. 
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Constants in the Hamiltonian (35) are model-dependent. In the remainder 
of this subsection we will show how to calculate them for some particular 
interactions and discuss the structure of the non- universal part H^^/s) 

At this point it is important to mention that universal Hamiltonian (35) is 
defined within the Hilbert space of one-electron states with energies of the 
order or less than Thouless energy. The matrix elements of this Hamiltonian 
are not just the matrix elements of the interaction potential (32). This poten- 
tial is defined in a much wider energy strip which includes one-electron states 
with energies ~ Ep. Virtual transitions between the low-energy Et) sector 
and these high-energy states renormalize the matrix elements of the universal 
Hamiltonian. It turns out [13] that the third term in the Hamiltonian (35) 
corresponding to the Cooper channel of interaction, can be strongly renor- 
malized by such virtual transitions. In order to avoid this complication, we 
will consider the case of the unitary ensemble in the discussion below, where 
the "bare" matrix elements corresponding to the Cooper channel are already 
suppressed by a weak magnetic field (it is sufficient to thread a flux of the 
order of unit quantum $o through the cross-section of a dot), and relegate 
the corresponding discussion of the interaction in the Cooper channel to the 
Appendix B. Omission of the Cooper channel here is a simphfication that will 
not affect our principal conclusions, as we will consider the case of repulsive 
interaction where the Cooper channel, if it were included, is renormalized to 
zero anyway. 

The statistics of the interaction matrix elements can now be related to the 
properties of the one-electron wave functions. The easiest way to study the 
statistics of the wave- functions 0a (r) is to relate them to the Green func- 
tions (5) and then use the diagrammatic technique for the averaging of their 
products. From Eq. (5) we have 

g^{e; fi, f2) - g'^ie; n, = 27ri ^ Mn)<l>:{r2)5{s - £«). (36) 

a 



At given energy e, at most one eigenstate contributes to the sum in Eq. (36). 
Furthermore, it is known that there is no correlation between the statistics of 
levels and that of wave functions in the lowest order in 1/g, see, e.g., Ref. [41], 
so we can neglect the level correlations and average the (5-function in Eq. (36) 
independently. As a result, we can estimate 

Mri)Mr2ye{6/2-\ea-e\)^—J ciei [e?(5i; n, fa)]. , 5 < A, (37) 

£-5/2 
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where we introduced the notation 



(38) 



Let us first calculate the average of the matrix element (32). Because of the 
randomness of the wave functions, the corresponding product in Eq. (32) does 
not vanish only if its indices are equal pairwise. It is then readily expressed 
with the help of Eq. (37), 



A2 



5/2 5/2 

dei / de2\6a5Sfs^{[g{ea + ei;fi,f4)]_[g{ep + e2;f2,r3)]_) 



47r2 



-5/2 -5/2 



(39) 



In deriving Eq. (39) ior a P we have used the relation 



{Mri)Mr2rMri)Mr2re{6/2 - |£„ - e\)9{S/2 - je^ - e\)) 
= {(l>a{fi)cl>a{f2yMri)Mr2r){0{S/2 - - £1)^(5/2 - \ep - e\)) 

= -^{(Pa{ri)(pa{r2)*(p0{fi)(pfs{f2)*). 

The result (39) can be justified also ior a — P and kplri — r2\ ^ ^ (with kp 
being the Fermi wavevector), see Ref. [68]. The corrections to Eq. (39) are of 
the order or smaller than and will be neglected since we are interested 

in the leading in 1/gf terms. 

In the leading approximation inl/g needed to derive Hi^^, the Green function 
entering into the products in the above formula may be averaged indepen- 
dently. Substituting 

{[Gisa + Si; n, f2)]_) = — {kp\n - rll) ; (40) 

with {■ ■ ■)fs denoting the average over the electron momentum on the Fermi 
surface, into Eq. (39), we find 

iVdf{Mn)M^2)<P;{r3Ws{n))^"'^ = [5a55/37-^14^23 + 5aT<^/35^13^24] , (41) 
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where we introduced the short-hand notation 

J^^. = J^{kF\f,-f,\), (42) 

and Vd is the volume of a (i-dimensional grain. Equation (41) does not depend 
on the disorder strength and corresponds to the universal limit. The coordinate 
dependence in function T indicates simply that all the plane waves that are 
allowed by the conservation of energy are represented in the wave function. 

2.3.1 Universal description for the case of short-range interaction 
We start with the model case of a weak short range interaction 

V{r) = AAVd5(r) (43) 

to illustrate the principle, and then discuss the realistic long range Coulomb 
interaction. In Eq. (43), Vd is the volume of the dot in dimension o? = 3 or its 
area for d = 2, and A 1 is the interaction constant. 

Averaging of the matrix elements (32) with the help of Eq. (41) allows us to 
derive the Hamiltonian from Eq. (31). Indeed, such averaging yields 

Now we substitute Eq. (44) into Eq. (31) and rearrange the summation over 
spin indices with the help of identity 

where a = {a^, a^, a^) and are the Pauli matrices in the spin space. As a 
result, the interaction Hamiltonian takes the universal form (34), (35), with 
the coupling constants 

£;e = |AA; J5 = -AA; = 0. (45) 

The vanishing of J,, is a feature of the unitary ensemble, see also Appendix B. 
We see that for a weak short-range repulsive interacton all the matrix elements 
in the universal part of the Hamiltonian are smaller than the level spacing. 
Later we will see that for the Coulomb interaction this is not the case, and the 
constant is large compared with A. 

We derived the universal Hamiltonian (35) in the leading approximation, 
which corresponds to the limit 1/g — > 0. In this approximation, only the 
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"most diagonal" matrix elements of the interaction Hamiltonian are finite. 
In the first order in there are two types of the corrections H^^^j to 
the universal Hamiltonian. First, the matrix elements of the diagonal part of 
the Hamiltonian, Eq. (35), acquire a correction oc 1/g. This correction ex- 
hibits mesoscopic fluctuations and has a non-zero average value. Second, the 
non-diagonal matrix elements of the interaction Hamiltonian for a mesoscopic 
quantum dot become finite, but with zero average value. 

We start with the evaluation of the average correction (H^^p^j) to the matrix 
elements H^^^g of the Hamiltonian (35), 

{H^afi$) = {Ha,,s) - H%. (46) 



To this end, we have to find the 1/g contribution to the averages of the prod- 
ucts (39) of wave functions. This contribution can be found from the corre- 
sponding irreducible product of the Green functions. Using diagrams for the 
diffusive systems, see Fig. 6, we find 



{[g{e + cu; fi, r2)]_ [G{e- fg, u)]^),, = 2Re(^;«(£ + cu; n, r2)g''{e; rg, r^))-,. 
Here 7„ and /„ are defined in Eq. (11), and normalization condition 

/ drfn{f^fm{f^ = Sum, (48) 



is imposed. Substituting Eq. (47) into Eq. (39), and using the condition 
A ^ 7„ we will obtain the 1/g correction to the average product of the 
wave functions (41): 



,(ri)0,(rl)0;(f3)0Kr4))(^/^) =5«.5/..^i3^244 E 7" ^""^^'lif'] 

^ '^'i 7n>0 (£a - +7n 



where we used the short-hand notation (42). 

We are interested in the interactions between the electrons in states sufficiently 
close to the Fermi surface (within an energy strip ~ Et around the Fermi 
level). It means that the difference of one-electron eigenenergies entering in 
the denominators of Eq. (49) are much smaller than 7„, and can be neglected. 
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Fig. 6. Diagrammatic expansion for the correlation of the Green functions of the 
disordered system {Q^Q-^). 




Fig. 7. Diagrammatic expansion for the fluctuations of the matrix element in the 
disordered systems. 

Substituting Eq. (49) into Eqs. (32) and (43), and assuming \ea — e^l ^ E^, 
we find 

(^a/fi) = ciA- (S^sSpy + S^^Sps) ; Cl = - ^ ^ ~ 1, (50) 

where the dimensionless conductance of the dot g is defined in Eq. (15). Thus, 
we have shown that in the metalhc regime, g ^ 1, the non-universal correc- 
tions to the average interaction matrix elements are small indeed. 

Now we turn to the mesoscopic fluctuations SHajj^ys of the matrix elements 
of the interaction Hamiltonian. With the help of Eqs. (32), (43), (37), and 
of the diagrammatic representation for the disorder average. Fig. 7, one ob- 
tains [67,69,66,68,70]: 

.^^E (g^ (51) 

(This result is for the generic case where all indices a, f3, 7, and S are differ- 
ent. If they are not all different, the fluctuations are bigger by a number of 
permutations which that do not change the matrix element. ) Equations (50) 
and (51) show that the typical deviation of all the matrix elements from their 
universal values (44) are indeed small as 1/g. 

All the above derivations of the (x 1/g corrections were performed for diffusive 
systems. The final estimates in terms of the conductance hold also for ballistic 
quantum dots with chaotic classical dynamics of electrons, see Appendix A. 
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2.3.2 RMT for the intra-dot Coulomb interaction 



So far, we have shown that the matrix elements of the Hamihonian of weak, 
short-range interaction are small compared to A. Therefore any observable 
quantity involving electron states within the energy strip ~ around the 
Fermi level, can be calculated by the lowest-order in A perturbation theory. 
The situation is different for the long-range Coulomb interaction 



where k is the dielectric constant of the medium containing the electrons of 
the dot.Q 

For the simplest geometry, the corresponding matrix elements Ha/s^yS are in- 
versely proportional to the linear size of the dot, rather than to its volume (or 
area in dimension d = 2). li the linear size of the dot exceeds the screening 
radius (or the effective Bohr radius = Kh^/e'^m* in d = 2), the matrix 
elements Hai3^s exceed A, and the lowest-order perturbation theory in the 
interaction Hamiltonian fails. 

In the theory of linear screening, it is well-known how to deal with this diffi- 
culty. When calculating an observable quantity, one should take into account 
the virtual transitions of low-energy electrons into the high-energy states. If 
the gas parameter 



(with vp being the electron Fermi velocity) of the electron system in the dot is 
small, then the random phase approximation (RPA) allows one to adequately 
account for such virtual transitions. Accounting for these transitions yields, 
in general, a retarded electron-electron interaction [18]. However, the char- 
acteristic scale of the corresponding frequency dependence is of the order of 
7i, see Eq. (11). That is why at the energy scale \b\ ^ Et we can consider 
the interaction as instant ant aneous one, and derive the effective interaction 
Hamiltonian acting in this truncated space. 

A modified matrix element in RPA scheme is shown in Fig. 8. It involves sub- 
stitution of the bare potential (52) in Eq. (32) with the renormalized potential 
(^i5 ^2)- This potential is the solution of the equation 



' The constant k does not account for the screening provided by the electrons of 
the dot. This screening will be considered in detail below. 



V{r) 



(52) 




(53) 
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Fig. 8. Matrix element in Random Phase approximation. Polarization operator 11 
involves all the transitions where at least one of the states 7, (5 is outside the energy 
strip of the width of the Thouless energy around the Fermi level. 

Vsc (n, f2) = V{fi-f2)~JJ dfsduV (fi - fa) 11 (fs, u) Vsc (u, , 

(54) 

where the integration over the intermediate coordinates is performed within 
the dot only. The polarization operator 11 should include only the states where 
at least one electron or (hole) is outside the energy strip of width e* ^ Et for 
which we derive the effective Hamiltonian, 

n(ri,r2) = 2 2^ — - — ■ e{-ea,ep). (55) 



Here all energies are measured from the Fermi level, and 9{x) is the step func- 
tion. The factor of two accounts for spin degeneracy. We can compare Eq. (55) 
with the usual frequency dependent operator (cu; fi, that involves all the 
electron states 

^R, ^ ^ ^ 0a(n)0;(ri)0/3(f2)0:(f2) 

[u; ri, r2 = 2 ^ sgne^, 56 



and obtain 

n (fi, f2) = -im / — {u- fi, f2) e{\uj\- s*) . (57) 



On the next step, we replace the polarization operator with its average value [18] 

^„>o + In 
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where 7„ and /„ for a diffusive system are defined in Eq. (11). The pref actor in 
Eq. (58) is nothing but the thermodynamic density of states per unit volume 
(area).Q 

Substituting Eq. (58) into Eq. (57) and taking into account that e* ~ ^71, we 
find 

n(r-'i,f2) = -r|r E /n(n)/n(r2). (59) 



Now we can use the completeness of the solution set of the diffusion equation, 

E fn{n)fn{r2) =S{fi- fa), 



and the explicit form of the zero- mode solution fo{r) = OAotif) / \/Vd for this 
equation, in order to present Eq. (59) in the form: 



n(fi,r2) 



5 in - r,) - 1 



^dot(r-*i)^^dot(f2). (60) 



Here 6'dot(^) = 1, if f belongs to the dot and 6'dot(^) = otherwise. The 
structure of Eq. (60) is easy to understand. The first term in brackets charac- 
terizes local screening in the Thomas-Fermi approximation. The last term in 
brackets subtracts the constant-potential contribution of the zero mode which 
cannot induce electron transitions between the levels and therefore can not be 
screened. 

Substituting the polarization operator (60) in Eq. (54), we find the equation 
for the self-consistent potential in the Thomas-Fermi approximation: 



Vsc{n,r2)=V{n-r2) (61) 

2 /■ r If 

j dr^V ifi - fs) V^c {rs, fa) - — / d^Vsc (fi, f2) 



The use of Thomas- Fermi approximation is justified if the linear size of the 
dot exceeds the screening radius. The integrals here are taken over the volume 
of the dot (or over the corresponding area in the 2d case) . We can rewrite the 
integral equation (61) in the more familiar differential form. 



^ The mesoscopic fluctuations of the polarization operator are related to the vari- 
ation of the density of states in the open sample of the size of the order of \ri — ral, 
corresponding discussion can be found, e.g. in Ref. [73]. The latter fluctuations 
are small as 1/g, and taking them into account here would be the overstepping of 
accuracy of 1/g expansion for the matrix elements. 
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-V^-Fsc(r,fi) 



47re^ 



K 



5{f- fi) - — -6'dot(r)V;c(f,fi) 



^dot(^) / C^^2Hc(^2,ri)6'dot(^2) 



(62) 



supplemented with the requirement that the solution vanishes at r — > cxd. In 
the case of a 2d dot, the operator is still acting in the three-dimensional 
space; the right-hand side of Eq. (62) must be multiplied by 5{z), where co- 
ordinate z is directed normal to the plane of the dot. The second term in the 
right-hand side of Eq. (62) represents the familiar Thomas-Fcrmi screening. 
In the three-dimensional case, this term is oc Hc(r)/''"i>' with the screening 
radius td- 

First, we find an approximate solution of Eq. (62) inside the dot. We consider 
only length scales larger than the screening radius. This enables us to neglect 
the left-hand side of Eq. (62) altogether. In addition, we will neglect the cor- 
rections of the order of the level spacing, which are generated by the last term 
in the right-hand side of Eq. (62). Then the solution takes the form: 



14c(r,fi) = 



— -d(r-ri) + V 



(63) 



where the constant V will be determined later. 

The easiest way to find V is to consider Eq. (62) for r outside the dot (keeping 
ri inside the dot). The right-hand side of Eq. (62) then vanishes, 

V?-14c(r,r'i) = 0. (64) 



The constant V defines for this Laplace equation the boundary condition at 
the surface of the dot, 

V,.{f,n%^s = y (65) 



(the surface is defined unambiguously in the limit of zero screening radius). 
After the solution Kic(^i ^i) of the Laplace equation (64) is found, the constant 
V can be obtained by integration of Eq. (62) along the surface of the dot: 

-> 47rp^ 
dSVV,lr,r\) = . (66) 

s 
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For the two dimensional case, Eq. (66) obviously reduces to 

d , 



/ 



(67) 



where the two-dimensional integration is performed over the area of the dot 
(z^O). 

Now the solution of Eqs. (64) and (65) can be formally found with the help 
of the Green function of the Laplace equation 



(68) 



describing the electrostatic field outside the dot. If the system contains metaUic 
gates, Eq. (68) should be supplied with additional Dirichlet boundary condi- 
tions on the surface Si of those gates 



V{fi,r2) 



0, 



(69) 



where index i enumerates the corresponding gates. 
One immediately finds from Eqs. (65) and (68) 



Kc(r,fi) = Vf dS2VV{f,f2) [1 - ^dot(r)] 
s 



(70) 



for the potential V"sc outside the dot. Substitution of Eq. (70) into Eq. (66) 
yields 



c 



(71) 



for the constant V. The geometrical capacitance of the dot is given by formula 
valid for any shape of the dot 



s s 



(72) 



In a two-dimensional system, Eq. (72) is replaced with 



C 



47r 



dfi / df2 dzidz2 'D{fi,f2) 



(73) 
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where the two-dimensional integration is performed over the area of the dot 
{z = 0), and operator dz is defined in Eq. (67). The concrete value of capac- 
itance C is geometry-dependent. If the size of the dot is characterized by a 
single parameter L, then C is proportional to L with some coefficient, which 
depends on details of the geometry. In the practically important case of a 
gated dot, the capacitance is C ~ L^/And, where is the area of the dot, 
and L is the distance from the gate to the plane of the dot. 

The charging energy, e^/2C, is the dominant energy scale for the matrix ele- 
ments of the interaction Hamiltonian, as we will now show. Corrections to this 
scale appear from terms of the order of A in the effective interaction potential 
V^c- In order to find the matrix elements Va/s-ys, we need to know the screened 
potential V^c(^i, "^2) within the dot. The constant part of this potential ~ e^/C 
contributes only to the "diagonal" matrix elements (32), with a = 5; (3 = 
This part does not contribute to any other matrix element, because of the 
orthogonality of the corresponding one-electron wave functions. 

In order to find the non-diagonal matrix elements, we need to account for 
smaller but coordinate-dependent terms in V"sc(ri, ^2). To accomplish that goal, 
we substitute the potential (63), (70) into the left-hand side of Eq. (62) and 
treat it in the first order of perturbation theory in AC/e^ <S 1. As the result, 
we obtain with the help of Eq. (71) 

V^n, ^2) = ^ + - ^2) + V(n)A + V(r2)A. (74) 

Here, the potential V{r) appears due to the finite size of the dot: a charge 
e "expelled" from the point f = ri 2 in the process of screening, can not be 
pushed away to infinity because of the dot's boundary. In the 3d case, this 
charge forms a thin layer near the boundary of the dot; in the 2d case, it 
creates an inhomogeneous distribution over the whole area of the dot. The 
exphcit form of this potential [in units of level spacing A, cf. Eq. (74)] is: 



V{r) = ^ / S{f- fi)d^iVi f dS2^2V{r,, r^); d = 3; (75) 

s s 

V{r) = / dfi d.d., V{r, fi); d ^ 2, 

where Green function V{f, fi) and the dot capacitance C are given by Eqs. (68) 

and (71), and dz is defined in Eq. (67). The importance of the additional po- 
tential (75) for the structure of the matrix elements of the interaction Hamil- 
tonian was first noticed in Ref . [74] . 
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The substitution of the first term of Eq. (74) into Eq. (32) immediately yields 
the new value of the interaction constant Ef. in the n-dependent part of the 
universal Hamiltonian (35). This constant in fact is the single-electron charging 
energy: 

E. = ^. (76) 

Note, that this value is much larger than the single electron level spacing: 

E 



because of a large factor kpL ^ 1. We intend to show now that this large 
scale appears only in the charge part of the universal Hamiltonian and it 
enters neither the spin channel nor the non-universal part 



(1/9) 



int 



The constant Js in Eq. (35) originates from the two-body interaction term of 
the screened potential. The use of the corresponding (second) term of Eq. (74) 
would yield Js = —A/2, which is an overestimation of the spin-dependent 
interaction term. The correct result for Js in the random phase approximation 
reads:[^ 

, A f ^^1^(1 + 2^); ^ = 3; 

-Js = J drJ^{kF\r\) Vsdr) = — x S / i+y/I 



In ! ; d = 2, 

(77) 

where the gas parameter is given by Eq. (53) and the function was 
introduced in Eq. (40). It is worthwhile to notice that at 1, the RPA result 
(77) gives Js = —A/2, which forbids, in particular, the Stoner instability. 
However, in the regime ~ 1, the RPA scheme becomes non-reliable and the 
constant (77) should be replaced with the corresponding Fermi-liquid constant. 
A universal description still holds in this case provided that the Fermi liquid 
description does not breakdown at length scales smaller than the system size. 



^ The additional smallness in the expression for Js arises because the main 
contribution to this constant comes from the distances of the order of 1/kp; at 
such distances the approximation of the local part of the potential by 5 -function 
in Eq. (74) is no longer valid and one should replace: 

AVd,,^ / dk V(k) 
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We now turn to the discussion of 1/g corrections. The first type of such cor- 
rections resuhs from the substitution of the third and fourth terms of Eq. (74) 
into Eq. (32). One thus finds [74]: 



(1/9) 



SasXl^^ + S^A Kp= I dfV{f)<t^a{r2W0{r^). (78) 



The matrix elements in Eq. (78) are random. Their characteristic values can 
be easily found from Eqs. (74) and (49), with the result 



r 1^ 



(79) 



where boo is a geometry dependent numerical coefficient, 



300 



1 ^ \l drV{r)fn{r) 



7i 



TT 



7n>0 



7n 



(80) 



and the potential V is defined by Eq. (75). The coefficient 6oo is of the or- 
der of unity. According to Eq. (78), the matrix elements \H^^^^^ of the 

Hamiltonian H^^^ with a — S or (3 = 'y are of the order of A/^, contrary to 
the assumptions of Ref. [75]. The rest of these matrix elements are smaller. 

One can use the two-body part of the screened interaction potential (74) to 
evaluate them [67]. This part coincides with the simple model (43), up to the 
constant A, which should be substituted by 1/2. After the substitution, one 
can use the results (50) and (51). 

In all the previous consideration, we assumed the potential on the external 
gates to be fixed, sec Eq. (69). Such an assumption was valid because ad- 
ditional potentials created by those gates were implicitly included into the 
confining potential U{f) from Eq. (1). One, however, may be interested in 
a comparison of the properties of the dot at two different sets of the gate 
voltages. To address such a problem one has to include gates into the electro- 
static problem, i.e., to find the correction dU{r) to the self-consistent confining 
potential due to the variation of the gate potentials. This requires that the 
Thomas - Fermi screening of the charge on the external gates by the electrons 
in the dot has to be considered. The resulting equation for this correction is 
similar to Eq. (62): 



- V'5C/(r) 



yd J 



(81) 
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supplemented with the requirement that the solution vanishes at r — > oo and 
has the fixed value on the surface of each gate, Si, [compare with Eq. (69)]: 

Sm\reSi--^V}'^^ i = l,2,..., (82) 

where indices i enumerate the gates, and V^''^ is the electrostatic potential on 
i-th gate. 

Solving Eqs. (81) - (82) involves the same steps as in derivation of Eq. (74) [the 
only difference is that now there is no 5-function in Eq. (63), the right hand 
side of Eq. (66) vanishes, and the boundary conditions on the gate surfaces are 
given by Eq. (82)]. The solution for the potential outside the dot is obtained 
in terms of the Green function (68) - (69), similarly to Eq. (70): 

SU{r) = 6U (f dS2VV{f, fs) - ^ eVg^''> i dS2VV{f, fs), (83) 



where the zero-mode part of the potential of the dot is given by 

5U^-2E,Af, Af^Y^^i^ eNi^CiV^\ (84) 



The geometrical capacitance of the dot C is defined by Eqs. (72) and (73) and 
the geometrical mutual capacitances between the dot and the gates are 



Si s 



(85) 



for three-dimensional dots. To obtain the result for the two-dimensional case, 
the surface integral over S here should be replaced as in Eq. (67). The resulting 
potential inside the dot is found similarly to Eq. (74) in the form 



5U{r} = -2EM - A 



(86) 



where V{r) is defined in Eq. (75) and potentials V^'^^ (r) are given by 



F®(r) 



Vd 
Vd 



6{f- fi)dSiVi f dS2V2V{n, f2); d = 3; 



Si 



dfi dzd,,V{f,fi)- d = 2. 



(87) 
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We see that the effect of the external gates has the same hierarchical structure 
as the electronic interaction inside the dot: The largest scale [the first 
term in Eq. (86)] corresponds to a simple uniform shift of the potential, while 
the much smaller energy scale A characterizes the changes in the potential 
affecting the shape of the dot. Therefore, the effect of the external gates can 
be separated into a universal part 5C/° and a fluctuating non-universal part 



5Ul^ = -2E,J\fSaf,; (88) 

U^^^'^^-AfX'^,-j:j\fjXi,. (89) 

3 

In Eq. (89), the are random Gaussian variables with second moments 

(Vi/3) = %y' (90) 

where the geometry dependent coefficient 6oo is given by Eq. (80) and all other 
coefficients are given by 



bij^-T.- [ dfV^\r)fn{f^ [ drV^^\r}fn(f^, (91) 

^ 7n>0 •' 



7n>0 T'" 

for i, J = 1, 2 In general all these coefficients are of the order of unity. 



2.3.3 Final form of the effective Hamiltonian 

To summarize this subsection, we note that in the absence of superconduct- 
ing correlations, the universal part of the interaction Hamiltonian, Eq. (35), 
consists of two parts. The dominant part depends on the dot's charge n. The 
corresponding energy scale Ec, Eq. (76), is related to the geometric capaci- 
tance of the dot C, Eq. (72), and exceeds parametrically the level spacing A. 
[There are corrections of order A to Eq. (76) arising from the potentials V in 
Eq. (74) and from the local part of the screened interaction.] The next part 

— * _ 

depends on the total spin S of the dot. The corresponding energy scale Js is 
smaller than A. If the level spacing did not fluctuate, then the smallness of Js 
would automatically mean that the spin of the dot S' = 0orS' = f/2 depend- 
ing on whether the number of electrons in the dot is even or odd, respectively. 
Fluctuations of the level spacing may lead to a violation of the strict period- 
icity in this dependence [76-78] . However, the Stoner criterion [79] guarantees 
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that the spin of the dot is not macroscopically large (i.e.,, does not scale with 
the volume of the dot). The strongest non-universal correction to the Hamilto- 
nian (35) comes from the effect of the boundary of the dot on the interaction 
potential, see Eq. (74), and is of the order of A/y/g, where g is the dimen- 
sionless conductance of the dot (15). In the presence of the gate, the universal 
part of the Hamiltonian (including the RMT representation (13) of its 
free-electron part) and the two largest non-universal corrections H^^^ have the 
following form: 



a,7 



(92) 



(93) 



The last term in Eq. (92) is a c-number and it can be disregarded in all 
subsequent considerations. The fluctuations of the matrix matrix elements 
xj^^p are given by Eq. (90), the elements H^a^^l can be estimated from Eqs. (50) 
and (51) with A = 1/2. Finally we should mention that for a two-dimensional 
dot the relation A. holds only if the distance between the dot and gate 

exceeds the screening radius in the dot. 



2.4 Inclusion of the leads 



So far we considered an isolated dot. We have established the validity of the 
universal RMT description both for the one electron Hamiltonian and for the 
interaction effects. Our goal now is to present a similar description for the 
connection of the dot with the external leads. 

We will assume that the leads contacting the dot are sufficiently clean, and 
neglect the mesoscopic fluctuations of the local density of states in the leads. 
In other words, we neglect closed electron trajectories in the leads that start at 
the contacts.[^ We also neglect the electron-electron interaction in the leads. 
In the case of two-dimensional leads, this approximation is justifled if the 
dimensionless conductance of the leads greatly exceeds the number of channels 
in the point contacts connecting leads to the quantum dot, or, equivalently. 



" Accounting for such trajectories would result in an additive contribution to the 
mesoscopic fluctuations of observable quantities; this contribution is negligibly small 
as long as the sheet conductance of the leads is large. 
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when the resistance of the leads is much smaller than the conductance of the 
point contacts. 

We wish to keep the RMT description of the dot. On the other hand, we saw 
that it is justified only for energy scales much smaller than Thouless energy 
Et = (?A, where g is the dimensionless conductance of the dot. The presence 
of point contacts with a total number of propagating modes Nch broadens 
each level by an amount of the order of A^chA. Hence, the use of the RMT 
description requires the condition 



A^ch < g, (94) 

where the dimensionless conductance of the dot is defined by Eq. (15). Since 
g ^ 1, the condition (94) is met even if more than one channel is open. 

The conductance of a point contact connecting two clean conducting continua 
can be related, by the Landauer formula [80-84] , to a scattering problem for 
electron waves incident on the contact. These waves can be labeled by a contin- 
uous wave number k and by a set of discrete quantum numbers j. The number 
k accounts for the continuous energy spectrum of the incoming waves, and j 
defines their spatial structure in the directions transverse to the direction of 
incidence. Although in principle the number of discrete modes j participat- 
ing in scattering is infinite, the number of relevant modes contributing to the 
conductance is confined to the A'ch modes that are propagating through the 
contact; all other modes are evanescent and hardly contribute to the conduc- 
tance. An adiabatic point contact [5,10,85] is an example adequately described 
by a model having A^ch niodes in a lead. 

Let us now turn to a quantitative description of the theory. The total Hamil- 
tonian of the system is given as a sum of three terms [86] , 



Ht^H + HL + Hld, (95) 

where the Hamiltonian of the closed dot, H, will be taken in the universal 
limit, see Eq. (92), describes the leads, and Hld couples the leads and the 
dot. A separation into three terms, like Eq. (95), is well known in the context 
of the tunneling Hamiltonian formalism [87], and in nuclear physics [88]. Its 
application to point contacts coupling to quantum dots can be found, e.g., in 
Refs. [14,86,89,90]. The Hamiltonian of the leads reads 



iVch 



Hl^Vj.±J -ki;}{k)ij,{k), (96) 
where we have linearized the electron spectrum in the leads and measure all 
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the energies from the Fermi leveL The vector k kp is the deviation of 
longitudinal momentum in a propagating mode from the Fermi wavevector 
kp. [For the sake of simplicity, we assume the Fermi velocity to be the same 
in all the modes; the general case can be reduced to Eq. (96) by a simple 
rescaling.] The Hamiltonian (96) accounts for all the leads attached to the 
dot; the lead index (in case the dot is coupled to more than one lead), the 
transverse mode index, and the spin index have all been combined into the 
single index j, which is summed from 1 to N^h, N^h being the total number 
of propagating channels in all the leads. (In this subsection, we reserve Greek 
and Latin letters for labelling the fermionic states in the dot and in the leads 
respectively.) 

The Hamiltonian Hpc in Eq. (95) describes the coupling of the dot to the 
leads. 



Hld = T.T. Uz [Waj^i^k] + h.c] . (97) 

j=l a=l •' 



Here the coupling constants Waj form a real M x N^h matrix W, and M — oo 
is the size of the random matrix describing the Hamiltonian of the dot. We 
emphasize that the matrix W describes the point contacts, not the dot. That 
is why this matrix is not random: all the randomness is included in the matrix 
n from Eq. (92). 

In the absence of electron-electron interactions in the dot, the Hamiltonian 
(95) of the combined system of dot and leads can be easily diagonalized and 
the one-electron eigenstates can be found. These eigenstates are best described 
by the Nch x Nch scattering matrix S'(e), which relates the amplitudes a°"* of 
out-going waves and a}- of in-going waves in the leads. 



(For a precise definition of the amplitudes a™ and a°"* in terms of the lead 
states ipj{k), and for a derivation of the formulae presented below, we refer 
to appendix C.) The matrix S is unitary, as required by particle conserva- 
tion. It can be expressed in terms of the matrices W and 7i that define the 
non-interacting part of the Hamiltonian and a matrix U that describes the 
boundary condition at the lead-dot interface 

The states Tpj (k) in Eq. (96) represent scattering states in the leads, which are de- 
fined with the help of a suitably chosen boundary condition at the lead-dot interface. 
These boundary conditions may lead to a nontrivial scattering matrix 5*0 = UU"^ 
even in the absence of any lead-dot coupling Hld, see App. C for details. 
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S{e) = U 



1 - 2muW^ {s-n + mvWW"^) 



-1 



W 



(99) 



where V = l/(27rfi?) is the one-dimensional density of states in the leads, the 
M X M matrix Ti. is formed by the elements Tia-y of the Hamiltonian of the 
dot (92), and f/f/^ is the scattering matrix in the absence of the coupling 
between the leads and the dot (i.e., when Hj^d = 0). The scattering matrix 
Sij{e) describes scattering of electrons from one lead to another, as well as 
backscattcring of an electron into the same lead. Equation (99) can be also 
rewritten in terms of the (matrix) Green function of the closed dot 



g{e)^{e-n)-\ (100) 



as 



The matrix S describes both reflection from the point contacts (i.e., scattering 
processes that involve the backscattcring from the contacts only, not the dot), 
and scattering that involves (ergodic) exploration of the dot. Alternatively, 
the backreflection from the point contacts into the leads can be described in 
terms of a reflection matrix Tc, which is related to the matrix W as 



A contact is called ideal if Tc = 0, or, equivalently, W'^W = n^u/MA. For a 
single mode contact, the transparency |tcP is given by |icP = 1 — l^'cP- 

For non-interacting electrons Eqs. (99) - (102) essentially solve the physical 
part of the problem since all the observable quantities are expressed in terms 
of the scattering matrix. We now consider three observables in more detail: 
the two-terminal conductance, the tunneling density of states, and the ground 
state energy of the dot in contact to the leads. 

Two-terminal conductance. The two-terminal conductance is defined for a 
quantum dot that is connected to electron reservoirs via two leads (num- 
bered 1 and 2) with A^i and N2 channels each (where N1 + N2 — A^ch), see Fig. 
9. The electron reservoirs are held at a constant chemical potential fi + eVi 
{i = 1,2). Labeling the group of channels belonging to lead 1 (2) with the 
index 1 < j < A^i (A^i + 1 < j < Ach), the total current / through, say, lead 
no. 1 is given by 
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Fig. 9. Setup for a measurement of the two-terminal conductance. The quantum dot 
(grey) is connected to two leads, numbered 1 and 2, that are connected to electron 
reservoirs. A current / flows through the dot as a function of the voltage difference 
Vi — V2 between the reservoirs. 

I^G{V2-Vi), 

where the conductance G is expressed in terms of the scattering matrix S by 
the Landauer formula [93,80-82] 

fpi^) = 1/(1 + e^''^) being the Fermi distribution function. 

It is convenient to use the unitarity of the scattering matrix and rewrite 
Eq. (103) as 



G 



27ih 



ch 



- de 



^ ) TrASAS^ 

OS J 



(104) 



where the A^ch x -^ch traceless matrix A is defined as 



Aij — 5ij X < 



, i = l,...,A^i, 

, i = Ari + l,...,Arch. 



(105) 



ch 



The advantage of Eq. (104) over the more conventional form (103) of the Lan- 
dauer formula is that it separates the classical conductance {2e'^ /h){NiN2/ Nch) 
and the quantum interference correction, the second term on the r.h.s. of Eq. 
(104). 

Tunneling density of states. An important particular case of Eq. (103) is the 
strongly asymmetric setup, where one of the contacts, say the right one (no. 
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1 in Fig. 9), has only one channel (A^i = 1), with a very small transmission 
amplitude, |tcP = 1 — |?"cP ^ 1- The corresponding coupling matrix W from 
Eq. (97) then acquires the form, see also Eq. (102) 



W, 



6ctj 



tr AM 



+ (1 - 5i,) 



(2) 
aj ■ 



(106) 



Here the channel j = 1 corresponds to the right contact, while the channels 
j = 2, . . . , A'ch correspond to the left contact, and we have chosen the basis 
of the states inside the dot such that the right contact is connected to "site" 

(2) 

a = 1. The matrix W^j describes the coupling to the left contact. Substituting 
Eq. (106) into Eq. (99) and using Eq. (106) one obtains, with the help of 
Eq. (103) 



G = GiAM J de 



' de 



(107) 



where the tunneling conductance of the left contact is given byp^ 



The quantity i^ri^) in Eq. (107) is the tunneling density of states and is given 
by (in the next two formulae we omit the superscript (2) in W^"^^): 



1 



1 



1 



e-n- mvWW^ 



11 



Im 

TT 

Im 

TT 



11 



11 



(108) 



a result which one can obtain also by treating the right contact in the tunneling 
Hamiltonian approximation. Here, we introduced the Green functions for the 
dot connected to the leads, cf. Eq. (100) 



qR{A) _ L ^ IT^vWW^ 



(109) 



The extra factor 2 accounts for spin degeneracy. Strictly speaking, for electrons 
with spin 1/2 one should set A^i = 2 and sum contributions to the tunneling density 
of states for the two spin directions. In the presence of spin-rotation symmetry, 
however, the result is the same as that of Eqs. (108) and (110) below. 
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where the "+" ( "— " ) sign corresponds to the retarded (advanced) Green func- 
tions respectively. 

An equivalent form of Eq. (108) is obtained by noticing from Eq. (100) that 
dGa/s/dHn — SaiSip- One immediately finds from Eq. (108) 

and with the help of Eq. (101) and SS'^ = 1, arrives at [94] 



1 (^S* 

= -^TV5t^. (110) 

This indicates that the tunneling density of states can be related to the para- 
metric derivative of the scattering matrix of the dot without the left (tunnel- 
ing) contact. 

Free energy for non-interacting electrons. The thermodynamic potential Q of 
the dot at chemical potential /j, can be found in terms of the Green functions 
as 



T / dsln 



1 + e 



-(£-M)/T 



where the total density of states {i.e., the local density of states integrated 
over the volume of the dot) u{e) is given by [compare with Eq. (108)] 



1 1 

u(s) = ImTr -. (112) 

Similarly to the derivation of Eq. (110), one can rewrite the latter formula in 
terms of the energy derivative of the scattering matrix as [94] 



1 (9'? 

The average number of particles in the dot {n)q = —dQ/dfi (here (. . .)q indi- 
cates the average over the quantum state without disorder average) can then 
be found as 



(n), = I dsMsMe), (114) 
where /^(e) is the Fermi function. 
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Statistical properties of S. The calculation of the statistical properties of the 
two-terminal conductance, the tunneling density of states, and the ground 
state energy is now reduced to the analysis of the properties of the scattering 
matrix for a chaotic quantum dot, which is a doable, though not straightfor- 
ward task. The available results are collected in an excellent review [14] (see, 
in particular, Ch. 2 of that reference). Here, we mention a few results that are 
needed in Sec. 4. For particles with spin, but in the absence of spin-orbit scat- 
tering, the scattering matrix is block diagonal, 5* = diag (S'°, S°), where each 
block 5*° is a matrix of size = Nchf^, the total number of orbital channels. 
In the presence of spin-orbit scattering such a block structure does not exist, 
and one commonly describes S* as a Ai'^j^-dimensional matrix of quaternions 
(denoted as S°), which are 2x2 matrices with special rules for complex con- 
jugation and transposition [54]. For reflectionless contacts all averages that 
involve 5* (or S^) only vanish. 



{Shji{^i)Si^j^{e2) . . .Si„j„{en)) =0, n = l,2,.... (115) 

For nonideal contacts, the average of S does not vanish, and is given by the 
reflection matrix of the contact, 



{Stdii^i)Si2j2i£2) . ■ ■ Si,^j,Xen)) = n n = l,2,.... (116) 

k=l 

Moments that involve both S and have a rather complicated dependence 
on energy (see, e.g., Ref. [95] for f3 = 1). Here, we mention an approximate 
result for the second moment {Si-^jj^{£i)Si*j^{e2)) in case Vc = 0, 



/QO \co* f^\\_ '^ili2^jii2 + ^)^iij2^jii2 

(^n.. (^1)^.2.2 (^2)) - No^ + 2/P-l + 27c^{e2-e,)/A' ^^^^^ 

where A is the mean level spacing of the dot.[^ The approximation (117) is 
valid for ei — 62 = and for large values of | £2 — £^1 1 ^ A. It is also valid 
for arbitrary |ei — £2! if ^ch ^ 1- This is sufficient to estimate average and 
fluctuation properties of the conductance for temperatures T ^ Nch'A (see 
Sec. 4). 

If Tc is nonzero, one finds 

For the symplectic ensemble, for which /3 = 4, S°j is a quaternion and the left hand 
side of Eq. (117) has to be interpreted as the quaternion modulus of (•S'f^ j^'S'i'Jjj)' 
where S"^j^ is the hermitian conjugate of Sf^j^ . The same interpretation applies to 
Eqs. (118) and (122) below. In the r.h.s. of Eq. (118), the complex conjugate {rc)*^j^ 

should be replaced by the hermitian conjugate (rc)!^,^ for /? = 4. 
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N^^ + 2//3 - 1 - IV TcA + 2m{e2 - Si) /A 
^ (2//3-l)(l-r,rt),,,.,(l-rtrJ,,,, 
A^5, + 2//3 - 1 - Tr rcrl + 27ri (£2 - £i)/A ■ 

This approximate result is for large energy differences or for the case A^ch — 
IVr^rJ > 1. 

In Sec. 4 we use the Fourier transform S{t) of the scattering matrix, 



00 

S{e) = J dtS(t)e''\ (119) 


for e in the upper half of the complex plane. The Fourier transform of the 
hermitian conjugate is defined with negative times and for e in the lower half 
of the complex plane, 



S\e) = j dtS\-t)e-''K (120) 


Statistical averages for the scattering matrix in time-rcprcscntation can be ob- 
tained by Fourier transform of Eqs. (115) - (117), recalling that ensemble and 
energy averages are equivalent. In particular, we find that the scattering from 
a nonideal point contact with energy-independent Tc 7^ is instantaneous. 



S{t) = Tc 5{t) + fluctuating part, (121) 

where the reflection matrix of the point contact Tc is defined in Eq. (C.9), 
and the ensemble average of the "fluctuating part" (and of its integer powers) 
vanishes, see Eq. (116). For the correlator of S{ti) and S\t2) for Vc — 0, we 
find from Eq. (117) 

{SljMi)S°*j^{t2)) = ^ + (2//9 - l)Si,j,Sj,i2] 

X 5{h + i2)e-(^<=^+2//3-i)t,A/(2.)^ ^^22) 

where (3 — 1(2,4) corresponds to the orthogonal (unitary, symplectic) ensem- 
ble. 

Statistical properties of Qg. We will see in Section 4.8, that the tunneling 
conductance in the presence of the interaction can not be expressed in terms 
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of the parametric derivative of S. as it was done in Eq. (110). The transport 
in this case will be related to the statistics of the Green functions for the open 
dot (109), and we give those properties below. It is more convenient to write 
the results in the time domain. For averages which include only retarded or 
only advanced components one finds 
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AM' 



(123) 



which means that the attachment of the leads does not change the average 
level spacing in the dot. For averages involving the retarded and advanced 
components, we have for the case of reflectionless contacts 



0:171 



27r 



a272 



ai72 '^0:271 



+ (2//? - l)(^ai7i<^a27iJ 



M2A 

X 5{h + i2)e-(^ch+2//3-l)tiA/(2.)_ (^24) 



The derivation of Eq. (124) may be performed similar to that of Eq. (25). 

This concludes our brief review of the two-terminal conductance, the tunneling 
density of states, and the ground state energy for a quantum dot in the absence 
of electron-electron interactions. The interaction between electrons leads to the 
Coulomb blockade, and simple formulas as Eq. (103), (110), and (113) cease 
to be valid. With interactions, the answer very substantially depends on the 
conductance of the dot-lead junctions. The consideration of this regime will 
be the subject of the two remaining sections. 



3 Strongly blockaded quantum dots 



In this Section, we discuss the regime of strong Coulomb blockade in quantum 
dots. To allow for this regime, the conductances of the contacts of the dot to the 
leads Gi^2 must be small in units of c^/tt/i. As it was already explained in the 
Introduction, junctions to a quantum dot formed in a two-dimensional electron 
gas of a semiconductor heterostructure are well described by the adiabatic 
point contact model. Normally, a small conductance is realized only in single- 
mode junctions. Since there are two leads connected to the dot, labeled 1 and 
2, we have total number of (orbital) channels A^^j^ = 2, and the transmission 
block of the S-matrix of the point contacts has the form 
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The transmission amplitudes for the separate left and right point contacts ti^2 
are related to the corresponding conductances by Landauer formula (103) 




(^1,2), (^1,2) = 1^1,2 



(126) 



where the extra factor of two comes from the spin degeneracy. We recall that 
the coupling coefficients ti^2 characterize only the point contact and are not 
random, see discussion after Eq. (97); they characterize the average conduc- 
tance of each point contact and not the total conductance of a given sample 
consisting of both the two contacts and the quantum dot. 

Since the conductances Gi and G2 are small, we can express the coupling 
matrix W relating the Hamiltonian of the closed quantum dot to its scattering 
matrix in terms of Gi and G2, see Eq. (102), 



Because the coupling is weak, it is possible to construct a perturbation theory 
of the conductance in terms of this coupling. This perturbation theory differs 
substantially for the peaks and the valleys of the Coulomb blockade. These 
two cases will be considered separately. 

3. 1 Mesoscopic fluctuations of Coulomb blockade peaks 

In the weak tunneling regime, the charge on the dot is well defined — quantum 
fluctuations of the charge are small — , and the particle number n is quantized. 
In the course of an electron tunneling on and off the dot, its charge varies by 
one. The activation energy for the electron transport equals the difference be- 
tween the ground state energies of the Hamiltonian (92) with two subsequent 
values of n. A peak in the conductance as a function of J\f occurs at those 
values of J\f where the two ground states are degenerate. If one neglects the 
level spacing A, the conductance peaks are equally spaced: Mn+i — J^n — 1- 
According to Eq. (92), the presence of a finite (and fluctuating) level spac- 
ing yields a random contribution to the spacings of the peaks. We will not 
discuss the mesoscopic fluctuations of the peak spacings, where a signiflcant 
disagreement between the theory and experiments still exists (see Section V.E 
of the review [15] for a discussion), and concentrate on the amplitudes of the 
conductance peaks. 

To describe the conductance peaks at low temperatures, T -C A, it is sufficient 
to account only for the two charge states mentioned above, in which the dot 



1_ / AM(gi) 
27rV ly 



+ SajSj2 — 




(127) 
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carries n and n + 1 electrons, respectively, and does not have any particle-hole 
excitations. Even under these simplifying circumstances the transport problem 
is not entirely trivial. Every time the equilibrium number of electrons on the 
dot, n, is odd, its state is spin-degenerate, which may result in the Kondo 
effect. The transition between the spin-degenerate state and the singlet state 
occurs just at the gate voltages J\f that correspond to the points of charge 
degeneracy n ^ n + 1. At these points, the quantum dot attached to the leads 
is in a state similar to the mixed valence regime of an atom embedded in a 
metaUic host material. We postpone the discussion of these dehcate many- 
body phenomena till Section 3.3, and confine ourselves here to the results 
of the theory of rate equations. The applicability of this theory sets a lower 
limit for the temperature in the discussion below: T must exceed the width 
{gi +92)^ of a discrete level broadened by the electron escape from the dot. At 
such temperatures, the quantum coherence between the electron states on the 
dot and within the leads is not important. In this section we will also restrict 
ourselves to the case that only a single discrete level in the dot contributes to 
the current, i.e., we restrict ourselves to the temperature interval 

(^1 + 52) A < T < A. (128) 

The case A <^ T <^ Ec is briefly discussed in Subsection 3.4. 

The upper bound in Eq. (128) ensures that only the last occupied discrete 
level in the dot contributes to the electron transport. In the absence of the 
interaction this level can be in four states: we denote by the probability 
of this level to be occupied with one electron of given spin projection a =t 
, I, by Po to be empty, and by P2 to be occupied with two electrons. Those 
probabilities are normalized as Pq + -P2 + -Pj + -Pj. = 1- If the level is adjusted in 
resonance with the Fermi level in the leads, then due to the electron transfer 
through the dot, the state of the level changes between all of those states. 

The picture is changed when the interaction between electrons is included. As 
can be seen from Eq. (93), only three states with the number of electrons on the 
dot differing by 1, can be in resonance. States with a different charge always 
have an activation energy of the order of the charging energy E^.. As long as 
Ec ^ T, these states do not participate in the electronic transport through the 
dot. A conductance peak is characterized by the condition \Af—Af*\^T/Ec <^ 
1. Depending on J\f*, one can distinguish two cases: 

(i) A/"* — 2j + 1/2, with j being an integer, which corresponds to states |0), 
I I), and I I) in resonance, and |2) having an extra energy E^, 

(ii) A/"* = 2j — 1/2, which corresponds to states |2), | f), and | |) in resonance 
and |0) having an extra energy Ec- These limitations on the Hilbert space of 
the level lead to the deviation from the simple Breit-Wigner formula even in 
the temperature regime (128). 
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The form of the conductance peaks can be calculated using rate equations. 
We will present the master equation describing case (i). Case (ii) is consid- 
ered analogously and we will state the result after Eq. (133). Considering the 
Hamiltonian (97) with the coupling matrix W given by Eq. (127), and in the 
Fermi Golden Rule approximation, one finds the following rate equations [96- 
981 



{gi [fiPo - (1 - fi)P^] + 92 [/2P0 - (1 - /2)A]} , (129) 
{g, im - (1 - h)Pi] + 92 [/2P0 - (1 - /2)PJ} , 



dP^ 


A 


dt 


2'Kh 


dPi 


A 


dt 


2Txh 


dPo 


A 


dt 


2Tih 



{^7i[(1-/i)(Pt + A)-2/iPo] 
+ 92 [(1-/2)(Pt + A)-2/2Po]}, 

where gi^2 are the exact conductances of the point contacts, 

g, = {g,)M\i^{l)\\ g2 = {g2)M\i^{2)\\ (130) 



and il){n) are the RMT components of the M-component one-electron wave- 
function of the partially occupied level responsible for the electron transport 
through the dot. It is readily seen that (/i 2 experience strong mesoscopic fluc- 
tuations, which will be important later .Q^ In Eq. (130), the normalization of 
the wavefunction ip is such that (|?/'(n)p) = M~^. 

Among the three equations (129) only two are independent, and we have to 
supplement those equations with the normalization condition Pq + P\^ + P^ = 1 
(recall that P2 = 0, since double occupancy has an energy cost T). In 

the presence of a small bias eV = fii — fj.2 applied between the leads, the 
occupation factors /i and /2 are given by the Fermi distribution functions of 
electrons taken at the energy E±eV/2, where E = 2Ec{N' — N'*) , which is the 
energy it takes to put an electron into the unoccupied state with the lowest 
energy in the dot. 

The current through, say, junction no. 1 is 

^ A(7i[2/iPo-(l-/i)(PT + n)]- (131) 



2nh 



The conductances gi ((72) are conductances for electron transport from reservoir 
1 (2) through the point contact into a dot state with wavefunction ip. They are 
fluctuating quantities, because they depend on the local density of states in the dot. 
In contrast, {gi) and ((72) in Eq. (126) are conductances of the point contacts if 
connected to a (fictitious) electron reservoir on either side. These conductances are 
a property of the contact only, and do not exhibit mesoscopic fluctuations. 
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The stationary solutions for Pq, P| and Pi can be found from the rate equa- 
tions (129). We substitute the result into Eq. (131), take the limit of low bias 
V, and find that the conductance of the system is given by 

r(xr)-^ - gig2 A dfp/dx 



The dependence of the conductance on comes through the energy depen- 
dence of the Fermi distribution function fpi^)'- 

/f(x) = -^, x=^(Ar-Ar*). (133) 
+ 1 1 



Note that the maximum of the conductance is slightly shifted (by ~ T/Ec) 
away from the point A/" = A/"*, and more importantly, the conductance peak is 
not symmetric. This is the consequence of correlations in transport of electrons 
with opposite spins through a single discrete state on the dot. Because of the 
repulsion, no more than one electron can reside in that state at any time (P2 = 
0). Case (ii) is also described by Eq. (132) after the replacement x —x. 

In the maximum, the function (132) takes the value oiP^ 



Equations (132) and (134) relate the conductance through an interacting- 
electrons system to the single-particle electron wave functions; the height of the 
conductance peaks is expressed in terms of gi^2- According to Eq. (130), these 
quantities experience mesoscopic fluctuations, which result in the fluctuations 
of Gpeak- We discuss these fluctuations next. 

Our goal is to compute the distribution function 

W^(G) = (5(G-Gpeak)). (135) 



Most easily this task is accomplished for the pure orthogonal /3 = 1 or uni- 
tary P = 2 ensemble [99]. In these cases, the wave-functions are distributed 
according to Eq. (21). Substituting Eq. (21) into Eq. (130) and the result into 

The numerical prefactor in Eq. (134) is different from that used in Refs. [43,62,99- 
102]. The differences arises, because these references did not properly account for 
the two-fold spin degeneracy in the conditions of the Coulomb blockade. However, 
this does not affect the results for the functional form of the distribution functions 
for the fluctuations of the conductance. 
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Eq. (134), we obtain Gpeak in terms of the average conductances {gi,2) and a 
single random variable a, 

Gpea. = « (3 - 23/^) 4 {^) '^'^^^'^^ (136) 



The distribution function for a is very sensitive to the presence (/? = 2) or 
absence (/3 = 1) of a time-reversal symmetry breaking magnetic field, 



W^=i(a) = ^; (137) 

Wp=2{ci)=a{l - af e-"(i+'^) ji^o [«(! - a)] + ["(1 - a)]} • 

Here -^'0(3^)5 Ki{x) are modified Bessel functions of the second kind [103], 
while the parameter a characterizes the asymmetry of the point contacts: 



(138) 



It is interesting that the distribution function in the orthogonal case (3 = 1 
remains universal no matter how asymmetric contacts are. In the unitary en- 
semble such a universality does not hold. In the case of symmetric contacts, 
{91) = {92), this result was first obtained in [99]. The general case was consid- 
ered in [104]. 

It is important to emphasize that the probability distributions (137) are 
strongly non-Gaussian and they are sensitive to the magnetic field. These 
two features were checked experimentally in [105,21]. The results were in a 
reasonable agreement with the theory,]^ see Fig. 10. 

It is also noteworthy that the average conductance depends on the magnetic 
field: 



(GU) = (3 - 2'/') 4 ( ^) ^-^^ 2F [2-(3+iP-l)a] 



(139) 



The symmetry breaking occurs due to the orbital effect of the magnetic field, and 
one can neglect the role of Zeeman energy. Indeed, in experiments [21] the crossover 
between the ensembles corresponds to a field B ~ lOmT. Zeeman energy at such 
a field would be comparable to temperature only at T ^ 2.5mK, which is 10 times 
lower than the base temperture in the measurements. 
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Fig. 10. Histograms of conductance peak heights for orthogonal (a,c) and unitary 
ensembles (b,d). Data (a,b) and (c,d) are taken from [105] and [21] respectively. An 
example of the raw data used for the graphs (c,d) is shown in Fig. 2 (b). 

where the asymmetry parameter a was defined in Eq. (138) and the dimen- 
sionless function F{x) is given by 

^W-^-^"-"^"!!^)^ ^« = '^ no) = I (140) 



We see from Eqs. (139) and Eq. (140) that the average conductance in the 
magnetic field (unitary case, /3 = 2) is larger than in the absence of the field 
(orthogonal case, (3 — 1). This phenomenon is analogous to the effect of neg- 
ative magnetoresistance due to the weak localization in bulk systems [106] 
with the same physics involved. For the average peak heights, the magnetore- 
sistance is sensitive to the asymmetry of the contacts (138) and in the limiting 
case of a very asymmetric contacts, a — > 1, the magnetoresistance vanishes. 

So far, we have presented results for the peak height distributions for ensembles 
of Coulomb-blockaded quantum dots with either fully preserved or fully bro- 
ken timc-rcvcrsal symmetry. To investigate how the crossover between those 
two ensembles happens for the temperature regime (128), one has to use the 
statistics of the wavefunctions in the crossover ensemble of random Hamiltoni- 
ans distributed according to Eq. (19). This statistics is described by Eqs. (22) 
and Eq. (23). To obtain the moments of the conductance at the crossover 
regime, is a lengthy, albeit straightforward calculation involving Eqs. (130), 
(134), (22) and (23). We refer the readers to the original papers [100,107] and 
review [15] for the results. (The peak height distribution function for the limit 
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a — > 1 corresponding to the case of strongly asymmetric contacts was obtained 
earlier in Ref. [61], see also Ref. [43]. 

Besides the distribution functions of the peaks heights at a fixed value of the 
magnetic field, another quantity of interest is the correlation function of the 
peak heights for different values of the magnetic field. It is defined as 



i{SGiB^y){5G{B,r)f 

SG{B) = G'peak(S) - (G'peak(S).) (141) 

Physically, the correlator C characterizes "how fast" the ensemble changes 
under the effect of the magnetic field. 

The peak-height correlation function (141) was studied numerically in Rcfs. 
[101,102]. For the unitary ensemble, N^^ :$> 1, the results were found to be 
well approximated as 

C(Bi, B2) « n, (142) 

^ ' ' (l + 0.257Vf)2' ^ ' 



where the parameters Nf^ ' are related to the difference and sum of the mag- 
netic fields by Eq. (20). Once again, the characteristic magnetic field is deter- 
mined by the condition Nh — 1. 



3.2 Mesoscopic fluctuations of Coulomb blockade valleys 



In this subsection, we study the mesoscopic fluctuations of the conductance 
in the regime where the gate voltage A/" is tuned away from the degeneracy 
point [108], \Af-Af*\ > T/E^. Equations (132) and (133) predict in this case 
the conductance to be exponentially small G{M) oc e-x.-p[—2Ec\M — M*\/T]. 
However, this expectation is not correct. The rate equations of Sec. 3.1 take 
into account only processes which are of the flrst order in the right or left 
conductance. In these processes a charge is transferred from a lead (initial 
state) to the dot (final state), or vise versa. In the valleys such processes are 
exponentially suppressed because of the high energy in the final state. On the 
other hand, higher order processes allow for both the initial and final states to 
be in the leads, so that the thermal exponent does not appear. Such processes, 
commonly referred to as co-tunneling [28], are suppressed only algebraically, 
G{J\f) oc 1/ [2Ec\J\f — N"*]], and dominate the transport already at moderate 
deviations from the conductance peaks. 
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Fig. 11. Schematic picture of the second order processes contributing to the con- 
ductance in the Coulomb blockade vaheys. Processes (a), (b) are electron- Hke and 
(c), (d) are hole-like. Elastic co-tunneling [28] corresponds to the processes (a) and 
(c), and inelastic co-tunneling [27,28] is described by (b) and (d). 

Such second-order processes are shown in Fig. 11. The processes (a),(b) and 
(c),(d) are electron-like and hole-like processes, respectively. In an electron- 
like process, the intermediate state corresponds to an extra electron placed on 
the dot. In this case the electrostatic part of the energy of the virtual state is 
given by [see Eq. (92)], 



where M* is the half-integer corresponding to the Coulomb Blockade peak. 
For a hole- like process, the intermediate state has one less electron on the dot, 
so that the electrostatic part of the energy of the virtual state is 



In writing Eqs. (143) and (144) we consider the dimensionless gate voltages 
U in the interval U* - I < U < U* . 

Even though all the processes of Fig. 11 are of the same order in the tunneling 
constants, there is an important difference between the transitions (a) and (c) 
on the one hand, and (b) and (d) on the other hand: For the cases (a) and 
(c), the very same electron (hole) that tunnels into the dot leaves the dot, 
whereas in the cases (b) and (d), a different electron (hole) tunnels from the 
dot. The result is that in the cases (b) and (d), the final state differs from 
the initial state by the excitation of an extra electron-hole pair (inelastic co- 
tunneling). The phase volume for such a two-particle excitation is proportional 



(143) 



Eh = 2E,{M -M* + 1). 



(144) 
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to {sehY (if £eh ^ A), where ^ch is the energy of the electron-hole pair. Since 
£eh can not exceed the temperature T, the inelastic contribution vanishes at 
low temperature as T^, even in the limit A — > 0. The elastic contributions (a) 
and (c), in contrast, are temperature independent. For that reason we here 
concentrate on the elastic co-tunneling and return to inelastic co-tunneling 
only at the end of the Section. 

The tunneling conductance can be calculated according to the Golden rule as 

G^2'^u^\A, + Ah\\ (145) 



where u is the density of states per one spin in the leads and the extra factor 
of 2 takes care of the spin degeneracy. The amplitudes Ae and Ah correspond 
to the processes of Fig. 11 (a) and (c), respectively. They can be calculated 
with the help of the Hamiltonian (97), with Eq. (127) for the matrix W, in 
second order perturbation theory in the tunneling probabilities. The result is 



where and ipai^) ^-re the exact eigenenergies and eigenvectors of the nonin- 
teracting Hamiltonian (13). and the argument i = 1, 2 of the wave function ijj 
labels the sites coupled to contacts 1 and 2, respectively, see Eqs. (127). The 
eigenenergies Ea are measured from the upper filled level in the dot, so that 
step function 9{±ea) selects empty (occupied) orbital states for electron (hole) 
like processes of Fig. lla(c). The relative sign difference in the electron, A^, 
and hole. Ah, amplitudes comes from the commutation relation of the Fermion 
operators. 

Substituting Eq. (146) into Eq. (145), we find 

G=i;^{92){9im + Fh\', (147) 



where the dimensionless functions F are given by 



(148) 
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It is also possible, and more convenient for some applications, to express F in 
terms of the exact one electron Green functions of the dot 



where the Green function is given by Eq. (24). 

Equations (147) - (149) solve the part of the problem depending on the in- 
teractions; the conductance is expressed in terms of the single electron eigen- 
f unctions and eigenvalues of the isolated dot. What remains is to perform 
the statistical analysis of the conductance in a fashion similar to what was 
done in Sec. 3.1 for the peak heights. Note, however, that there is an impor- 
tant difference with the case of the peak height statistics: while the height 
of a conductance peak depended on the wavefunction of one level only, many 
levels with energies of the order of £'e(/i) contribute to the tunneling. The su- 
perposition of such a large number of tunneling amplitudes significantly affects 
the conductance fluctuations [108]. Despite these complications, expressions 
for the distribution function of the conductance in the regime of elastic co- 
tunneling, and the correlation function of the conductance in a magnetic field, 
can be found in closed form. We now describe their calculation. 

The cases of the pure orthogonal {f3 = 1) or unitary ensembles {f3 = 2) are the 
easiest to investigate. We use the fact that the eigenvectors and eigenvalues of 
a random matrix in these cases are independent of each other. Moreover, the 
eigenvectors of different eigenstates are also Gaussian variables independent 
of each other, see Eq. (21). 

Using Eq. (148) and Eq. (21), we immediately find 



Since the averaged density of states docs not depend on energy, and since we 
consider Ee, Eh A, we can replace the summation A J2a by the integration 
/ dsa in Eq. (150). The resulting expressions are 




(149) 



— oo 




(150) 
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«f;> = |, (FA) = I - 1 A; 



where all the averages are found in the same fashion as (150). 
Using Eq. (151), we can find the average of the conductance (147): 



(151) 



{G) 



{92) {gi) 



A ^ A 

Ee Eh 



(152) 



This expression was first obtainedp^ [28] by a different method of the en- 
semble averaging. (To calculate the quantitites of the type of Eq. (151), the 
authors of Ref. [28] used the technique developed in the semiclassical theory of 
inhomogeneous superconductors [109,110].) Notice that a magnetic field has 
no effect on the average conductance; (G) remains the same for the unitary 
and orthogonal ensembles. We remind the reader that Eq. (152) is obtained 
within the RMT theory, and therefore assumes Ec <C Et- 

However, as we show below, the fluctuations of the conductance are of the 
order of the average conductance, so that the average (152) is not a good rep- 
resentative of the ensemble. In order to characterize the conductance fluctua- 
tions, we need higher moments of the factors and Fh than those calculated 
in Eq. (148). Using Eq. (21), we find 

{{FeF:r) = {A-f3)A'Y. 0{ea)e{e,) ( ^ ^ ) . (153) 

a,7 \ {£a + J^e) [^-y + ^e) I 



Further simplification follows from the condition i^g ^ A. It enables us to 
neglect the level repulsion and substitute A^ Y^a^ ~^ I dsade^. One eas- 
ily checks that the corrections to this approximation are of the order of 
(A/Ee)^ln(Ee/A) < 1. As a result, we have 

mF:f) = 2{F,F:f + {f:f:){f,f,), (154) 



where the irreducible averages are given in Eq. (151). Equation (154) indicates 
that the amplitudes F entering the conductance are random Gaussian vari- 
ables with zero average and variance (151). (This can be proven by explicit 
consideration of the higher moments.) This observation immediately allows 



Equation (152) is different from Eq. (13) of [108] by a factor of 2, because of 
algebraic error in this reference. 
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one to compute the fluctuations of tlie conductance in tlie valley, 

- ^ (155) 



With the breaking of time reversal symmetry the fluctuations are reduced 
by a factor of 2 similar to the conductance fluctuations in the bulk systems, 
(G) ^ e^/ (nh). For the cases of Ee ^ Eh or E^ ^ Eh (vicinity to the Coulomb 
blockade peak), Eq. (155) was first obtained in [108]. 

It is clear from Eq. (155) that the fluctuations of the conductance are of the 
order of the conductance itself, despite the naive expectation that the large 
number of the contributing states produce the self-averaging. The reason for 
this, is that we have to add amplitudes, not probabilities, to compute the 
conductance .p^ 

Because the fluctuations of the conductance are large, its distribution func- 
tion is not Gaussian. However, it can be found analytically with the help of 
Eq. (147) and the established fact that F's are Gaussian variables with the 
variances (151). Straightforward calculation yields that the distribution of G 
normalized to its average is given by 

\ ^ [e{a)e~^ (3 = 2 



which coincides with Porter-Thomas distribution [59]. This was to be expected 
from the central limit theorem, because the conductance is determined by a 
large number of random amplitudes, which is exactly the assumption behind 
the Porter-Thomas distribution. 

Now, we turn to the study of the effect of the magnetic field on the valley 
conductance. Our purpose is to find the correlation function for the conduc- 
tance fluctuations, similar to Eq. (141), and the distribution function on the 
crossover between orthogonal and the unitary ensembles. Here we discuss only 
the situation in the vicinity of the peaks, which corresponds to the condition 
Ef. ^ Eh or Ef, ^ Ehi so that only one contribution is dominant. The results 
for the general case are presented in Appendix D. Because a large number of 

This statement is valid for the case when the quantum dot is coupled to the 
electron reservoirs via tunneling point contacts, which we consider here. In case of 
a wide junction (i.e., a contact with many contributing channels with small trans- 
parency each), the conductance is a sum over all channels in the contacts of the 
transmission probabilities for these channels. In that case, the relative size of the 
fluctuations is decreased with respect to the average. 



58 



levels contribute to the transport of electrons through the quantum dot, we can 
use Eq. (149) and then apply Eq. (25) to average the product of the Green 
functions. We thus find that the amplitudes F are still Gaussian variables, 
while their averages are modified in comparison with Eq. (151): 



^ ' Ee \27iEeJ' ^ " Ee \27rEj' ^ ' 

(FhF*) = —A , (FhFh) = —A . (158) 

Here the dimensionless function A (a;) is given by 



K{x) = - 



lnxln(l + x^) + Trarctana; + -Li2(— x^) 



(159) 



with Li2(a;) being the second polylogarithm function [103]. The asymptotic 
behavior of function A is k.{x) = 1 + (a;lna;)/7r, for a; <^ 1, and K[x) ~ 
(7ra;)"Mn^ x, for x » 1. The parameters Nh and are defined in Section 2.2. 
The limits of the pure orthogonal (unitary) ensembles (151) are recovered by 
putting A^;P = and = 0(oo) in Eq. (157). 

The correlation function acquires a universal form [108] (ie., all the depen- 
dences for different Ef.[h) can be collapsed to a single curve upon rescaling of 
the magnetic field): 



C(5i,52) 



$1 + $2' 2- 



= A 



+ A 



^1 - $2 



-\\2 



(160) 



where the scahng function A{x) is defined in Eq. (159), $i is the magnetic fiux 
through the dot due to the field Bi, and the correlation magnetic fiux $c is 
controlled by the charging energy 



^ (^) ^ ^ = "^i" ^'^) ' 

Vxg V A y 



(161) 



where $o is the fiux quantum, g is the dimensionless conductance of the closed 
dot, see Sec. 2, and the dimensionless coefficient x is discussed after Eq. (20). 

The dependence of the correlation flux on the charging energy is caused by 
the fact that the transport is contributed to by a large number of states in the 
energy strip of the order of E (see the discussion in Sec. 2.2), which depends 
on the gate voltage. 
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Fig. 12. The correlation function C{AB = B1 — B2) for tlie conductance fluctuations 
in the elastic co-tunneling regime (solid line) and for the peak height fluctuations 
C = [1 + (AS/Sc)^]~^ (dashed line) in the regime $i,$2 3> $c- For the elastic 
co-tunneling, C(A5) is non-analytic at Ai? — > 0, see the inset, and Eq. (160). 

We emphasize that the functional form of C{AB) is different from the results 
for the peak heights fluctuations (142), see Fig. 12. It is worth noticing from 
Eq. (152) and Eq. (161) that in this regime, the correlation magnetic flux 
$c drops with approaching a charge degeneracy point, whereas the quantity 
remains invariant. 

The distribution function of the conductance in the crossover regime can be 
found from Eq. (147) and the fact that the amplitudes F are Gaussian with 
the correlators (157). The calculation yields, instead of (156), a distribution 
function characterized by a single crossover parameter A, 



60 



here lo{x) is the zeroth order modified Bcsscl function of the first kind. In 
the hmiting cases A = 1 and A = 0, the distribution function P{g) coincides 
with the Porter-Thomas distribution (156) for the orthogonal and unitary 
ensembles respectively. 

The dependence of the crossover parameter A on the magnetic flux $ in our 
case is quite nontrivial [108] 



where the correlation flux is given by Eq. (161), and function A{x) is deflned 
in Eq. (159). The crossover in the mcsoscopic fluctuations of conductance in 
the valleys far from the peaks is discussed in Appendix D. 

The interference effects in the Coulomb blockade valleys were studied experi- 
mentaUy in Ref. [111]. In agreement with the theory of this section it was found 
that the average conductance does not show a weak locahzation correction (it 
remains the same for the unitary and orthogonal ensembles). The statistics 
of mesoscopic fluctuations of the valley conductance also agrees reasonably 
well with the theory. Extracted with the help of Eq. (160), the dependence 
of the correlation magnetic fleld on the gate voltage across a valley is shown 
in Fig. 13. One can see that the correlation magnetic fleld increases with the 
deviation from a conductance peak, in agreement with Eq. (161). However, 
the value of the ratio of the correlation flelds for peaks, see Sec. 3.1, and for 
the valleys was found somewhat smaller than the theoretical prediction. 

We recall that a relatively large correlation magnetic field in a valley results 
from a wide energy band (~ Ec/A) of virtual discrete states in the dot. By the 
same token, the conductance in adjacent valleys is also correlated, because all 
but one virtual states participating in transport are the same for such valleys. 
In fact the conductance remains correlated over a large number ~ Ec/A of 
valleys. The corresponding corellation function was calculated in Ref. [112], 
and the valley- valley correlation function for the differential capacitances was 
found in Ref. [113]. 

Considering electron transport through a blockaded dot, we have used so far 
finite-order perturbation theory in the dot-lead tunneling amplitudes. If the 
dot carries a non-zero spin, like in the case of an odd number of electrons on 
the dot, such a treatment may miss specific effects caused by the degeneracy of 
the spin state of an isolated dot. Tunneling results in the exchange interaction 
between the spins of the dot and leads. The exchange, in turn, leads to a many- 
body phenomenon, the Kondo effect. This effect results in an unexpected 
temperature dependence of the conductance across the dot at temperatures 




(163) 



r< A. 
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Fig. 13. Ensemble-averaged characteristic correlation field (solid) and average con- 
ductance (dashed) across peak-valley-peak for ~ 14 independent data sets. 



3.3 Kondo effect in a strongly blockaded dot 



The Kondo effect is one of the most studied and best understood problems 
of many-body physics. Initially, the theory was developed to explain the in- 
crease of resistivity of a bulk metal with magnetic impurities at low temper- 
atures [114]. Soon it was realized that Kondo's mechanism works not only 
for electron scattering, but also for tunneling through barriers with magnetic 
impurities [115-117]. A non-perturbative theory of the Kondo effect has pre- 
dicted that the cross-section of scattering off a magnetic impurity in the bulk 
reaches the unitary limit at zero temperature [118]. Similarly, the tunnel- 
ing cross-section should approach the unitary limit at low temperature and 
bias [119,120] in the Kondo regime. 

The Kondo problem can be discussed in the framework of Anderson's impurity 
model [121]. The three parameters defining this model are: the on-site electron 
repulsion energy U, the one-electron on-site energy sq, and the level width F 
formed by hybridization of the discrete level with the states in the bulk. The 
non-trivial behavior of the conductance occurs if the level is singly occupied, 
r < l^ol < ^, and the temperature T is below the Kondo temperature 

Tk VC/r exp I ^ \ , (164) 
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where £o < is measured from the Fermi level [122]. 



Before considering the Kondo effect in a quantum dot, we first review it for 
tunneling through a single localized level with on-site repulsion U using the 
Hamiltonian of the Anderson model, 



q,(7 a 



(165) 

(166) 



q,(7 



Here a|g^, alq^r, and aj^ arc the electron creation operators in the left and 
right leads (1 and 2), and on the locahzed level, respectively; are the corre- 
sponding energies in the electron continuum. For brevity, the tunneling matrix 
elements ti and t2 connecting the localized state in the dot with the states in 
the leads are taken to be g-indcpcndent. The widths Fj arc related to these 
tunneling matrix elements: Fj = 27rz/i|tip, where Ui is the density of states of 
lead i = 1, 2. At first sight, an Anderson impurity in a two-band model (166) 
may be associated with a two-channel Kondo problem [123]. However, it is 
easy to show that in our case the parameters of this two-channel problem are 
such that it can be reduced to the conventional single-channel one. Indeed, by 
a unitary transformation 



a 



qa 



= uaiqa ± va2qa, with 



u 



qa 



tl\^+\t2\^ I t2 



, (167) 



The Hamiltonian (166) can be converted [120] to the conventional one-band 
model with Anderson impurity [121]. The localized state and band described 
by the fermion operators form the usual Anderson impurity model, which 
is characterized by three parameters: U, Eq, and F = Fi -|- F2. The band 
described by the operators Pq^ is entirely decoupled from the impurity. The 
scattering amplitude .4(1, g 2, g') between two states in the opposite leads 
is related to the scattering amplitudes A{a, q ^ a, q') and A{I3, q ^ fi, q') 
within the bands a and /3, respectively, by relation 

.4(1, g ^ 2, q') = uv[A{a, q ^ a, q') - A{I3, q^/3, q')]. 



The scattering amplitude in the «-band is directly related to the scattering 
phase Sk for the conventional Kondo problem, A{a,q a,q) = exp(2i5;^); 
the scattering problem for the /3-band is trivial, A{P,q 13, q) — 1. Using 
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Eq. (167), we finally find 



1^(1, g ^ 2,g)p = 4|Mt;|2sin2 5 



4^9192 




i9i + 92) 



This way, the "Kondo conductance" Gk associated with the spin-degenerate 
localized level can be expressed in terms of the problem of one-channel scat- 
tering off a single Kondo impurity in a bulk material,P^ 



Here Tk is the characteristic Kondo temperature, and f{x) is a universal 
function. This function, found with the help of numerical renormalization 
group in Ref. [125], is plotted in Fig. 14. A remarkable property of scattering 
on a Kondo impurity, is that the corresponding cross-section approaches the 
unitary limit at low energies, /(O) = 1. The low-temperature correction to the 
unitary limit is proportional to and described by Nozieres' Fermi-liquid 
theory [118], 



So far we considered the Kondo-enhanced tunneling through a single level. 
The same situation can be realized in tunneling through a quantum dot. The 
number of electrons in the dot is controlled by the applied gate voltage. If 
this number is an odd integer, 2n + 1, i.e., the dimensionless gate voltage 
N = CgVg/e lies in the interval 



the Kondo effect leads to a dramatic increase of the conductance at low tem- 
peratures. The consequence is an even-odd alternation of the valley conduc- 
tance through the quantum dot: In "even" valleys, the low temperature con- 
ductance is due to elastic co-tunnehng, while in "odd" valleys, the presence of 

This mapping is possible due to the complete decoupling of mode (3 from the 
Anderson impurity spin. Mapping also works for quantum dots which have spin 
1/2 and are described by the universal model (92). Corrections to this model, see 
Eq. (93), violate the mapping, but do not change qualitatively the behavior of the 
Kondo conductance in the case of S = 1/2. If the spin S > 1/2, then the temperature 
and magnetic field dependence of Gk may be quite different from the predictions 
of Eq. (168), even for a dot described by the universal model, see Ref. [124]. 




(168) 




(169) 



\M-(2n + l)\<- 



(170) 
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Fig. 14. Plot of the universal function y = f{x) vs. x = T/Tk from Ref. ([125]). 

an unpaired spin in the highest occupied level leads to the Kondo effect, which 
causes a remarkable increase of the low temperature conductance compared 
to the even valleys. 

The advantage of using the quantum dots for the experimental studies of the 
Kondo effect stems from an opportunity to effectively control the parameters 
of the system, which is hardly possible for a magnetic impurity embedded in 
a host material. 

Quantitatively, the low-temperature limit for the conductance in the odd val- 
leys is given by Eqs. (168) and (169). The Kondo temperature T^, Eq. (164)) 
depends on the gate voltage through parameter Sq: 

So = 2E^{M - AA*) < 0, U = 2E^, (171) 

where A/"* = n + | is the degeneracy point at which two different charge states 
of the dot have the same energy. 

At a first glance, it appears that the Kondo temperature for a quantum 
dot may be obtained by the substitution of the parameters Eq. (171) into 
Eq. (164). However, it is not exactly the case, as we explain below. Unlike the 
single-level Anderson impurity model, the discrete energy spectrum of a dot 
is dense, A <^ Ec- Still, if the junctions' conductances are small, Gj <^ e^/^, 
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2n+\ ln+1 J\f 




(b) 



Tf, 



Fig. 15. (a) Linear conductance at T » Tk (solid line) and Tk (dotted line). 
Due to the Kondo effect, the function G{M) develops plateaux in place of the "odd" 
valleys of the Coulomb blockade, (b) Sketch of the temperature dependence of the 
linear conductance, G{T), in the "odd" valley. The conductance decreases with the 
decrease of temperature from T ^ Ec down to T ~ T^i ~ \/EcA, see, e.g., [28,108], 
and Eq. (181). At very low temperatures, T ^ Tk, the conductance grows again. 
If the singly-occupied level, which gives rise to the Kondo effect, has equal partial 
widths, Fi = F2, the conductance reaches the unitary limit G = e^/nh at low 
temperatures. 

(strong Coulomb blockade), those levels of the dot which are doubly-filled or 
empty are important only at scales larger than one-electron level spacing, and 
only the one single-occupied level contributes to the dynamics of the system 
at energy scale smaller than A. As a result, the model of the dot attached to 
two leads can be truncated to the Anderson impurity model, however the high 
energy cut-off is determined by A rather than by the charging energy. (See 
Sec. 4.3 for more discussions). This change in the energy cut-off replaces the 
first factor U in Eq. (164) by A. The level width F can be related to the sum 
91 + 92 for a given level in the quantum dot. Finally, Eq. (171) establishes the 
relation between Eq and U on the one hand, and the charging energy Ec and 
gate voltage Af on the other hand. The resulting Kondo temperature (for the 
interval < A/" - M* ~ 1/2) is found as 



Tk ~ AW(^i + ^2)^exp 



271 E, 
91+92 ^ 



(172) 



Equations (168) and (169) tell us that upon sufficiently deep cooling, the 
gate voltage dependence of the conductance through a quantum dot should 
exhibit a drastic change. Instead of the "odd" valleys, which correspond to 
the intervals of gate voltage (170), plateaux in the function G{M) develop, 
see Fig. 15a. In other words, the temperature dependence of the conductance 
in the "odd" valleys should be very different from the one in the "even" val- 
leys. The conductance, determined by the activation and inelastic co-tunneling 
mechanisms, decreases monotonously with the decrease of temperature, and 
then saturates at the value of Ggi in the even valleys, see the previous Section. 
On the contrary, if the gate voltage is tuned to one of the intervals (170), the 
G{T) dependence is non-monotonous, see Fig. 15b. When temperature is low- 
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ered, the conductance first drops due to the Coulomb blockade phenomenon. 
A further decrease of temperature results in the increase of the conductance 
at T ^ Tk- Its T = saturation value depends on the ratio of the partial level 
widths Fi and V2 of a particular discrete level. The partial widths are related 
to the tunneling matrix elements, Fj = 27ri/|tj|^ oc gi, where the matrix ele- 
ments tj, in turn, arc proportional to the values of the electron eigenfunctions 
within the dot, ti oc ipnii), where n labels the eigenstate, and the argument i 
of the wavcfunction tpn refers to the site in the dot adjacent to junction i, see 
Eqs. (126) and (130). In a disordered or chaotic dot, the electron wave func- 
tions are random quantities, described by the Porter- Thomas distribution [59]. 
As we already discussed in the previous section the randomness of the wave 
functions results in mesoscopic fluctuations of the elastic co-tunneling con- 
ductance [108]. The T = value of the conductance in the Kondo regime 
should fluctuate as well, see Eq. (168). Experiments are performed, however, 
at a flnite temperature T, where fluctuations of conductance occur mostly due 
to the fluctuation of the Kondo temperature (172), which has an exponential 
dependence on the level width Fi -|- F2. 

The experimental search for a tunable Kondo effect brought positive results 
[126-128] only recently. In a number of experiments the difference in the low- 
temperature conductance behavior in the "odd" and "even" valleys was clearly 
observed, see Fig. 16. In retrospect it is clear, why such experiments were hard 
to perform. The negative exponent in Eq. (172) contains a product of two large 
parameters, Ec/A and e^/h{Gi + G2), leading to a very small value of Tk- 

To bring Tk within the reach of a modern low-temperature experiment, one 
may try smaller quantum dots in order to decrease Ec/A; this route obvi- 
ously has technological limitations. (In the experiments [126] this ratio was 
Ec/ A^7.) Another, complementary option is to increase the junction conduc- 
tances, so that Gi^2 come close to e^/irh, which is the maximal conductance of 
a single-mode quantum point contact. The junctions in the experiments [126] 
were tuned to G ~ (0.3 — 0.5)e^/7r^. Under these conditions, one would expect 
Tk ~ lOOmK only ii Af - Af* ^ 0.07, see Eq. (172); in the middle of the valley 
{J\f — J\f* = 0.5) the Kondo temperature is unobservably small. Indeed, an 
evidence for the Kondo effect was found [126] in a relatively narrow interval 
interval, jV — J\f* ^ 0.15. Only in this domain of gate voltages the anomalous 
increase of conductance G(T) with lowering the temperature T was clearly 
observed; the anomalous temperature dependence of the conductance was not 
seen in the middle of the Coulomb blockade valley. 

In conventional Kondo systems, a signature of the effect is in the characteristic 
minimum of the resistivity. The minimum comes from the competition of the 
two contributions to the resistivity: the phonon contribution is decreasing 
with lowering the temperature, while the Kondo contribution has the opposite 
temperature dependence. A similar feature in tunneling through a quantum 
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Fig. 16. The gate voltage dependence of the conductance at various temperatures 
(A), and the temperature dependence of the conductance in the valleys labelled 
4 and 5 (B), see [127]. Reprinted with permission from Science 281, 540 (1998). 
Copyright 1998 American Association for the Advancement of Science. 

dot would be a minimum in the conductance in an "odd" valley. However, such 
a minimum is quite shallow at Gi,G2 /%. Indeed, the minimum comes 
from the competition of the temperature dependencies of the inelastic co- 
tunneling and of the Kondo contribution. In the domain T <S A the inelastic 
contribution to the co-tunncling is exponentially small. Gin oc cxp(— A/T). 
The proper expansion of the function /{T/Tk) in Eq. (168) at T ^ and 
the use of (172) yield the following temperature dependence of the Kondo 
contribution to the valley conductance 



(G(T)) = (Gel) 



1 + 



h{Gi + G2) 



(173) 



where (Gel) = {h{Gi) {G2) /4:7re'^){A/ Ec) is the average elastic co-tunneling 
conductance, see Eq. (152). As one can see from Eq. (173), the Kondo correc- 
tion to the conductance remains particularly small everywhere in the temper- 
ature region T ^ Tk- 

In order to increase the Kondo temperature and to observe the anomaly of 
the temperature-dependent conductance G(T) under these unfavorable condi- 
tions, one may try to make the junction conductances larger. However, if G12 
come close to e^/vr^, the discreteness of the number of electrons on the dot is 



68 



almost completely washed out [39]. Exercising this option, therefore, raises a 
question about the nature of the Kondo effect in the absence of strong charge 
quantization. We address this question later on, in Section 4.3. 



3.4 Overall temperature and gate voltage dependence of the conductance 



Closing this Section, we briefly discuss the gate voltage dependence of the 
linear conductance through a strongly blockaded dot in various temperature 
intervals. 

As we have already mentioned in the very beginning of Subsection 3.2, rate 
equations predict that the peaks in the G{J\f) function should have width ~ 
T/Ec] the conductance in the valleys is exponentially small, with the activation 
energy set by the charging effect. In the temperature interval Ec ^ T ^ A, 
the rate equations yield [129] for the peak shape 

- 2^^sinh[£;,(AA-AA*)/r]' ^^^^^ 



where A/"* is a half-integer number, and the parameter g is determined by 
the conductances of the point contacts {gi,2), the statistical properties of the 
wave-functions of the dot, and by the rate of inelastic processes. 

If the rate of inelastic processes is small compared to the inverse electron 
dwell time in the dot {gi + g2) A, then the parameter g is given by 



Here, /3 = 1, 2 in the limits of no and strong magnetic fields respectively; the 
parameter of asymmetry of the point contacts a and the function F are defined 
in Eqs. (138) and (139). The conductance (174) has the same dependence 
on magnetic field as the averaged peak conductance in the low-temperature 
(T <S A) regime, see Eq. (139). Notice, that unhke the low-temperature peak 
conductance, the mesoscopic fiuctuations on top of result (174) are suppressed 
by a large factor of A/T due to the large number of discrete levels contributing 
to the transport; all the peaks have the same height, which equals the half of 
the high-temperature (T 3> E^.) value of conductance through the dot. This 
is a very straightforward consequence of the Coulomb blockade: out of the 
two elementary processes of charge transfer (one is increasing, and another is 
decreasing the initial charge of the dot by e), only one process is allowed. 

The important difference between the case of sequential tunneling through 
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the quantum dot mentioned above and the process of elastic co-tunnehng 
considered in Section 3.2, is that for sequential tunneling each level contributes 
to a real process, rather than serves as a virtual state for tunneling as is the 
case for elastic cotunneling. Therefore, for sequential tunneling not amplitudes, 
but the probabilities of tunneling through each level are added. 

In the case of rapid inelastic relaxation, l/T^p^ {gi + (72) A, the time elapsing 
between the events of an electron tunneling into and out of the dot is suffi- 
ciently long that an equilibrium electron distribution can be established inside 
the dot. If, in addition, the temperature T ^ A, an electron entering the dot 
uses one of many (~ A/T) discrete levels in it. Similarly, an equally large 
number of levels participates when the electron exits the dot. Because the en- 
ergy level used for the tunneling from the dot is most probably different from 
the one used for the tunneling into the dot, the dimensionless conductances 
gi and g2 defined in Eq. (130), has to be independently averaged over many 
levels, yieldingP^ 

— (176) 



{91) + {92)' 



The two cases of rapid and slow inelastic relaxation can be distinguished 
through their magnetoconductance, which, if normalized by the average peak 
height, is temperature independent as long as inelastic processes are not slow, 
I/t^ ^ {gi + (72)A, but quickly disappears for A when ^ {gi + g2)A 
[130]. Recent measurements of the magnetoconductance of Coulomb blockade 
peak heights show values between the two extremes of rapid and slow inelastic 
relaxation [131]. 

At lower temperatures T ^ A, the peaks heights are proportional to 1/T, and 
exhibit strong mesoscopic fiuctuations, see Subsection 3.1. The peak height 
reaches 

'^P-k ~ Itt-^tZ^, (177) 



when T becomes of the order of the discrete level width [gi + g2)A. At the 



'^^ Higher-order in gi and 52 corrections to Eq. (176) were considered in Refs. [30,40] 
for the cases of single-channel and multichannel junctions respectively. The latter 
case (51,52 <C 1 at Nrmch) is typical for the metallic Coulomb blockade systems. 
Summation of the leading series in (51 -|- 52) ^^(Ec/T), performed in the limit Nch — > 
00 results [30] in the replacement g 2g/[2+{{gi) + {g2)) ln(£'c/r)]. It is noteworthy 
that the increase of the contacts conductances result in the decrease of the peak 
conductance G through the device. 
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same time the dimensionless peak width (57Vpcak decreases down to the valueP^ 

5AApeak~ ('7l + (72)^. (178) 

When the temperature is further lowered, the peak continues to grow, and si- 
multaneously broadens again and becomes asymmetric. However, these effects 
become noticeable at much lower temperature, T ~ Tk, when the Kondo effect 
starts playing role. 

In the valleys, upon lowering the temperature, the exponential decrease of the 
conductance that was mentioned above is replaced by a weaker, power-law 
(cx T^) temperature dependence provided by the mechanism of inelastic co- 
tunneling [27,28]. In the middle of the valley this mechanism results in the 
conductance 




(179) 



One can see that the crossover between the two regimes occurs at temperatures 

Tc^T,^ = EJ\\n{{g^) + {g2))\. (180) 

A comparison of Eq. (180) with the result (152) for elastic co-tunneling shows 
that the crossover to the temperature-independent conductance occurs at 

T ~ Tel = ^/EA. (181) 

If the dot carries a finite spin, then at much lower temperatures, T ^ T^, the 
conductance increases again due to the Kondo effect, see Subsection 3.3 and 
Fig 15. 

In a typical experiment with a quantum dot formed in a semiconductor het- 
erostructure, the five relevant energy scales are related as: 

T/c i$ A i$ Tq\ ^ Tin ~ Ec- (182) 

Temperatures T easily accessible; at such temperatures the con- 

ductance fluctuations are of the order of the average conductance. At higher 
temperatures, when the main conduction mechanism switches to the inelas- 
tic co-tunneling, fluctuations are still determined by the elastic mechanism as 

Width of the peak at larger values of conductances is discussed in Ref. [132] 
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long asT < (E^Ay. At even higher temperatures, the inelastic contribution 
dominates also in the fluctuations. Their relative magnitude, however, be- 
comes small, {SGf^) / {Gi-a)"^ ~ A/T. Simultaneously, the correlation magnetic 
field becomes independent of the gate voltage, as in the inelastic co-tunneling 
regime it is controlled by the temperature rather than by the charging energy. 



4 Weakly blockaded dots 

Charge quantization deteriorates gradually with the increase of the conduc- 
tance of the junctions connecting the quantum dot to the leads. The charac- 
teristic resistance at which a substantial deterioration occurs can be estimated 
from the following heuristic argument. The energy of a state with one extra 
electron in the dot is Ec. Dispensing with the discreteness of the quasiparti- 
cle spectrum (A —>■ 0), one can estimate the lifetime of a state with a given 
charge as the i?C-constant of a circuit consisting of a capacitor with capaci- 
tance C = e'^/2Ec, imitating the dot, and a resistor with resistance R= 1/G, 
imitating the dot-lead contact. According to the Heisenberg uncertainty prin- 
ciple, this state is well-defined as long as -Ec~ TiG/C . Therefore charge quanti- 
zation requires a small lead-dot conductance, G ~ /%. We refer to dots with 
G^e^ /h, for which the charge on the dot is quantized, as "strongly blockaded" 
dots. Dots with G ~ e^//i, for which charge is not quantized, are referred to 
as "weakly blockaded" . While the constraint on G involves only the universal 
quantity e^/^, the crossover from "strong" to "weak" Coulomb blockade oc- 
curring at G ~ /h is not universal. The detailed behavior of a partially open 
dot depends on the number of modes propagating through the junction and on 
their transparency, even in the limit A ^ 0. In the case of a tunnel junction 
with large number of modes each of which is characterized by a small trans- 
mission coefficient, charging effects vanish gradually at G~ /ith. At large G, 
the remaining periodic oscillations of the average charge of the dot with the 
varying gate voltage were shown to be exponentially small, proportional to 
exp(— vr/iG/e^) [133] (a controllable calculation of the pre-exponential factor 
still remains an open problem). 

In the case of a single-mode junction. Coulomb blockade oscillations of the 
average charge vanish exactly at G = e^/vr/;,, see Ref. [38,39]. Here we concen- 
trate on the case of a narrow channef^ supporting a few propagating modes 
and connected adiabatically [10] to the lead and quantum dot at its ends. 
This is the conventional geometry of a quantum dot device formed in a two- 
dimensional gas of a semiconductor heterostructure [8]. The theory for such 



We use the word "channel" instead of "contact" to emphasize that the structure 
of the wavefunction in and near the contact is one-dimensional, see Sec. 4.2. 



72 



a model is quite well developed, but the calculations are very involved.P^ In 
this review, our aim is twofold. On the one hand, we want to present the de- 
tailed derivations that transform the problem into a one that can be solved 
using well-developed formalisms. On the other hand, we also want to give the 
general picture of the thermodynamic and transport properties of a weakly 
blockaded dot, as it emerges from the theory. Therefore this Section is orga- 
nized as follows. We start with Subsection 4.1, where an overview of the final 
results is given. For those readers who wish to enter into the more detailed 
derivations, this Subsection may also serve as a guide to the more detailed 
presentations, which can be found in Sections 4.5-4.8 and Appendices F-H, 
and use a combination of the bosonization and perturbation theory methods. 
Before presenting a rigorous treatment of the mesoscopic fluctuations, we de- 
scribe a simplified version of the theory in Sections 4.2-4.4. This version is 
still sufficient for performing the estimates of the most important observables. 
Going beyond the estimates, however, requires a rigorous theory. It is neces- 
sary not only for fixing the numerical values of the coefficients in formulas 
for various measurable characteristics, but also for finding their dependence 
on control parameters. One can view the result for the correlation function 
of mesoscopic conductance fiuctuations (367) as an example. Formalism and 
theory presented in Sections 4.5-4.8 were needed to predict the four- fold re- 
duction of this function in the magnetic field [70]; this prediction was later 
checked experimentally [154]. 

4-1 Main results for the conductance and differential capacitance of an almost 
open quantum dot 

Conductance through a quantum dot. 

We first discuss the case of a dot connected to two leads by almost refiectionless 
single-mode junctions. We assume that the refiection amplitudes vi and r2 in 
the two contacts connecting the dot to the leads are of the same order, 

ri ^ r2 = r, (183) 

and obey the condition 

E^\r\^<A. (184) 

For such an open dot the conductance G depends only weakly on temperature 
and gate voltage as long as T ^ A. In the regime ^ T A, the disorder 

Coulomb blockade at G ~ e'^/nh in the case of large number of channels, A'ch ^ 1, 
in our opinion, is less understood. General scaling ideas regarding the evolution of 
the Coulomb blockade with the increase of G in this case are presented in Ref. [134] 
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averaged conductance (G) is 



(G) 




(185) 



where C 0.577 is the Euler constant, and T{x) is the Gamma function. This 
result summarizes Eqs. (346) and (349) below, see also Ref. [40]. The symmetry 
index ^ — 1 (2, 4) for the orthogonal (unitary, symplectic) ensemble. The two 
temperature-independent terms in Eq. (185) correspond to the classical resis- 
tance of the quantum dot and the interference (weak-localization) correction 
to it, and are known in the context of mesoscopic systems without electronic 
interactions. The temperature-dependent terms appear as a result of electron- 
electron interactions. However, they are small compared to the temperature- 
independent terms at temperatures A. As a matter of fact, the amplitude 
of "usual" mesoscopic fluctuations characteristic of a non-interacting system 
exceeds the interaction correction to the average conductance, and dominates 
the variance and the conductance correlation function {6G{J\fi)SG{Af2)), as 
can be seen from Eqs. (349) and (351) below [jV is the dimensionless gate 
voltage, see Eq. (92)]. 

The temperature dependence of the interaction correction to the conductance 
saturates at T ~ A, when its magnitude is comparable to the weak-localization 
correction. At these low temperatures, interaction effects become important. 
Under the condition Eq. (184), Coulomb blockade does not manifest itself 
in the form of oscillations of the ensemble averaged conductance {G{J\f)) as 
a function of the dimensionless gate voltage jV. Instead, Coulomb blockade 
affects the mesoscopic fluctuations of G, as described by the conductance cor- 
relation function {SG{J\fi)SG{J\f2)), which acquires an oscillatory dependence 
on A/i — A/2. The amplitude of oscillations is of the order of e'^/nh. The os- 
cillations remain correlated over an interval — A/2I ~ E^/A. Analytical 
formulae for this regime are not available yet. 

In the case of stronger reflection. 



the above consideration and Eq. (185) is applicable, though in a narrower 
temperature interval. 



Ec\r\^ > A, 



(186) 



(187) 
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At lower temperatures and for contacts with equal reflection probability, |ri| = 
\r2\, the conductance shows peaks with maxima of order G ~ e'^/nh, see 
Ref. [40]. The maxima are well pronounced if Eq. (186) is satisfled: in the 
valleys of Coulomb blockade the estimate for the conductance is 

A 

Gmin ~ ^ I (188) 



At even lower temperatures, 



rn <. T T 

1 ~ Iki J-k 



Er\r\ 



exp 



-a (AT) 



EM 



;i89) 



the Kondo effect develops, and the conductance in the "odd" valleys ap- 
proaches the unitary limit. The function a (A/") is positive with values ~ 1; 
the pre-exponential numerical factor is ~ 1. While the estimates Eqs. (188) 
and (189) are easily obtainable by the combination of methods reviewed in 
Sections 3.2 and 4.3, the explicit dependence of the conductance on gate volt- 
age across a valley and the explicit form of the functional dependence on M 
in the exponential of Eq. (189) are not known. 

If the reflections in the junctions are of different strengths, then several tem- 
perature intervals appear, and the behavior of the conductance is signiflcantly 
different in each interval. We will review here only the case of a strongly asym- 
metric setup: in one channel, backscattering is strong, Gi = 1 — |rip ^ 1, 
while in the other channel, backscattering is weakj^f] but still not inflnitesi- 
mally small: 

A/E^ < |rp < 1. (190) 



(Here we wrote r for the backscattering amplidude r2 in the weakly backscat- 
tering channel.) In this conflguration, the flrst contact can be treated within 
the tunneling Hamiltonian formalism, and all the results of Sec. 4.8 are appli- 
cable. We flnd from Eqs. (358) and (359) of that section 



(G(T)) ^ Gi 



TT 



3T 



['ie^EMMy 



T->E,- 
E,->T->T^{J\fy 

T<To(Ar), 



(191) 



24 



The case of |r2| ^ |ri| ^ 1 is considered in Ref. [135]. 
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where To(7V) is the energy scale at which crossover to the Fermi hquid behav- 
ior, G oc occurs [cf. Eq. (221)], 

To{Af) = —E,\r\^ cos' ttAT. (192) 

TT 

The transport processes contributing to Eq. (191) are inelastic. At tempera- 
tures T <C To(A/') the inelastic mechanism results in the Coulomb blockade 
oscillations. However, this contribution to the conductance vanishes at T — 0. 
At very low temperatures the elastic co-tunneling becomes dominant in the 
electron transport. This mechanism yields, see Eq. (366), 



(Gel) 



Gi In 

2Er 



E, 



To (AT) 



T« JATo(Ar). 



(193) 



Here, in accordance with the left inequality of Eq. (190), we assumed To{Af ) 3> 
A. In this temperature regime, the mesoscopic fluctuations of the conductance 
are of the order of conductance itself. 



The result (193) does not depend on the temperature and one might think 
that this is the limiting low temperature behavior of the conductance. This 
is indeed true if the average charge of the dot (Q) = eAf, with 2j — 1/2 < 
A/" < 2j + 1/2 ("even valley"), where j is an integer number. In this case, the 
ground state of the closed dot is non-degenerate and there is no new physics 
at low temperatures. 

If however 2j + 1/2 < J\f < 2j + 3/2, ("odd valley"), the conductance at 
low temperatures shows the re-entrance behavior due to the Kondo effect, see 
Eq. (240), 



G(T) = Gi 



■ Ae-c, 

ai In 

2E, 



' E, \ 32 To (AT) / T ' 



. (194) 



Here, the first term represents elastic cotunneling, cf. Eq. (193), while the 
second term is due to the Kondo effect. Further, / is the universal scaling 
function [125], see Fig. 14, and the Kondo temperature is given by 



The parameters ai, a2 and as are independent random quantities!^ obeying 

note, that in the case of asymmetric setup, the AA-dependence of the exponential 
in Eq. (195) is known, unlike the case |rip = 1^2^, cf. Eq.(189). 
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Porter- Thomas statistics [59], see Eq. (156). Therefore, the Kondo contribu- 
tion also experiences mesoscopic fluctuations that are of the same order as 
the average conductance. At T <C Tk the Kondo contribution exceeds the 
elastic co-tunneling contribution to the conductance by the large parameter 
To{^f)/A > 1. 

The above formulas for the co-tunneling contribution to the conductance, sec 
E's. (193) and (194), are invalid in narrow intervals around the conductance 
peaks, where To{M) becomes smaller than A. In these intervals, the factor 
To{Af) in the argument of logarithm in Eqs. (193) and (194) should be replaced 
by A. The schematic temperature dependence for the conductance in even and 
odd valleys is similar to the one shown in Fig. 15. 

Differential capacitance of an open dot. 

In addition to complications stemming from quantum fluctuations of charge, 
the transport phenomena discussed here involve one more important part of 
many-body physics: the Anderson orthogonality catastrophe [136]. This ele- 
ment becomes especially important in the analysis of conductance in a strongly 
asymmetric setup, see the discussion following Eqs. (293) and (356). Thermo- 
dynamic quantities are free from that complication. Therefore measurement of 
a differential capacitance provides information additional to that found from 
the transport measurements. Experimental investigation of the capacitance of 
a quantum dot is quite a challenging task [137], and by now there were no 
systematic study of mesoscopic fluctuations of this quantity. 

In presenting the theoretical results, we concentrate on the case of a dot con- 
nected to a single lead. We also assume that the junction allows for propaga- 
tion of a single mode of electrons carrying both directions of spin (small or 
no Zeeman splitting) . Reflection of the mode is characterized by the reflection 
amplitude r. 

At a sufficiently high temperature, 

T:$> Ec\r\^,A, (196) 

the result for the average capacitance can be found in Ref . [39] , 

~ ln{EjT)\r\'cos{27rAf). (197) 



The derivation of that result and its generalization to the case of many modes 
(see also Ref. [138]) is given in Section 4.6, see Eq. (323) and the discussion 
around it. 

Upon lowering the temperature it becomes important, whether the energy 
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scale £^c|^P is larger or smaller than the mean level spacing A. If the condition 

To (AT) > A (198) 

with To (A/") from Eq. (192) is satisfied, then oscillations of the differential 
capacitance follow the law [39] 

^^^^ - In ( , „ \ J \r\^cos{2nAf). (199) 
C \\r\^ cos^ ttM J 

In the opposite case, 

A > £;e|r|^ (200) 

mesoscopic fluctuations of the differential capacitance exceed the amplitude 
of oscillations of the average (C(jiff(7V)). 

To characterize the mesoscopic fluctuations of the differential capacitance, we 
calculate the correlation function K{s) — {SCdis{J^ + s)SCdis{J^)) /C^, see 
Section 4.6. For temperatures T ^ A, we find for the leading term in K{s) 

^W-y^-(2..)(|)'^ln^(f), (201) 



see Eq. (321) in Section 4.6; here P = 1(2) for a system with (without) time- 
reversal symmetry. The periodic oscillations of K{s) are due to the Coulomb 
blockade; they remain correlated over an interval s ~ Ec/A. Note, that for 
the case of one single- mode junction we consider here, the periodic part of the 
correlation function exceeds the non-periodic part [see Eq. (337) with A^^h = 2] 
by a large factor ~ 10(T/A) ln^{Ec/T). If the temperature becomes lower than 
A, the periodic part of K{s) can be estimated by Eq. (201) with T replaced 
by A under the logarithm, and for the estimate of the non-periodic part one 
may replace T by A in Eq. (337). The periodic part of the correlation function 
dominates at the lowest temperatures as well. With the increase of the number 
of channels N^h, however, the relative magnitude of the periodic part drops 
rapidly (see Section 4.6 for details). 



4-2 Finite-size open dot: Introduction to the bosonization technique and the 
relevant energy scales 



We start with a single-junction system, see Fig 17. The relevant electron states 
which facilitate the electron transfer through the junction and thus deter- 
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Fig. 17. Schematic view of a quantum dot with a single junction 



mine the charge fluctuations, can be labeled by the momenta measured in 
the narrow region at the contact, referred to as channel. This way, the prob- 
lem of charge fluctuations in the system of interest is reduced to an effective 
one-dimensional problem [39] . The corresponding Hamiltonian consists of two 
parts. The first part describes non-interacting fermions moving in one spatial 
dimension. Without backscattering in the channel, and in the limit A — > 0, 
this part takes form 



E 

(7 



r 1 



(202) 



Here ipl{x) and ipcr{x) are the creation and annihilation operators of the one- 
dimensional fermions in spin state cr; the half-axes (0; oo] and [— oo; 0) corre- 
spond to the dot and the combined channel plus lead system, respectively. We 
are interested in the quantum dynamics of the system at energy scale ~ 
which is small compared to the Fermi energy. This allows us to linearize the 
spectrum of the fermions. Writing ipa{x) = e~'^^''^ipLu{x) +e^^^^ipR(j{x)-, where 
'ipLcr and 'ipRa are left- and right- moving fermions respectively, we obtain from 
Eq. (202) 

oo 

^0 = ivFY^ i dx (ijl^d^ci'La - i/jR^d^ipRcr) . (203) 



Here Vp is the Fermi velocity of one-dimensional fermions. The possibility of 
reflection in the channel, which results in a deviation of the channel conduc- 
tance from its perfect value e^/irh, can be taken into account by addition of 
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a backscattering term to the Hamiltonian (203), 

i^bs = Irli;^^: (^L(0)^^.(0) + 4,(0)^^.(0)) , 



(204) 



where |rp <^ 1 is the reflection coefficient .P^ 

The second part of the effective Hamiltonian represents the charging energy. 
Until now, we were expressing it in terms of the number of electrons inside 
the dot. Using the conservation of the total number of electrons in the entire 
system, it is convenient here to express this energy in terms of charge outside 
the dot, 



Hn = Er 



dx 



(205) 



(Here : . . . : denotes normal ordering of the fermion operators.) We would 
like to treat the charging energy (205) non-perturbatively, while developing a 
perturbation theory in ifbs- Refs. [38,39] suggest a scheme in which this can 
be done, using the bosonized representation for the one-dimensional fermions 
in the channel. The boson representation makes the interaction Hamiltonian 
(205) quadratic in the new variables, thus removing the main obstacle in 
building the desired perturbation theory in the backscattering Hamiltonian 
(204). 

In the boson representation, the two parts Hq and -ffbs of the Hamiltonian 
take form: 



-(V0,)2 + 2(V^,)^ 



TT 



|L'cos[20F^p(O)] cos[2v^^,(0)]. 



(206) 
(207) 



where D is the energy bandwidth for the one- dimensional fermions, which are 
related to the boson variables by the transformation [139]: 



(x) 



D 



2txvf 



exp 



ix) + ejx) + aes{x] 



2^ The operators ^(0) and ^"^(0) in Eq. (204) are regularized as / dxf{x)ip{x) and 
/ dx f {x)'il^^ {x) , respectively, where / is a positive function, symmetric around x = 0, 
with / dxf{x) = 1, and with support close to the origin. This same regularization 
is used throughout the remainder of this section. 
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D 



(x) 



2ttvf 



exp < l\'TT 



Here we introduced Majorana fermions 't)±i here to satisfy the commutation 
relations for fermions with opposite spins, {57+1,^7-1} = 0, f(j^^ = 1. Anti- 
commutation of electrons of the same spin (a = lor — 1), is ensured by 
the following commutation relations between the canonically conjugated Bose 
fields: 



= —i6yY6{x — x'); 7,7' = p, s. 



(208) 



The interaction term (205) becomes also quadratic in the boson representation: 

1 2 



He 



26,(0) 



TT 



(209) 



The operators (2e/ ^)V6p{x) and {2 / y/Ti)V 6 s{x) are the smooth parts of the 
electron charge (p) and spin (s) densities, respectively.Pl 

The average charge of the dot can be determined from the AA-dependence of 
the thermodynamic potential Vt of the full Hamiltonian. Taking into account 
that only one of its parts, the charging energy, depends on A/", we find: 



1 dVL 

2E-W' 



(210) 



(hereinafter (. . .)g indicates the average over the quantum state without disor- 
der average). In the absence of backscattering, r = 0, the Hamiltonian of the 
system Hq + Hq is quadratic in the spin and charge densities. The explicit M- 
dependence of the Hamiltonian can be easily removed by the transformation 
dp{x) 9'p{x) + (l/2)Ar7ri/2. Hence, in this case the thermodynamic potential 
fl, and hence the ground state energy, obviously have no AA-dependence, the 
average charge is linear in A/" and is not quantized, and 

^(^,(0)),= (iV),=Ar. (211) 



Note that scattering between the right- and left-moving states occurs only 
at the point of a barrier in the channel x = 0. The Hamiltonian (202)-(204) 

Equation (209) contains no contribution from the boson fields at —00, because of 
the order of limits implied in the interaction Hamiltonian (205). First, the range of 
the interaction is taken to infinity, then the length of the one-dimensional channel. 
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completely ignores the fact that the dot has a finite size, and a particle that 
entered it, eventually ought to exit. However, these two events are separated 
by the dwell time ~ /i/A, which makes the model (206)-(209) applicable at 
A — > 0. To illustrate the effect of a finite dwell time, we replace the infinite 
interval (0, oo] by a finite segment (0, L], choosing the length of the effective 
one-dimensional system in such a way that the time delay of a particle entering 
this segment at a; = and bouncing off the "wall" at a; = L equals /i/A. In 
other words, throughout the remainder of this section we replace the upper 
limit of the integration in Eq. (206) by L, 



^0 = Y j dx J2 



(212) 



where L ~ nhvp/A, and set zero boundary conditions for the displacement 
fields in the spin and charge mode: 

e,{L)^e,{L)^o. (213) 



When the refiection amplitude r is finite, the Hamiltonian of the system is 

H = Ho + Hc + i^bs- (214) 



The largest energy scale appearing in the Hamiltonian (214) is the charging 
energy E^. At energies smaller than Ef. the charge is pinned to the vahic of 
gate voltage J\f, and does not exhibit quantum fluctuations. Alternatively, in 
the language of Eqs. (206)-(209), the charge mode 9p is pinned at x = by the 
charging energy (209). The backscattering term (207) induces only negligible 
perturbations in this mode, so that the dynamics of the charge mode 6p is 
described by the quadratic Hamiltonian + -^C; as in the absence of the 
backscattering term. This observation allows us to reduce the Hamiltonian 
(214) to an effective Hamiltonian for the spin mode 9s, by averaging (207) 
over the Gaussian fluctuations of 9p around its average (211). 

This averaging yields [39] 



{cos[2V^9p{0)])g = cos[2V^{9p{0))g] exp 

= cos TTjV' 



-2n {{9pm) 



/2eC£;, 



ttD 



where C ~ 0.577 is the Euler constant. As a result the effective Hamiltonian 
for spin dynamics is given by 
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-(V0,)' + 2(V^,)2 



H^^ = -JDE^\r\ cos(7rA/') cos[2v^^,(0)]. 



(215) 
(216) 



The Hamiltonian (215), (216) represents a one-mode, g = 1/2 Luttinger liquid 
with a barrier at x = 0. We will now discuss the effect of backscattering at 
the barrier, as described by the Hamiltonian H^^. 

First, at not too low temperatures (to be defined more precisely below) the 
backscattering Hamiltonian (216) can be considered perturbatively. In the 
lowest non- vanishing order one finds for the backscattering-induced correction 
to the thermodynamic potential: 



6Q 



l/T 



Tr 






Tr 







(217) 



Calculation of the trace in the numerator of this formula reduces to the eval- 
uation of the correlation function: 



(cos[2v^e,(x = 0, r)] cos[2v^^,(a; = 0, 0)])^ 



D sin ttTt ' 



which is readily performed with the help of quadratic Hamiltonian (215). The 
result of this calculation is 



4e^ 



I^EcCos^A/'ln ' 



T 



(218) 



At low temperatures this correction diverges, which signals the breakdown of 
the simple perturbation theory in the backscattering Hamiltonian. In order 
to find the low energy cut-off of the logarithm, we invoke the following renor- 
malization group arguments.p^ The Hamiltonian Eqs. (215), (216) represents 
spin fluctuations in a wide energy band D. We can integrate out the high- 
energy part of this band, thus deriving an effective Hamiltonian acting in a 
reduced band of width D. The backscattering at the barrier, described by the 
Hamiltonian H^^, is known to be a relevant perturbation [140]. Indeed, aver- 
aging of the cosine term (216) over the "fast" degrees of freedom with energies 



The model is exactly soluble [39], however, the arguments we use seem to be more 
physically transparent and suitable for further purposes. 



83 



£^ D yields 



(cos[2v/^e,(0)]) 



D>e>D 



exp 



1 ^dk 



D 



cos[2V7r^,(0)] (219) 



(we use the same notation, Og, for the fields defined in the initial band and in 
the truncated band). Therefore, the amplitude of the effective backscattering 
potential U^s scales with the band width D as 

Wbs oc -/e^ |r| cos{TTAf), (220) 



and the effective scattering amplitude U\,s/D grows proportionally to Ec/D 
with the band width reduction. Eventually, at sufficiently small band widths, 
the amplitude Wbs exceeds D, signalling the crossover from weak backscattering 
at higher energies, to rare tunneling events at lower energies. The crossover 
energy Tq{N') can be found from the condition Z^/bs ~ Z). It yields, up to a 
numerical factor, the following energy scale: 

To [N) = {%e^/'K)E^\r\^ cos^ ttN . (221) 



The numerical coefficient in Eq. (221) could not be found from the RG scheme 

but results from the exact consideration. We will not describe the calculation 
of this coefficient. The renormalization procedure must be stopped when the 
band width D reaches Tq^M). 

At lower energies the dynamics of the spin mode is described by the renor- 
malized Hamiltonian — i^l^^} + U{Os) with the "potential energy" 

t^{^.(^)} = -A2To(AA)^^^^^cos[2v^^,(0)]+^,. / dxiVOsf .{222) 

cos 7r7V J 



The spin mode 6*^(0) is strongly pinned by the potential energy. Physically, 
it means that the fluctuations of the spin of the dot are rare events of hops 
between the local minima of the potential (222) [the numerical coefficient in 
Eq. (222) is A2 ~ 1]. 

The "potential energy" (222) can be easily minimized with respect to the 
possible spatial configurations of the field 9s{x), which satisfy the boundary 
condition (213), and have the value ^s(O) at x = 0: 

C/min[^.(0)] = -A2ro(AA)^^4^ cos[20F^,(O)] + OM^N^- (223) 
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Fig. 18. The potential C/min[^s(0)] at various gate voltages, (a) At cos(7rAA) > the 
potential has one minimum (spin singlet state), (b) At cos(7rAA) < the minimum is 
double-degenerate (spin doublet state). Tunneling between the two minima corre- 
sponds to the hybridization of the spin states of the dot and lead. This hybridization 
and consequent lifting of the degeneracy essentially is the Kondo effect. 

In the limit L — oo {i.e., A — 0) the second term of Eq. (223) does not 
contribute, and the minima of the potential energy correspond to 



2v^^,(0) = 2nn , if cos{nN') > , (224) 
2^^,(0) = {2n + 1)71 , if cos(7r7V) < 

with integer n. In each of these two sets, the energy is the same within a set 
for all n. At large but finite L, when the condition 

To (AT) > A (225) 



is satisfied, the minima of energy (223) hardly shift, but their degeneracy is 
reduced strongly. If the number of electrons on the dot is close to an even 
number, cos(7rA/') > 0, then the degeneracy is removed completely. If the 
charge is closer to an odd number of electrons, cos(7rA/') < 0, the energy 
minimum preserves double degeneracy {n = —1 and n = 0), see Fig. 18. 

The "kinetic" part K{(j)s} of the Hamiltonian is quadratic in V0s, see 
Eq. (215); it does not commute with the "potential" part, and causes tun- 
neling between the potential minima. The tunneling amplitude is energy- 
dependent [141], and small at E Tq^M). To describe the low-energy dy- 
namics of the spin mode, it is convenient to project out all the states of the 
Luttinger liquid that are not pinned to the minima of the potential (223). 
Transitions within one of the sets (224) then can be described [140] by a 
tunneling Hamiltonian 

^± = - OTrTfA/-^ {v^[0.(+O) - 0s(-O)]} , (226) 
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which operates in an energy band D -C Tq{M). Here a discontinuity of the 
variable (psix) at x = is allowed, and the point x = is excluded from the 
region of integration in Eq. (215). The Hamiltonian Hxy, which is a sum of two 
operators of finite shifts for the field 6*3(0), represents hops 6*5(0) — > 6*^(0) ± 
between pinned states. In other words, H± is the finite shift operator with 
respect to the variable 6 as can be checked with the help of the commutation 
relations (208). 

These hops correspond to a change by 1 of the z-projection of the dot's spin. 
At energy scales D ^ A the hops are predominantly inelastic. The variation 
of the transition amplitude with energy can be found by continuing the same 
RG procedure, until the running cut-off D reaches the value A. At energies 
below A, fiuctuations which involve the excitations of spin mode in the dot 
(0 < X < L) are suppressed. If the minimum of energy (223) is non-degenerate 
(spin-singlet state of the dot), then only one spin state in the dot remains avail- 
able. However, in the case cos(7rA/') < 0, there is an exact degeneracy between 
two lowest-energy states, see Fig. 18; the ground state is a spin-doublet. The 
Hamiltonian (226) hybridizes the spin of the dot with the continuum of spin 
excitations in the lead. The Kondo effect consists essentially of this hybridiza- 
tion, which lifts the degeneracy and ultimately leads to the formation of a 
spin singlet in the entire system. The energy scale at which the hybridization 
occurs, is the Kondo temperature of the problem at hand. 

Note that Eq. (226) does not contain any unknown numerical factor. Finding 
the multiplicative factor in (226) cannot be done within the RG treatment we 
presented above. To find the numerical constant, one has to use [135] the exact 
result [140,40] for the linear response of an unconstrained (L 00) Luttingcr 
liquid with g = 1/2. Matching the low-temperature asymptote of the exact 
result with the result of the calculation based on the tunneling Hamiltonian 
formalism, we were able to find the constant in Eq. (226). Also, we would 
like to notice that the initial Hamiltonian (206), (207), (209) commutes with 
the total spin of the system: scattering in the junction, represented by the 
Hamiltonian (207), is clearly spin- conserving. The spin Hamiltonian (215)- 
(216), obtained after the charge mode was integrated out, also preserves the 
total spin. However, the tunneling Hamiltonian (226) has the symmetry of the 
operator S'2:(— 0)S'2:(-|-0) + Sy{—0)Sy{+0), and represents the in-plane part of 
exchange between the spin densities of the lead and dot [S'(+0) and S'(— 0), 
respectively] at the point of contact. To restore the SU(2) symmetry of the 
Hamiltonian, we must supplement H± with another term, which represents 
the Ising component of the exchange interaction: 

= t^tt^Lt V^s(-0)Ve,(+0). (227) 
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This term is omitted usually in the DC tunneling problem for a Luttinger 
liquid [140]. In agreement with the general idea about SU{2) symmetry of 
the g = 1/2 Luttinger liquid [142], both operators (226) and (227) have the 
same scahng dimension. To find the numerical factor in Eq. (227), again a 
comparison with the exact solution is necessary. 

We have seen in this Subsection, that there is a hierarchy of the energy scales 
characterizing an open dot. The charge of the dot fluctuates only at energy 
scales exceeding E^, at lower energies these fluctuations are suppressed. How- 
ever, strong fluctuations of the spin of the dot continues down to lower energies, 
To(A/') see Eq. (221). If To{J\f) exceeds the level spacing A, then at energies 
A ~ e ~ To (A/") the spin of the dot is quantized, and may change only in 
quanta of SS = 1 due to the exchange interaction with the electrons in the 
lead; at the same time, presence of many intra-dot excitations (electron-hole 
pairs) is allowed. Finally, at energies below A the intra-dot excitations are also 
suppressed. Under these conditions, the spin of the dot is 5* = or 5* = 1/2, 
depending on the sign of cos ttA/". If this cosine is positive, then the dot is in a 
singlet state, and there are no peculiarities in the temperature dependence of 
observable quantities at T ^ A. However, if costtA/" < 0, then the dot carries 
a spin, and the Kondo effect develops at sufficiently low temperatures. In the 
next Subsection we will study the temperature domain A, find the Kondo 
temperature, and discuss the transport properties of the dot in Kondo regime. 
Starting from Subsection 4.4, we return to the intermediate temperature scale, 
in order to discuss the dominant transport mechanisms there. 



4-3 The limit of low temperature: The effective exchange Hamiltonian and 
Kondo effect 

We could continue using the same RG method applied to the bosonized rep- 
resentation of the Hamiltonian, down to scale D ~ A; it would allow us to 
estimate the order of magnitude of the cff'cctive exchange constant J responsi- 
ble for the hybridization. It is possible, however, to develop a different method 
applicable in the energy domain To{J\f) ^ A, which allows us a more ac- 
curate determination of J and of the Kondo temperature Tk- 

At energies E <^ To(A/'), the spin field 9s{0) is pinned at the point contact. 
RecaUing that 9s{L) — 0, we see that the spin of the dot indeed takes discrete 
values only, as was mentioned above. At such energy scale, it is instructive to 
return to the fermionic description of the problem. After the two parts of the 
Hamiltonian, Eqs. (226) and (227), are found, one can explicitly see that the 
initial SU{2) symmetry of the problem is preserved. Therefore, the effective 
Hamiltonian controlling the spin degrees of freedom of the system dot-|-lead 
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at low energies, corresponds to the isotropic exchange interaction, 

4, = J{D)^^Sd. 



(228) 



Here Sd = ^xU{Ro)^aia2Xa2{Ro) and S = |^|^(^o)^ai(72^a2(-Ro) are the op- 

— * 

erators of spin density in the dot and in the lead respectively, at the point Rq 
of their contact; (To-i(T2 are the Pauli matrices. In what follows, we adopt this 
point to be the origin, Rq = 0. The dimension of the exhange constant J{D) 
depends on the way we normalize the electron wave functions; dimension-wise, 
it can be written as J{D) oc [i/di/To]"^, where and v are the densities of 
states in the dot and lead, respectively. 

In fact, the form of the Hamiltonian (228) could be derived without resorting 
to the explicit calculation performed above within the bosonization technique. 
The only important notion needed for the derivation, is the low-energy decou- 
pling of the dot from the lead. Pinning of both charge and spin modes is the 
manifestation of this decoupling. The decoupling is complete, however, only in 
the absence of excitations. To describe the behavior of the system at low but 
finite temperatures, one needs to establish the least irrelevant operators that 
violate the complete decoupling. The decoupled dot and lead are described 
by independent Fermi liquids. The number of electrons in the dot can not be 
changed, as such a variation of charge would result in a large energy increase 
of order Ec- That is why the leading irrelevant operators preserving the initial 
SU (2) symmetry are the interaction in the charge and spin channels between 
the dot and the lead at the point of contact, 

Hirr = A{D)pPd + ^ex, 

where ^ex is defined in Eq. (228) and pa = xl{Ro)Xa{Ro) and p = ipl{Ro)ipa{Ro) 
are the operators of the charge density in the dot (a; > 0) and in the lead 
(x < 0) respectively, at the point Rq of their contact. The first term here 
represents the interaction in the singlet channel (inelastic scattering between 
an electron in the dot an electron in the lead) . Because the charge fluctuations 
are suppressed at energy Ec which exceeds significantly the pinning energy of 
the spin mode r'^Ec, one finds A/ J ~ <^ 1. Therefore, we can neglect the 
interaction in the singlet channel. Moreover, even if this mode is taken into 
account, it does not lead to any non-trivial effect at low energies. 

The Hamiltonians (226), (227) and (228) describe, in the bosonic and fermionic 

representation respectively, the low-energy properties of the very same sys- 
tem. Correlation functions of observable quantities, calculated with the help 
of Eqs. (226), (227) should coincide with the results for the same functions 
obtained from the calculation which uses Hamiltonian (228), at temperatures 
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A <^ T <^ T(}{Af). This fact will enable us to find the exact value of the 
coupling constant in the Hamiltonian (228). The advantage of the fermionic 
representation is that it allows for the description of the system at T ~ A 
(ultimately leading to the Kondo effect). We turn to this problem now. 

We first rewrite the Hamiltonian (228) in the basis of electron eigenfunctions 
(y9„(R) and (pk(R) in the dot and lead respectively. In this basis, the Hamilto- 
nian takes the form 

-^ex — ^ Jnin2kik2 ip^ aia'i'^ ascri) ^nia\^n'2<J2^kicz^k2aA (229) 



with the matrix elements Jnin2kik2 = J {D)(p'^_^{0)(pn2{^)Vkii^)'^k2{^) ■ The de- 
tailed structure of the one-electron functions in the continuum of states {k} 
is not important in the discussion of the Kondo effect. An examination of the 
standard calculation [114] shows that only the densities |99fc(0)p, averaged over 
the continuum of states k, matter. Therefore we will suppress the indices ki 
and /c2 in the exchange constants hereinafter, assuming that the functions (pk 
are normalized so that 

^ {\Mm - (230) 



If the dot is in the spin-doublet state, the exchange constant gets renormalized 
at low energies. Unlike the "ordinary" Kondo model with only one localized 

orbital state involved, here the renormalization occurs due to virtual transi- 
tions in both the continuum spectrum of the lead and the discrete spectrum of 
the dot. To see the signature of the Kondo effect, we calculate the second or- 
der correction to the exchange constant Jqo where labels the singly occupied 
level of the quantum dot. 

This second order correction, diagrammatically shown in Fig. 19, has the form 
6ME,D)^j\D)\M0)\'J'^d^ E ,„ ^Zw^^\io)r (231) 

where the integral is taken over the continuum energy spectrum ^ in the leads, 
u is the corresponding density of electron states. The sum over the discrete 
energy levels of electrons in the dot includes the term involving the singly 
occupied level / = 0. This term, which corresponds to a spin-flip within the 
same orbital level, C,f^ = 0, is responsible for the logarithmic singularity in 
the "ordinary" Kondo effect. The rest of the terms in the sum correspond 
to the virtual transitions onto (from) the partially occupied level from (to) 
doubly-occupied (empty) electron levels. 
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G4,ki G",k" G2M2 ^2X2 ^"^k" a4,ki 



G4,ki ^ a2,^2 



— T ^ ^ 



Fig. 19. The diagrammatic representation for the second order correction to the 
exchange interaction constant. The electron states in the intermediate hnes include 
excitations within energy strip D 

At energies E ^ A the discreteness of the spectrum is not important. If one 
neglects this discreteness all together the correction to the exchange constant 
is ultraviolet divergent oc D. 



D D 



5J=J^ J ud^ 

\E\ 



-D 



\i\ + \u 



^> A. 



(232) 



Here Vd = 1/A is the one particle density of states in the dot. 

This divergent correction is similar to the irreducible vertex function in Fermi 
liquid theory and can be incorporated into a new constant J which can be nu- 
merically different from the bare constant. It is this "new" constant J through 
which all observables at low energies are expressed. This condition allows us to 
exclude a large portion of the quasi-continuous spectrum of states in the dot 
at energies E ^ D, and to relate the "new" constant J{E, D) in the fermion 
Hamiltonian, which is valid at energies E <^ D, to the constant in the boson 
Hamiltonian (226), (227), which does not have free parameters. To establish 
such a relation, one needs to consider a quantity that does not exhibit strong 
mesoscopic fluctuations, and therefore can be expressed in terms, insensitive 
to the particular values of the wave functions V5„(0). 

The observable quantity to look at, is the spin magnetization of the dot caused 
by a magnetic field applied only to the lead. Application of such a field re- 
sults in a spin polarization, due to the Pauli susceptibility, of electrons in the 
lead. This polarization, via the contact interaction, induces a finite spin po- 
larization of the dot. Such a calculation can be performed both in the bosonic 
and fermionic variables. In the boson variables, the Hamiltonian causing the 
perturbation of electrons in the lead has the form dx hzV9s{x). This form. 



90 



together with the Hamiltonian of the free spin mode, Eq. (215), allows us to 
find the induced non-zero value of (V^s(O)) ~ hz/vp- This polarization, in 
turn, creates a perturbation [vFhz/TQ{N')]V6s{+Q) for the spin mode in the 
dot, as seen from the Hamiltonian (227). Finally, the spin polarization of the 
dot S = JodxVes ~ h,/To(X) can be evaluated by considering the just 
mentioned perturbation together with the unperturbed Hamiltonian (215). 
The final result is L-indcpcndcnt and thus preserves its meaning in the limit 
A — > 0. In a Fermion representation, one finds S ~ hzJi^a^ for the spin of 
the dot by a straightforward first-order perturbation theory in the exchange 
Hamiltonian (229). Comparing the two results for S, we find 

J{E, b)^J= [vavn{N)]-^ , L> > > A. (233) 

With the help of Eq. (233) we can relate the constant J{D) in the Hamiltonian 
(228) with To(AA). 

After we incorporated all ultraviolet contributions into the definition of the 
exchange constant, we are ready to study the Kondo effect per se. Taking the 
difference of Eqs. (231) and (232) to find the second order correction to the 
"new" coupling constant (233), we find that the correction now only involves a 
contribution from the difference of the true (discrete) spectrum of the dot and 
the approximation where the density of states is constant (and continuous): 




(234) 



where t/^ = 1/A is the one particle density of states in the dot. The right hand 
side of Eq. (234) is now perfectly convergent in the high energy. On the other 
hand, the contribution of the term with I — gives a logarithmic divergence 
at i? ^ 0. This is precisely the term which describes the Kondo effect, since it 
corresponds to a spin flip of an electron on an upper level without excitations 
of electron-hole pairs in the dot. 

We obtain from Eq. (234) 



ME, b) = J|(/^o(0)r + J""^ (ko(0)|Mn ^ + a) , (235) 
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A = jim lim / ud^ 



D->oo 



E 



D 



Id" I 



(236) 



The limit E = exists in Eq. (236), and the constant A is of the order of unity 
and fluctuates from sample to sample. We will not consider these fluctuations, 
because fluctuations of the wave function v'o(O) controlling the logarithmic 
Kondo correction to Joo{E,D) lead to much stronger effects. Equation (235) 
demonstrates that in the strong tunneling regime the bandwidth for the effec- 
tive Kondo problem at hand is A rather than Ec- Once Eq. (235) is established, 
one can obtain, following the hues of [122], the result for Kondo temperature 



exp 



aA 



(237) 



with To{Af) defined in (221) and a = \ipo{0)\^ / {\^o{0)\'^) . The parameter a is 
a random quantity obeying Porter- Thomas statistics [59], see Eq. (156). 

The above consideration of the Kondo effect in a partially open dot shows 
us that the low-temperature behavior of the blockaded dot has no qualitative 
changes with respect to the case of a dot weakly coupled to the leads, as long 
as 

Ec|r|^ > A. 



The only difference developing at larger conductances, is that the parameters 
defining the low-temperature behavior differ from the bare conductances of the 
junction. Comparison of the results for Kondo temperature (172) and (237) 
illustrates this point. In both cases, Tk is controlled by the value of Ju, 

Tk ^ A(Jz/)i/2g-i/J._ 



In the case of strong Coulomb blockade, this parameter is simply related to 
the conductance of the junctions and to the ratio Ed A by second-order per- 
turbation theory in tunneling amplitudes, see Eq. (172). In the case of strong 
charge ffuctuations, we needed a more sophisticated procedure of rcnormal- 
ization in order to relate Jv to the bare parameters of the quantum dot, see 
Eq. (233). The two expressions for the Kondo temperature Tk match each 
other at refiection coefficient |rp ~ 1/2. 

The Kondo effect in a single-junction system results in a specific behavior 
of the spin polarization. If the dot is in a singlet state, the gap for its spin 
polarization is ~ A. In the doublet state, the contribution of the dot to the 
susceptibility at low temperature and fields, T^jisH <^ A, is identical to 
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that of a Kondo impurity [118] with Tk given by Eq. (237); here /is is the 
Bohr magneton for the electrons of the dot. The manifestation of the most 
interesting effect, the enhanced low-temperature conductance, requires a two- 
junction dot geometry, and we turn to the corresponding discussion now. 

Our first goal is to write a fermionic Hamiltonian describing the system at en- 
ergies so small that, due to the pinning of all charge and spin modes, the Fermi- 
liquid behavior is already restored. Similar to our discussion after Eq. (228), we 
obtain such a Hamiltonian using solely the symmetry arguments, and relegate 
the explicit model calculation to Appendix E. 

Once again, the interaction of the charge degrees of freedom of the dot and the 
lead are suppressed, so we can restrict ourselves to the spin-spin interaction 
only. The difference from the one-contact geometry, is that now processes 
which involve two contacts, as well as the ones involving a single contact, 
are possible. The most general form of the Hamiltonian describing all such 
processes is 

^ex = 1 1] Jija<,,<,2^<,3^4t/'i^^i(^i)xj,3(^j)x<,4(^i)t/'^2j(^i), (238) 



where indices i,j = 1,2 label the contacts, Jij — Jji are coupling constants, 

— * 

and Ri describes the location of the point contact. 

The exchange interaction constants Jij should be found from the calculation 
of a set of low-energy correlation functions. We can do the explicit evaluation 
of the constants for the strongly- assy metric set-up, |ri| — > 1, |r2| <^ 1. In this 
case, the largest exchange constant J22 is determined by the single-junction 
result, J22 = J ^ where J is defined by Eq. (233) with r2 = r. The constant J12, 
which characterizes the non-local exchange accompanied by the charge trans- 
fer, is proportional to and the exchange constant at the weaker junction 
is Jii oc where gi ~ 1 — |ri p is the dimensionless conductance of the weaker 
junction. As a matter of fact, in the limit of strong asymmetry the numerical 
coefficients can be found also, with the help of the exact solution [40] available 
for the case A = 0. To find J12, we can calculate the linear conductance G{T) 
through the dot using the Hamiltonian (238), and then compare it with the 
exact result [40]. In the lowest order of the perturbation theory, transitions 
across the dot are induced by only the terms cx J12 in the Hamiltonian (238). 
The conductance can be calculated by applying the standard Fermi golden 
rule to the problem. Because we neglect the level spacing at the moment, and 
because of the four-fermion structure of the Hamiltonian (238), the electron 
transition rate and the conductance G at low temperatures is proportional to 
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G{T) 



The comparison with the exact result [40], see also Eq. (359), yields: 

32 Tih 



'J 1 



12 



37r^e 



2^C 



2 2 



(239) 



The smallest exchange constant in the problem, Jn, does not affect the con- 
ductance in the strongly-assymetric set-up which we consider here. For calcu- 
lating the conductance in the absence of strong asymmetry, one should know 
the value of Jn. For that purpose, some other response function should be 
calculated with the help of Hamiltonian (238) and then compared with an ex- 
act result; we will not address this problem in this Review. Instead, we utilize 
Eq. (239) to determine the conductance through the dot in the limit of low 
temperatures. If the gate voltage is close to an odd integer, costtA/" < 0, the 
lowest discrete energy level in the dot is spin-degenerate. At T ^ A, only this 
level remains important. This way, the initial problem of the dot, which has 
a dense spectrum of discrete levels, and is strongly coupled to the leads, is 
reduced to the problem of a single-level Kondo impurity considered in Sec- 
tion 3.3. We can adapt the result (168) to find the conductance. First, we 
express the temperature-independent factor in terms of the exchange con- 
stants!^ rather than in terms of the conductances gi and g2'- 

(91 + 92^ ^ (^11 + ^22)2' 

Next, we neglect the small term Ju in the denominator here, and find: 




Now we are ready to find the Kondo conductance through the dot in the regime 
of strong charge fluctuations which we consider in this Section. For that, we 
substitute the exchange constants (233) and (239) into the above formula. 



The following formula is valid only for the single-channel Kondo model, i.e., if 
the condition = J11J22 holds. Within the y^~accuracy of our calculation this 
relation trivially holds. Checking it to higher order in gi requires calculation of Ju. 
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Here, Tk(A/") is given by Eq. (237), and f{x) is the universal scaling function 
plotted in Fig. 14. Note that the Kondo conductance (240) in a strongly asym- 
metric set-up is proportional to the product of the small conductance Gi and 
small refiection coefficient |r2p and therefore is significantly smaller than the 
conductance quantum e^/Trh even at T — 0. 

Similar to the case of the closed quantum dot, the dominant mechanisms 
switch from inelastic co-tunneling at high temperatures to the elastic tempera- 
ture-independent mechanism at lower temperatures, and finally to the regime 
of the Kondo effect at the lowest temperatures. We postpone a more detailed 
discussion of those crossovers until after we have constructed a theory of elastic 
co-tunneling in open quantum dots. We turn to the outline of this theory in 
the following subsections. 



4-4 Simplified model for mesoscopic fluctuations 

As was already mentioned in the beginning of this section, the main difficulty 
for building a theory of interacting electrons in an open dot stems from the 
need to treat both the Coulomb interaction and the dot-lead conductance non- 
perturbatively. In the limit A ^ 0, the bosonization procedure allows one to 
treat exactly the interaction in a dot connected to a lead by a perfect quantum 
channel. Then, the effect of a finite reflection amplitude r in the channel can be 
treated by a systematic perturbation theory [39,40]. An important conceptual 
conclusion drawn from this theory is that in the low-energy sector (defined by 
e ^ Ec\r\'^), the excitations inside and outside the dot are independent from 
each other. This understanding allowed us to consider the low-temperature 
effects (T ^ A) in a partially open dot with £'c|rp ^ A, because in this case 
one is able to map the properties of a partially-open dot onto the "conven- 
tional" case of strong Coulomb blockade (see the previous subsection). Such 
a mapping allows one to consider the Kondo effect, but is insufficient for eval- 
uating the co-tunneling contribution to the conductance, for which a broad 
range of energies up to Ec is involved. (This was shown, e.g., in Sec. 3.2, where 
we considered the valley conductance through a strongly blockaded quantum 
dot.) The aim of this and the following Sections is to build a theory capable of 
describing simultaneously many-body effects and mesoscopic fluctuations at 
the intermediate range of energies, i?c~ ^~ m£ix{Ec\r\'^ , A). This section deals 
with a qualitative picture of the role of the mesoscopic fluctuations; a rigorous 
calculation to verify the qualitative picture below is given in the following 
sections. 

Let us first consider a completely opened channel (r = 0). In the limit A — 0, 
the electron charge of the dot varies with the gate voltage as {Q) = ej\f, to 
ensure the minimum of the electrostatic energy [39] . The interaction term in 
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Eq. (92) depends only on the number of electrons crossing the dot-channel 
boundary. Since this information is represented as well in the asymptotic be- 
havior of the wave-functions far from the entrance to the dot, such as the 
scattering phase shift, the properties of the ground state can be characterized 
by an analysis of those asymptotes. For the sake of simplicity, we consider 
the case of a dot connected to a lead by one channel, and consider spinless 
fermions instead of real electrons with spin. The ground state properties at 
low energy can then be understood from the following argument. 

The entrance of an additional electron with energy e (all the energies will 
be measured from the Fermi level) into the dot requires the charging energy 
Ec- Therefore, the electron may spend a time of the order of TijEc in the 
dot, after which the extra charge accumulated on the dot has to relax. One 
can distinguish two processes that lead to charge relaxation: (i) elastic process 
where the same electron leaves the dot, and (ii) inelastic process where another 
electron is emitted from the dot. At low energies the probability of the inelastic 
process is small as [ejE^'^ ^ by virtue of the smallness of the phase volume. (The 
last statement assumes Fermi liquid behavior at low energies and, as we will 
see later, is vahd only for the spinless one-channel case.) Therefore, for e <^ £'c, 
we may conclude that an incoming electron spends a time of the order l/£'c 
inside the dot, after which it bounces back into the channel without creating 
excitations. The same consideration is applicable to an electron attempting 
to leave the dot. The fact that within the energy range \e\ % E^ from the 
Fermi level the electron processes are predominantly elastic implies that the 
low energy properties of the system can be mapped onto a dot effectively 
decoupled from the channel, that there is a well-defined scattering amplitude 
from the entrance of the dot, and that the phase of this scattering amplitude 
is given by the Friedel sum rule 

5 = 7r{Q)/e = nX. (241) 

Equation (241) can be applied to electrons incident from inside the dot, as 
well as to electrons incident from the channel. The description outlined here 
resembles closely Nozieres' Fermi liquid description of the unitary limit in the 
one-channel Kondo problem [118]. 

We now demonstrate that the above considerations are sufficient to reproduce 
the result 

Eg(A/') ~ |r|EcCos27rA/', (242) 

obtained by Matveev [39] for the ground state energy Eg{M) of spinless 
fermions in the limit of zero level spacing A in the dot and with a finite 
reflection amplitude r in the contact, and then apply the same scheme to flnd 
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the corrections to the ground state energy arising from a finite A, but without 
r. Those corrections will result in the mesoscopic fluctuations of the ground 
state energy. 

We start with the limit A = 0. First, we put r = and calculate the density of 
electrons in the channel p{x). As wc discussed, the Coulomb interaction leads 
to a perfect reflection of electron at low energies; the wavefunctions have the 
form ipkix) = cos {k\x\ — 6), with phase shift S given by Eq. (241). This form 
of the wavefunctions leads to the Priedel oscillation of the electron density 
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VF\k-kF \ ^ Ec 

where vp and kp are the Fermi velocity and Fermi wavevector respectively. 
(Here we omitted the irrelevant constant part of the electron density; the 

restriction of wave vectors k to the range E^/vp around kp is because only for 
those wavevectors elastic reflection from the contacts takes place.) We thus 
obtain 

Ec 

— cos{2kp\x\ — 25), \x\ < vp/Ec, 

Pi^) = {7r.(oi..M-OM (243) 

\x\ > Vp/Ec. 
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Next, we take into account the effect of a scattering potential V{x) in the 
contact that generates a finite refiection amplitude r 7^ 0. In the first order 
of perturbation theory, the presence of the potential V gives a shift to the 
ground state energy given by 

Eg{M) = J dxp{x)V{x). (244) 



Substituting Eq. (243) into Eq. (244), assuming that the effective range of 
the potential around a; = is smaller than vp/Ec, and using the standard 
expression |r| = \V{2kp)\/vp, we obtain Eq. (242). Here V{k) is the Fourier 
transform of the potential V{x) . 

Having verified the relation between scattering phase shifts and the ground 
state energy for the case A = 0, we proceed with evaluation of the ground state 
energy of a finite dot connected to a reservoir by a perfect channel. According 
to the discussion preceding Eq. (241), the channel is effectively decoupled from 
the dot due to the charging effect even though r = 0. We now find the ground 
state energy of the system by relating the ground state energy of the closed 
dot to the scattering phase S of Eq. (241). For a chaotic dot, this problem 
is equivalent to finding a variation of the eigenenergies by introduction of an 
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arbitrary scattcrcr [143] with the same phase shift 5. The relevant contribution 
to the ground state energy is given by 

Eg^ E [m+^^]: (245) 

-Ec < ii<0 

where are the eigenenergies measured from the Fermi level /x. (Again, we 
have restricted the summation to a window of size ~ around the Fermi level, 
since the phase shift 5 applies only to particles in this energy range.) As soon 
as the scattering phase changes by tt, one more level enters under the Fermi 
level. Evolution of the energy levels with changing 6 is shown schematically in 
Fig. 20. The position of the level ^i{5) satisfies the gluing condition 

US + 7r)^Ci+iiS). (246) 

From Eqs. (245) and Eq. (246) we see that the ground state energy depends 
almost periodically on 5: 

E,{6)=E,{6 + n) + 0{A). (247) 

As we will sec below, the amplitude of the oscillation of the ground state 
energy with the variation of 6 is much larger than the mean level spacing, so 
that we can neglect the last term in Eq. (247). 

In order to estimate the magnitude of the oscillations of the ground state 
energy, we recall that the correlation function of the level velocities is given 
by [47] 

{dsei) = -, (dss^dssj) = S^jl(-) , (248) 

TT p V TT / 

where [3 = 1,2 for the orthogonal and unitary ensembles respectively, and (. . .) 
stands for the ensemble averaging. Formula (248) can be easily understood 
from the first order perturbation theory: At 5 <^ 1, we have £i(5) ~ £i{0) + 
{S/7ru)\ipf{0)\. Using the fact that ipi{0) is a Gaussian random variable, see 
Eq. (21), one obtains Eq. (248). 

Estimating the mesoscopic fluctuations of the ground state energy (245) with 
the help of Eq. (248), we obtain for 5 < 1 

{[E,{S) - E.m') - A^(^)' E 1%-^- (249) 
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Phase shift, 5 



Fig. 20. Evolution of the energy levels of the quantum dot with the scattering phase 
6. 

As we have already explained, the energy (245) is a periodic function of 6 
with period vr. On the other hand, for 5^1, Eq. (249) is valid. Therefore the 
characteristic amplitude of the oscillations is of the order of \/EcA, and it is 
plausible to assume that the correlation function of energies at two different 
parameters A/i, M2 takes the form 



where we use Eq. (241). It is important to notice that the variation of the 
energy of the ground state is much larger than the mean level spacing A. 

A different, though equivalent way to arrive at Eq. (250) is as follows. We have 
seen that because of the presence of the charging energy Ec, charge excitations 
in the dot relax elastically after a time ~ h/E^. If the dot size is infinitely big 
(corresponding A = 0), no backscattering from the dot's walls occurs within 
this short time interval, and the only source of elastic scattering is an eventual 
reflection amplitude r in the contact. However, if the dot has a finite size, there 
is a finite probability that the incoming particle will scatter from its walls, or 
from impurities inside the dot, and exit elastically before the time h/Ec- The 
corresponding amplitude r for reflection from the dot in a time shorter than 
h/Ec can be estimated as 
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h/Ec 




(251) 
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where S{t) is the scattering matrix of the dot in time representation, cf. 
Eq. (119). Hence, in order to inchide elastic scattering from the dot, we have 
to add the amplitude (251) to the reflection amplitude Tc from the contact 
in Eq. (242). The mesoscopic fluctuations of Eg can then be found from Eq. 
(122) for the mesoscopic fluctuations of S{t). One verifles that Eq. (250) is 
quickly recovered. 

From the present derivation, it is clear that the charging energy Ec enters into 
the flnal estimate (250) as an energy scale below which the dot is effectively 
closed. All higher lying energy levels are not sensitive to the charging energy, 
and hence not to the gate voltage N". A similar picture holds also for "real" 
electrons with spin. As we discussed in Sec. 4.2, see Eq. (221), in that case 
the dot effectively closes at the much lower energy scale Ec\r\'^. As a result, 
the amplitude of fluctuations for spin- 1/2 fermions acquires additional factor 
r^, compared to Eq. (250). Such an estimate is vahd unless the reflection 
coefficient |rp becomes smaller than A/E^. We will see later that in the case 
of a perfectly transmitting channel (r = 0), the energy scale £^ckP should be 
replaced by ~ A In (Ec/A). 

A magnetic field applied to the system affects the states in the dot. Therefore, 
the ground state energy and, as we will see later, conductance through the 
dot varies with the field. The corresponding correlation magnetic field fiux 
$c is determined by the number of levels that fit into the energy strip of one- 
particle levels "blocked" inside the dot. An estimate of $c can be obtained from 
Eq. (161) with E replaced by Ec or Ec\r\'^ for spinless or spin-1/2 fermions, 
respectively. 



4-5 Rigorous theory of mesoscopic fluctuations: formalism for multichannel 
junctions 

4.5.1 Derivation of the effective action 

We now proceed with an exposition of the general formalism used to calculate 
the capacitance, tunneling density of states, and conductance of an interacting 
quantum dot, without neglecting the effect of mesoscopic fluctuations of the 
wavefunction in the dot. (That is, without neglecting the possibility that an 
electron that has entered the dot will return coherently to the point contact.) 
The flrst step is to reformulate the problem in terms of an effective action that 
accounts for these coherent returns. 

Our formal starting point is the same as in Sec. 2.4: The Hamiltonian of the 
dot connected to the leads is written as the sum of a Hamiltonian of the closed 
dot, a Hamiltonian of the leads, and a term Hld coupling the lead and the 



100 



dot, 



(252) 



Here, the Hamiltonian of the dot is, in turn, written as a sum of the Hamilto- 
nian Hd of the closed dot in the absence of electron-electron interactions and 
the charging energy He [see also Eq. (92); here spin is included in the indices 
a and 7] 

Hd^YI ^"tV^I^V^t^ ^E,{n- Nf . (253) 

The Hamiltonian Hj^ of the leads is given by Eq. (96); the term Hld that 
couples the dot with the leads is given in Eq. (97). 

The operator h in Eq. (253) is for the total charge eh on the dot. Conven- 
tionally, it is defined in terms of the creation and annihilation operators ij)]^ 
and ipa of fermions in the dot, as in Eq. (34). For the purpose of the formal 
manipulations of this section, it is more convenient, however, to change the 
definition (34) of n, and reformulate it in terms of the operators '0](/c) and 
i'j^k) for the channel. This is possible, because the total number of particles 
in the entire system (leads and dot) is an integer number which can be added 
to the parameter H at no cost, so that we can write 

^^-Ey ^i^jmAk)- (254) 

Strictly speaking, the use of Eq. (254) instead of (34) is only permitted for a 
canonical description, in which the total number of particles in the system is 
held fixed. Below, however, wc use Eq. (254) in a grand canonical description. 
The difference between the two descriptions is only important if non-periodic 
fluctuations as a function of A/" are considered. The modifications of the theory 
in this case are discussed in Appendix F. 

We now discuss the derivation of the effective action separately for the three 
cases of interest. 

Free energy and differential capacitance. To calculate the ground state energy, 
we start with the thermodynamic potential, 

= -Tln(Tre-^/^) . (255) 

We evaluate the trace in two steps, Tr... = Tr/^Tr^) . . ., where L and D 
indicate the fermionic operators belonging to the leads and dot respectively. 
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With the redefinition (254) of the charge n, the interaction Hamiltonian He 
is now attributed to the leads. As a result, the Hamiltonian of the system 
becomes quadratic in the fermionic operators of the dot, so that this part of 
the system can be integrated out: 

Tro e-^/^ = Tr,, e-(^L+Hc+Ho+HLo)/T ^256) 
Here Hld{t) is the interaction representation of the dot- lead coupling operator 

Tj. denotes time ordering for the imaginary time r, flo is the thermodynamic 
potential of the closed dot, 

Qd = -TlnTre-^^/^, (257) 
and the average {■ ■ ■)d over the Hamiltonian Hd of the dot is defined as 

Note that, as the interaction has been shifted to the leads, the thermodynamic 
potential oi the closed dot is formally independent of the gate voltage M . 

It remains to compute the average {Hld{ti)Hld{t~2)) d in Eq. (256). Hereto 
we use the explicit form (97) of Hld and the Green function Qapij) of the 
closed dot, 

g^p{T) ^ -{T^MSpmD = T f: e— f . \ ] , (258) 

n=-oo yi^n ~ J^dJ a/3 



where a;„ = 7rT{2n + 1) is a fermionic Matsubara frequency and iP{t) 
■0'^(— t). Hence we find 



{HLD{n)HLD{r2))D = -2j2^ni'^u^) W^g{T,-T2)W . . ^^,,(t2,0), 

31,32 

(259) 

where the operators '?/'j(T, x) are the Fourier transform of '4^j{T, k), 

dk 
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It is convenient to introduce a new matrix, Lj^j2, by the following relation: 

4L(za;„),-,,, = -^-^^^5j,j, + \w^g{iujn)w] . . , (260) 



J1J2 



where u = 1/2itvf is the density of states in the leads. With this definition, 
the dot Green function Q and the couphng matrix W can be eliminated from 
Eq. (259) in favor of the matrix L, so that we obtain 



2T - ^ 

- I E 'Ai.(n,0)iju,(Tl - T2)4(T2,0). (261) 

The first term in the r.h.s. comes from the diagonal part of the Green function 
Q {i.e., the density of states in the dot); this term can be obtained by multiply- 
ing the ensemble-averaged Green function by the product M^^l^ = n'^u/MA 
characterizing an ideal dot-lead coupling. The second term, which contains the 
matrix Lj^ j,^, also has components off-diagonal in ji and j2- It accounts for a 
combined effect of electron trajectories backscattcrcd into the lead. Backscat- 
tering may be caused by electron reflection off a barrier in the contact, impu- 
rities in the contact and the dot, and the boundaries of the dot. This term is 
responsible for the mesoscopic fluctuations of the observable quantities. 

Except for the average density of states, the first term in Eq. (261) does not 
contain information about the dot. We could have obtained the very same 
term if the coupling to the dot would have been replaced by an ideal coupling 
of each chiral fermion field '4^j{x) in the lead to a free chiral fcrmion field bj(x). 
Writing these fields explicitly, we find that Eq. (256) takes the form 

Tr^e-^/^ oc Tr,, e-^'^^'^e-^-^/^T^e'^o^ (262) 



where Tr^ denotes a trace over the chiral fermion fields bj, H^s is an effective 
Hamiltonian that contains the fermion fields il^j{x) in the lead, the fictitious 
fermion fields hj{x), and their coupling at x = 0, 



He^ = ivFY / dx '4j]{x)dx'il)j{x) + h]{x)dj)j{x) 

+ Ec(^J dx ■.i^]{x)iPji,) : +Af^ , (263) 
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and <Seff is an effective action that contains the fields ipj{x) at x — only, 

l/T 1/T 

5eff = 4^ I dn i rfr2V5,,(ri;0)L,,,,(Ti-r2)4(T2;0). (264) 

jl,j2 

The equivalence of the representations (256) and (262) - (264) can be easily 
checked by tracing out the fermions bj in Eq. (262), with the help of the 
relation [cf. the first term on the r.h.s. of Eq. (260)] 

{Trbj,{T; 0)bjM 0)) = i7nyT5j,j, ^ e-^'^"^sgnu;„. 

The operators ■0^(0) and bj{0) appearing in the effective Hamiltonian (263) and 
in the effective action (264) are the fermion fields evaluated at the origin x = 0. 
To simplify notation, from now on we omit the coordinate in the argument 
of the one-dimensional operators ipj (or, later, of the bosonic operators (pj), 
when they are taken at the origin, 

^|;J{T) = iPj{T;x = 0), 0,(t) = 0,(t;x = 0). (265) 



It can be checked that the operator in the second line of Eq. (263) connecting 
the fermion fields ijj and 6 at a; = corresponds to an ideal coupling between 
those fields. This couphng means that a ^-particle is scattered to a 6-particle 
with unit probability. This motivates the following substitution of right and 
left moving fermion fields. 



1pj{x)^1pL,j{^)0{-x) + 1pRj{-x)0{x), 

bj{x)^i '4)Rj{-x)9{-x) - '4;L,jix)9{x) 



(266) 



[Here e{x) = Iff a; > 0, if x < 0, and ^(0) = 1/2.] With this substitution we 
obtain from Eqs. (263) - (264) 



oo 

1 — oo 

+ ^c|^E / dx:^l^l^^'4^L,j + '4^i,j'4^Rr.+M^ , 
1/r i/T 

»5efr = j dTi j dT2 Lj^j^ (ti - T2) 



JIJ2 



X 



V'LJi in) + IpRj, (ti) V'L,J2 (7-2) + -^i? J2 (T2) 



(267) 



(268) 
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where the time dependence of any operator A is defined as 



A{r) 



(269) 



In this formulation, all backscattering from the dot and from the contact is 
characterized by the kernel L. To make contact with the theory of Sec. 2.4 of 
a non-interacting quantum dot, it is advantageous to express L in terms of 
the scattering matrix S of the non-interacting system. Comparing Eqs. (101) 
with Eq. (260), one immediately finds 



l/T 



S{ilxln) 



27rii/ 1 -|- S{iujn) ' 
I 27rwl + ,5t(ia;„)' 



UJn > 0, 
UJn < 0. 



(270) 



This relation allows us to deduce the statistical distribution of the kernel L[t) 
from the better known statistical distribution of the scattering matrix S, see 
Sec. 2.4. For ideal contacts the ensemble average of S"" (with integer n) is zero, 
cf. Eq. (115), which implies that the ensemble average of L is zero as well. 

Equations (267) and (268) provide an effective description of the combined sys- 
tem of the lead and the dots in terms of fermions moving in one-dimensional 
wires, where the interaction only exists on one end of the wire. The effec- 
tive action <Sefr provides a coupling between the wires and between left and 
right moving fermions. The one dimensional description for the charging ef- 
fect, which was already introduced in Sec. 4.2, was first proposed by Flensberg 
[38] and Matveev [39] . Their theory can be interpreted as an approximation of 
the effective action theory presented here, corresponding to the replacement 
of S (and L) by its ensemble average. 



where Vj is the backscattering amplitude in the contact for the j-th channel, 
see Sec. 2.4. At weak backscatterring, <^ 1, Eq. (270) gives, cf. Eq. (121), 

L,^{t)^vf5{t)5,^t^, (272) 



so that the action loses its time dependence and acquires the form of a 
backscattering Hamiltonian. The approximation (271) loses all the informa- 
tion about the trajectories where the electron returns to the entrance of the 
dot, after having been reflected by the dot's walls or by impurities inside the 
dot. Such returns are the cause of mesoscopic fluctuations of the capacitance 
and other observable quantities. Mesoscopic fluctuations are present in the 
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full action (268), which is non-local in time. This action was first derived in 
Ref. [70]. 

From Eqs. (255) and (262), we find the expression for the thermodynamic 
potential of the system 

= - T In Tr {e~^^«/^T,e"'^^ff } (273) 



where fi^i is defined by Eq. (257), and the effective Hamiltonian HeS and the 
action iSefr are given by Eq. (267). 

Equation (273) allows the calculation of all thermodynamic properties of the 
system within the framework of the effective action. Such quantities involve, 
among others, the specific heat and the differential capacitance, 

Ca„(Af) = c(l-^0). (274) 



We will perform the explicit calculation of this quantity in the next subsub- 
section.p"] 

It is not difficult to extend the effective action approach to the calculation 
of kinetic quantities, such as the two-terminal conductance and the tunneling 
densities of states. We will sketch these extensions below. The corresponding 
expressions at = were discussed in Sec. 2.4. 

Two-terminal conductance, the general formulation. We will use the usual 
Kubo formula to express the two-terminal conductance G in terms of the 
Matsubara correlation function 



^ UUnn 00 + iO) - UUnn ^ 00 - iO) 

G = lim , 

uj^o 2iuj 

l/T 

n(zfi„) = j t^rn(r)e*^"^ n(r) = (r./(r)/(0))„ (275) 



where f2„ = 27mT is the bosonic Matsubara frequency, and 

/(r) = exp(^r)/exp(-i7r) (276) 



In Eq. (273) the thermodynamic potential of the dot Qd does not depend on gate 
voltage Af and, thus, does not contribute to the differential capacitance. More careful 
analysis, see Appendix F, reveals the dependence, which gives a contribution small 
as A/Ec to Cdiff. We relegate the corresponding discussion to the end of Sec. 4.6 
and Appendix F. 
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a) 



® 



b) 



® 



Fig. 21. The integration contour used in the evaluation of the conductance, see 
Eq. (275), for (a) 0„ > 0, and (b) for fi^ < 0. Branch cuts of the analytic continu- 
ation of n(r) are shown by thick lines. 

is the operator of the electric current in the Matsubara representation. Because 
the total current through any cross-section of the system is conserved we can 
define the current operator in the leads. For future technical convenience, we 
define the current / through the dot as a weighed difference of the currents Ii 
and I2 entering the dot through the two leads 1 and 2, 



where A^i (A^2) is the number of channels in the corresponding lead and 



X 



, a = 1,2 (278) 



x-*0 



is 



the operator of the current in lead a {a — 1,2). 



The fact that the current operator / is defined inside the leads allows us to 
copy the derivation of the effective action from the thermodynamic potential: 
we shift the interaction from the quantum dot to the lead and then take a 
partial trace over all dot states. As before, the result is represented by means 
of the same effective action for the one-dimensional fermions ipL, i/jr. We thus 
find 



n(r) 



Tr 


e-^«ff/^T^7(T)7(0)e 


-5eff 


Tr 


g-^eff/T7^^g-5eff 





(279) 



where the effective Hamiltonian and action are defined in Eq. (267), and the 
time dependence of the operator / is given by Eq. (269). 

Equations (275), (279) and (267) constitute the complete formulation of the 
two-terminal conductance problem within the effective action theory. The 
practical execution of the analytic continuation in Eq. (275) is more easily 
achieved in the time domain. As we will see below, the function 11 (r) can be 
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analytically continued from the real axis to the complex plane, so that the re- 
sult is analytic in a strip < Rer < 1/T, and has branch cuts along the lines 
Rer = 0, 1/T. It allows one to deform the contour of integration as shown in 
Fig. 21, thus obtaining 

oo 

u{inn) = i J dte-^"' [0{nn)e{t) - e{-nn)e{-t)] [u{it + o) - u{it - o)] . 

— oo 



Now, the analytic continuation (275) can be performed, with the result 

oo 

G^^ j dtt [U{it + 0) - U{it - 0)] . (280) 



Next, we use the analyticity of II{t) in the strip < Rer < l/T, and shift the 
integration variable t ^ t — i/2T in the first term in brackets in Eq. (280), and 
t t + i/2T in the second term. Bearing in mind that n(r) = II{t + 1/T), 
we find 

oo 

/ dm{tt + l/2T). (281) 



This formula is the most convenient for practical calculations. 

Tunneling density of states - In the case of a very asymmetric setup, [Gi <^ 
e^/irh and arbitrary G2], the general formula for two-terminal conductance 
(281) can be further simplified. One can expand the scattering matrix in terms 
of small transmission amplitude corresponding to Gi, like it was done for the 
case of non-interacting electrons, see Section 2.4, Eqs. (106)-(108). This way 
the problem of calculation of the two-terminal conductance is reduced to the 
one of evaluation of the tunneling density of states for a dot strongly coupled 
to one lead, see Figure 22. 

For the non-interacting system, the tunneling density of states is defined in 
Eq. (107). Because the transmission coefficient |tip <^ 1, we can consider 
the tunneling current / as the function of applied voltage V in second order 
perturbation theory in the tunneling Hamiltonian [first term in Eq. (106)]. 
This gives us the standard result [87] 

I{eV) = i [J {iQn ^eV + iO)-J (iQ„ ^ eV - iO)] , (282) 



where Q„ = 27rTn is the bosonic Matsubara frequency, and the Matsubara 
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Fig. 22. Schematic view of the asymmetric two-terminal setup. The left point contact 

has one channel almost open and the conductance of the right point contact Gi is 
much smaller than e^/27r?i. One of the electron trajectories contributing to elastic 
co-tunneling is also shown. 

current J is defined as 

l/T 

J{ifln) = MAp- I cire-""-ntu„(r)^^^, (283) 
27re J sihttTt 



where the last factor corresponds to the one-particle Green function of the 
electrons in the right lead (see Fig. 22; the right lead is not affected by inter- 
actions or mesoscopic fluctuations), and Htun corresponds to the one-electron 
Green function of an electron in the dot at the point of the tunnel junction, 

ntun(T) = (r,|i(T)^l(0)),. (284) 

Here we wrote ip^ and ipi for the fcrmion creation and annihilation operators 
in the dot at the location of the tunnel junction, in keeping with the RMT 
formulation of Sec. 2.4, see Eq. (106). 

Analytic continuation in Eq. (282) is performed similarly to the derivation of 
Eq. (281), and one obtains 

/ = frsinh ^) MAG, / (285) 
V 2T) J coshnTt ^ ' 

— oo 

Here ntun(T) must be analytic in the strip < Rer < l/T. The linear con- 
ductance G is therefore given by 

G^GiMA / dt \ ^ 286 

J 2cosh7rri ^ ^ 

— oo 
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One can check that at Ec = 0, Eqs. (286) and (284) are equivalent to Eqs. (107) 
and (108). The Green function (284) is dramatically affected by the interac- 
tions, and we would like to construct an effective action theory to describe 
these effects. 

Once again, we wish to get rid of the fermionic degrees of freedom of the dot. 
Similar to Eq. (254), it is convenient to rewrite the charge operator en in terms 
of the variables of the channel. However, here we have to keep in mind that 
the tunneling events described by operators ■01 and ipi change the charge in 
the dot by an amount +e and — e. This is not taken into account in the sim- 
ple redefinition (254) of n that we used for the effective action theory for the 
ground state energy and the two-terminal conductance. Instead, to account 
for the fact that the total number of particles in the system changes upon a 
tunneling event, one has to introduce three additional operators [40]: A Hermi- 
tian operator m and two unitary operators F and F^, with the commutation 
relations 



,^^1 = (287) 



These operators serve to keep track of the total number of particles in the 
system. They commute with all the fermionic degrees of freedom. We now 
include m into the definition of the charge operator. 



n^-Ej^^]ik)i^Ak) + rn, (288) 

and rewrite Eq. (284) as 

ntun(T) = {TrF{r)$, (r) F(0)V'i (0)). (289) 

It is easy to see from Eq. (287) that the operators F\F in Eq. (289) change 

the charge, as defined by Eq. (288), by +e and — e respectively, in accordance 
with the initial definition (34) of charge in terms of the fermionic operators of 
the dot. 

After this manipulation, the Hamiltonian of the system becomes quadratic 
in the fermionic operators of the dot, so that this part of the system can be 
integrated out. Manipulations analogous to the derivation of Eq. (273) give 



ntun(r)=ni„(r)+nei(r), (290) 
nin(r) ^ - g"^7| (T.e-^-«F(r)F(0))„ (291) 



110 



1/T 1/T 



1,1 ,V 



9 



X 



(292) 



where the effective action is given by Eq. (267), the Green function of the dot, 
Gapi^), is defined by Eq. (258), and the coupling matrix W^j corresponds to 
the stronger (left) junction connecting the dot and a lead. Note also that the 
scattering matrix S entering in the expressions for the effective action yS'efr, 
see Eqs. (268) and (270), includes this stronger junction only. The averages in 
Eqs. (290) - (292) are performed as 

(. . = IVe-^-/^ . . . , 
with the effective Hamiltonian 



oo 

3 — oo 

+ E,(^ J dx:7pl^7PLj + ^ij^^Rj:+J^-n)j . (293) 



The interaction representation of the operators is defined in Eq. (269), with 
Hes from Eq. (293). The difference between Eqs. (293) and (267) is in the 
different definitions of the charge operators in Eqs. (254) and (288). 

There is a good physical reason to distinguish the two contributions to the 
tunneling density of states in Eq. (290). The first contribution in (291) is 
inelastic: this term does not allow the introduced electron to leave the dot; 
the charge of the dot at the moment of tunneling suddenly changes by +e 
and all the other electrons have to redistribute to accommodate this charge. 
We will see, that the logarithmical divergence of the imaginary time action 
corresponding to such evolution (orthogonality catastrophe) completely sup- 
presses this contribution at T — > 0. Conversely, the second contribution, U^i, 
cf. Eq. (292), contains the kernel R{t) which promotes the tunneled electron 
through the dot into the channel (left lead) . Because the very same tunneling 
electron is introduced to and then removed from the dot, there is no need in 
the redistribution of other electrons, so no orthogonality catastrophe occurs. 
As a result, the elastic contribution survives at T — > 0, analogously to the 
elastic co-tunneling contribution to the conductance of a strongly blockaded 
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dot, see Sec. 3.2. 

It can be shown, that if the charging energy vanishes, Ec = 0, all physical 
results of the effective action theory contained in Eqs. (267) - (270) and (293) 
are equivalent to the non-interacting theory of Eqs. (103), (108) and (111). The 
advantage of the effective action representation becomes clear when one has to 
treat the effects of charging. The effective Hamiltonian (267) or (293) can then 
be diagonalized exactly and the effective action S^s treated perturbatively. The 
technique for such treatment is described in the following subsection. 



4-5.2 Bosonization of the effective Hamiltonian 

The interaction is treated by bosonization of the chiral Fermions. After bosoniza- 
tion, the Hamiltonian H^s is quadratic in the boson fields 0. The price we pay, 
however, is that the effective action S^s is not quadratic in the boson fields. 

We introduce boson fields 0Lj and 0F{j for left moving and right moving par- 
ticles by[3 

= -^e-^'^-^ ^^,. = -^e^^«^ j = l,...iV,h. (294) 



Here N^h is the total number of channels (including the spin degeneracy), and 
fjj = r/j is a Majorana fermion, = 26ij. Since the Majorana fermions 

do not enter into the Hamiltonian, their average is given by 

{%flj{Ti),7]i{T2))q = SijSgn (Ti - T2) . (295) 



The boson field 0Lj and 0Rj satisfy the commutation relations 

[(pLjix), (pLiiy)] = -iiT5ijSgn{x - y), (296) 
[ipRj{x),(fRi{y)]=iTc5ijSgn{x - y), (297) 
[(fRj{x),(fLiiy)] = -irt5ij. (298) 

The high energy cut-off 1/ A is of the order of 1/Ai? and is chosen consistently 
with the high-energy cut-off in the bosonic correlation functions, see below 
Eq. (303). 



These bosonic fields are related to the fields of Eq. (208) as 6j{x) = {(pRj — 
Lpij)/y/2T^, (pj{x) = {(pRj + ipLj) / The linear scale A is related to the cut-off D 
there by A = D/vp 
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The products 'ipij'^Lj ipj^jipj^j read 

Id ^ ^ Id 

■i^Ljix)'^Ljix)-= ^g^VLjix), ■i^Rj{x)i!Rj{x): = ——^bA^)- (299) 

It is easy to check, using Eqs. (296) - (298), that the fermionic fields (294) obey 
the standard anti-commutation relations. The commutation relation (298) 
guarantees the anticommutation relation jV'Lj, "^ij.j} = 0, | V'lj j} = 
and the correct commutation relations between the operator of the number of 
particles in the dot, see Eqs. (254) and (299), 



n ■■ 



i=i-oo 

1 iVch 



TT ["PrA^ = 0) + <^Lj(x = 0)] 



27r 



(300) 



and the fermionic operators. 



In terms of the bosonic fields, the effective Hamiltonian (267) acquires the 
form 



47r Y 

J — oo 

E, 



dx 



'd0L,j 

dx 



'dip 



R,j 



dx 



+ 



(301) 



27rAr 



The bosonized form of the Hamiltonian (293) is obtained by the replacement 
Af^ (Af-m). 

As the Hamiltonian (301) is quadratic in the bosonic fields, it can be easily 
diagonalized. One can immediately notice [38,39] that the eigenenergies of 
(301) do not depend on the gate voltage Af because the latter can be removed 



The absence of a term proportional to the bosonic fields at x = — oo in Eq. (300) 
follows from the order of limits taken with respect the range of the interaction and 
the system size: First, the interaction is defined in terms of the charge of the entire 
length of the channel for x < 0, and only then the length of the channel is sent to 
infinity. 
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by the transformation ^L{R),j VL{B),j " T^-^/Nch- On the other hand, the 
effective action (267) becomes a non-hnear functional of the bosonic operators 



l/T l/T 



SeS = X] / / '^^2 Lij (ri - T2) f]i{Ti)T)j{T2) 



V 



i'PLj{r2) _|_ Q-i'PRj{r2) 



(hereinafter we use the convention (265) for the operators taken at x = 0). 
Our general strategy will be to expand the quantity of interest in powers of the 
effective action, and then to average over the quadratic Hamiltonian (301). All 
the relevant averages have the form exp(a), where a is a linear combination 
of the boson fields ipij and (pjij, j = 1, . . . , N^h- Such an average can be found 
according the rule 



(exp(a))q = exp 



(302) 



The corresponding bosonic correlation functions can be found from the Hamil- 
tonian (301), see e.g. Ref. [70]. For readers wishing to reproduce the calcula- 
tion, we notice only that the commutation relation (298) should be taken into 
account to obtain the Green functions with correct analytic properties. The 
result is 



(A ttT \ / 1 \ 

^|sin7rTVi-r.)|j+^+(^-~^)'' 

{Tr(f>Ri{Ti)(f!Rj{T2)) = {Tr^Li{Tl)^Lj{r2)), 

%1T 1 

{TripLi{,n)ipRj{T2)) = 5^j—sg■^l{Tl - Tg) - — ln/(r2 - Ti) - A. (303) 

^ ^Vch 

Here all of the bosonic operators are taken at the origin x = 0. One should 
take the limit A +00, B +00 at the end of the calculation of each 
physical propagator.p^ 

The function ln/(r) has the integral representation 

XEcNche^ 7 _^ sinh[(ir + 7rx/^cA^ch)vrT] 

In / r = In dxe In [\ ' \ \ 304 

nvp J smh(7r^i x/Eci\ch) 



The precise meaning of these quantities is NchA/2 = J2ij{V'L,i{T)'PL,j{T))q, 
(iVch - 1)B + A/2 = T.^{HAr)Mr))q 
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where C ^ 0.577 is the Euler constant. We note that /(r) is analytic for 
Imr < 0. For T -C Ec, to a good approximation, we set 

^^^^^i^^sin[(r-iio)7rr]' = ^B^A^" 



This equation is vahd both for r = and for r ^ 1/A^ch-E'c- K also provides 
the correct upper cut-off for the arising logarithmic divergences. The approxi- 
mation (305) fails when contributions from times r ~ l/N^hEc are important, 
as is the case for spinless fermions with Net = 1- 

Using Eqs. (294), (295), and (303) one can find all the relevant correlation 
functions of the fcrmionic operators. We will discuss some of them here in 
order to demonstrate the analytic structure of the perturbation theory. 

First of all, from Eq. (303) we immediately find that any Green's function 
which involves only left (L) or only right (R) moving fermions is not sensitive 
to the interaction, i.e., 

TlT 

{2T:vF){T^i)L,i{n)'^m{r2))q = '^'J sin[7rr(Ti -T2)] ' 



and all the averages are given by the Wick theorem. All such Green functions 
may have only pole type singularities in the plane of the complex variable Tj. 
The insensitivity of such type of the Green functions to the interaction can be 
understood from the chirality of the fermions and the causality principle. 

Let us now consider a more complicated average involving left movers together 
with right movers. 



{2nvFf{Tri^LA^SLAr2)i^RA^^)'^RAn)),= (307) 

ttT ttT / /(r4-ri)/(r3-r2) y/^-' 

sin[7rr(Ti - T2)] sin[7rr(T3 - T4)] \/(t3 - Ti)/(t4 - T2) J 

We fix T2,T3,T4 to be real numbers and consider the correlator (307) as a 
function of the complex parameter ti . It has a pole singularity at Ti = T2 and 
branch cuts along the lines Reri = t^, Imri < 0, and Reri = T4, Imri < 0, 
see Fig. 23. It is readily seen, that the contribution to an integral over ri from 
the pole does not depend on the interaction (it does no longer involve the 
function /): the fractional-power factor in Eq. (307) turns into unity at the 
pole. Thus it is the branch cut contribution that comes from the interaction. 
The same is true for any higher-order correlation function. 

Imagine now that we expand an observable quantity in power series in L^j 
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Fig. 23. Analytic structure of the correlator (307) as a function of ti. The dot 
indicates the pole singularity at ti = T2; the thick lines the branch cuts along the 
lines Reri = rs, Imri < and Reri = t^, Imri < 0. 

and account in each order of perturbation theory only for the pole contri- 
butions in the integrals over the intermediate times ri,r2, . . .. The pole con- 
tribution in each order does not depend on the interaction, and constitutes 
just a Green function for non-interacting fermions. Such contributions can be 
easily summed up, resulting in a geometric series in L.\^ Since a geometric 
series of the matrix L is precisely the scattering matrix S, cf. Eq. (270), we 
thus arrive at the following simple result: in order to take into account the 
pole contributions in all subsequent orders of perturbation theory, one has to 
simply replace L in the first non-vanishing contribution by 



(Here Un is the fermionic Matsubara frequency.) 

Consider now a correlation function which includes a branch cut in the leading 
order {e.g., as function of the times ri, T2, r^, T4 in second order in the effective 
action). Then, in each of the subsequent higher orders we select the contri- 
butions including only one branch cut, all other singularities in these terms 
being single poles. Since such a contribution is factorised into a branch cut fac- 
tor involving two times, and a product of Green functions for non-interacting 
fermions for the remaining times, we can again sum up all such contributions 
as a geometric series. As before, the result is that one has to take the first non- 
vanishing term in an expansion in powers of L that contains a single branch 
cut and then make the replacement (308) to account for all higher orders in 
L with the same branch cut singularity. In such a way, we take into account 
all the orders of perturbation theory in L that contain only one branch cut 

For the calculation of the ground state energy, the geometric series appears only 
when one calculates the thermodynamic potential (273) of the entire system, in- 
cluding the contribution Q.^ from the closed dot. 
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singularity in each order. 

The same analytic structure occurs in all higher orders of perturbation the- 
ory: poles give a "reducible" average, whereas the "irreducible" part of the 
propagator has only branch cut singularities. This leads to a very important 
statement: the perturbation theory in the action (267) can he constructed as 
an expansion in powers of the scattering matrix S rather than of the kernel 
L. In this expansion all the reducible averages (poles) are included automat- 
ically in all the orders of the perturbation theory, whereas the branch cuts 
are considered perturbatively; each order of of the perturbation theory in S is 
characterized by the number of branch cuts it includes.P^ 

The correlators (306) and Eq. (307) do not depend explicitly on the gate volt- 
age A/". One can see from Eq. (303), that a correlator having such a dependence 
must include at least creation operators of left (right) moving particles and 
annihilation operator for right (left) moving particles in all channels. For a 
total number of A^ch channels, the lowest order A/'-dependent correlator is 



X 



^ fjrf - r/) 



/(O) 



The expression for the correlator of an arbitrary number of the fermionic oper- 
ators as well as other necessary correlation functions are given in Appendix G. 

To conclude this subsubsection, let us prove the assumption of Sec. 4.4 about 
the Fermi liquid behavior of the one-channel system A'^ch = 1 at low energies. 
To this end we substitute Eq. (305) into Eqs. (307) and (309), and discover 
that the Green functions now have only the pole- type singularities, which is the 
hallmark of the Fermi-liquid. Moreover, one finds from Eq. (309) at t ^ 1/Ec 



which is exactly a Green function of one channel system with unpenetrable 



Regarding this statement, we point to the fact that for noninteracting particles, 
the ground state energy, Eq. (Ill) and the conductance, Eq. (103), are of second 
order in the scattering matrix S. In the effective action theory, this result follows 
from the summation of the pole contributions to all orders in L, which produces 
a single term proportional to S'^. The interaction corrections, however, can be of 
higher order in S, depending on the number of branch cut singularities that are 
included in the perturbation theory. 
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barrier and the scattering phase S — n {J\f +1/2), see also Eq. (241). 



4.6 Results for the differential capacitance 

As we have aheady seen in the discussion of weak tunnehng, the existence 
of periodic oscillations of measurable quantities with the gate voltage A/" is a 
signature of Coulomb blockade. In this section we consider such oscillations, 
which occur due to the incremental charging of the dot, in the case of the 
differential capacitance of a partially open dot. 

First we derive the relation between the AT-dependent thermodynamic poten- 
tial and the scattering matrix S for an interacting system, see Eqs. (313)-(316) 
and (331)-(332) below. These formulae replace the known results for the non- 
interacting electrons, see Eqs. (111)-(114) of Section 2.4. 

Using the general result, we then find the ensemble- averaged differential ca- 
pacitance. The average capacitance oscillates periodically with A/" as long as 
some reflection in the channels connecting the dot with leads exists. The am- 
plitude of oscillations, however, decays with the number of channels as |r|^'='', 
where r is the reflection amplitude in a channel, and the number of channels 
A^ch accounts also for the spin degree of freedom. The detailed formulae for the 
oscillations of the average capacitance are given below in Eqs. (122), (321), 
and (323). 

In the absence of reflection in the channels, the average differential capacitance 
does not oscillate with J\f. However, mesoscopic fluctuations of the capacitance 
still bear a signature of charge quantization. The general formulae Eqs. (313)- 
(316) and (331)-(332) allow us to examine the two-point correlation function 



The flrst part, Kq{s), originates from the discreteness of charge, and can be 
represented as the product of an oscillatory function with unit period and an 
envelope function Sjv^^(sA/£'c) which decays on a scale |s| ~ Ej A, 



K{Ui-J\f2) 



(^Cdiff(A/'i)'^gdiff(.A/'2)) 

(J2 



(310) 



and show that it consists of two additive parts. 



K{s) = Kq{s) + K,{s). 




(311) 
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The explicit results for the special cases of a junction with one propagat- 
ing mode of spinless fermions and real spin 1/2 electrons are given below in 
Eqs. (319) and (321), respectively. The second part, K^^s), comes from the 
mesoscopic fluctuations of the capacitance as a function of the chemical po- 
tential (which itself depends on J\f) . This part is identical with the one that 
was found earlier in Ref. [144] within the RPA approximation. 




The derivation of K^^s) within the general formalism is also outlined at the 
end of this Section. 

With the increase of the number of channels Nch connecting the dot and 
leads, the oscillatory part of the correlation function quickly decays. It is 
experimentally detectable only for a small number of channels. 

The oscillatory part Kq{s) of the capacitance autocorrelation function 
As wc have seen in the previous subsection, the A/'-dcpcndent contribution 
to the thermodynamic potential appears in A/'chth order perturbation theory 
{Nch is the total number of channels): 

6Q{Af) = ^ , {[Se^f''),. (313) 



Substituting Eq. (267) into Eq. (313) and using Eq. (309), we obtain 



SQ{Af) = -2 (-1)^^'' Rc er'^^^T 



X 



l/T 

J det[S] 



iVch 

n 



7(rf-r/ 
. 27rto/(0) 



(314) 



where the function /(r) is defined by Eq. (304). The determinant of the scat- 
tering matrix in the imaginary time domain is defined by 



det[S] = 



'52i(r^-rf) 



'^iVehl(Ti^,h 



^1 ^iVr> 



Sin, 



'S'ATehAfchl'^JVeh 



R 



(315) 



The kernel L was substituted with the scattering matrix S in Eq. (315), in 
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accordance with the previous subsubsection, and 



S{t) =T J2 e- 

a;„=7rT(2n+l) 



S{iu;n)0{u}n) - S\iuJn)0{-UJn) 



The latter formula can be also written in a form of the Lehmann representation 

°° rri 

Sir)-^Jdt[S(t)^SHt)] (316) 

— OO 



Equations (314) and (315) express the gate dependent interaction to the ther- 
modynamic potential in a most general form. It is valid for any scattering 
matrix of the dot. In what follows we apply those equations for studying the 
different limiting cases. 

The oscillatory pari of differential capacitance and Kq{s) for the one-channel 
geometry (spinless fermions). - In this case the scattering matrix is formed 
by its only element S{t). Using the fact that the function /(r) is analytic in 
the lower complex semiplane, see Eq. (304), we find 

E 




With the help of Eq. (304), one obtains for T <^ Ec 



E^ 
X Re 



oo / oo 

j S{t)exp I- J dxe-^'ln 



X + tEc/ 7T 



dt 



X 



■ (317) 



One can use Eq. (317) to analyze the statistics of the differential capacitance 
(274). In the approximation (271), one obtains [39] 



c 



2e^|r| cos27rA/', 



(318) 



where (...) stands for the ensemble averaging and r is the reflection amplitude 
of the contact between the dot and the lead. 

The statistics of the mesoscopic fluctuations of the differential capacitance for 
reflectionless contacts was analyzed in Ref. [70] with the help of Eqs. (115), 
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(122), and (317). The result is 



Kq{s) ^ ^-cos2ns. (319) 



The correlation function for the capacitances at different magnetic field can 
be found in Rcf. [70]. It is noteworthy that the characteristic amplitude of the 
fluctuations is in agreement with the estimates of Sec. 4.4. 

The oscillatory part of differential capacitance and Kq{s) in the one-channel 
geometry ( electrons with spin). - In the absence of spin orbit interactions, the 
2x2 scattering matrix is given by S{t) = diag{S'°(t), S"{t)}. In this case, the 
main contribution comes from times larger than E^^, so that one can use the 
approximation (305) for the function /(r) in Eq. (314). With the help of the 
Lehmann representation (316), one finds 



27r2 



Re 



^-i2nM J dtidt2S{ti)S{t2)F{ti,t2), 



(320) 



F(ti,t2) = - 



1/2 



1/T 
X / dT3 



[sinh(7rr(ti + to)) sinh(7rT(t2 + to)] 

ttT 

[sinh7rr(ii +to- ir^) sinh7rr(i2 + io + ir^)]^'^ ' 



TT 



Similarly to the model of the spinless fermions, Eq. (320) can be applied 
to study of both the averaged differential capacitance and the mesoscopic 
fluctuations. Using Eq. (271) and Eq. (115), one obtains 



-7 — — — |r| cos27ryV m — 

C TT \T 



(321) 



Kq{s) = ^ cos(27r.)|^ In^ ( \^ In f §U 4e 



The result for the averaged capacitance was first obtained in Ref. [39]. The 
correlation function of the mesoscopic fluctuations in magnetic field can be 
found in Ref. [70]. The origin of the logarithmic divergence in Eq. (321) and 
the cure for it is discussed in Sec. 4.2. 



Capacitance oscillations in the case of multichannel geometry with reflection 
in the junction. - For the case of an ideal point contact with many channels 
between the dot and the lead, the periodic-in-jV capacitance fluctuations given 
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by Eq. (314) will be dominated by non-periodic fluctuations, see the following 
Subsection. Periodic capacitance fluctuations, however, remain important if 
the point contacts are non-ideal. In that case the scattering matrix contains a 
non- fluctuating direct component, given by the reflection matrix of the point 
contact, cf. Eq. (121), 

Sij{t) — rijd{t) + fluctuating part. (322) 

Substitution of Eq. (322) into Eq. (314) gives a periodic contribution to the 
thermodynamic potential, 

dn{M) = ^c^,. In (1^) Re det r, (323) 

where the numerical coefficient cat^^ depends weakly on the number A^ch of 
channels; C2 = 1/4, cg = (9/8)7r-=^/22i/3r(i/3)r(7/g)_ 

Equation (323) demonstrates that the dependence of the amplitude of the 
Coulomb blockade oscillations on the contact's conductance G is not universal; 

more detailed characteristics of the contact are important. The determinant 
det r in Eq. (323) can be evaluated in terms of the reflection eigenvalues Rj of 

1/2 4- 

the contact, | detr| = Yij R/ (here Rj are defined as the eigenvalues of rr^). 
The distribution of Rj is known for a number of relevant models, including the 
model of a short disordered wire [14] . With the help of that distribution, one 
can find from Eq. (323) the amplitude of the oscillations of Q, and of the dif- 
ferential capacitance, 5Cdis{J^) ~ ln(£'c/T) exp(-27r^^G/e^) cos(27rA^). This 
exponential dependence, and the more general formula expressing the ampli- 
tude of oscillations in terms of Yij Rj (but both witout the logarithmically 
divergent prefactor) were found in Ref . [138] . There a different method, based 
on an extension of the instanton technique of Ref. [145] to the case of an 
arbitrary (though energy independent) S- matrix of the contact was used. 
The same results were reproduced within a non-linear cr-model for interacting 
systems in Ref. [146]. 

Non-periodic interaction corrections to the capacitance, K^ls). - In addition 
to the periodic-in-A/" mcsoscopic fiuctuations of the capacitance, resulting from 
charge discreteness, there arc also non-periodic mcsoscopic fluctuations of the 
capacitance. These non-periodic fluctuations arise through mcsoscopic fluctua- 
tions of the thermodynamic potential Q as a function of the chemical potential 
II in combination with the jV-dependence of the chemical potential /i. In the 
absence of direct reflection in the contacts and for temperatures T ^ NchA, 
the non-periodic A/'-dcpcndcnt fluctuations dominate over the periodic ones 
that occur in higher (A'ch) order perturbation theory in the effective action. 
[Higher orders in perturbation theory are small by the quantity (A^ch^/TT^T).] 
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Mesoscopic fluctuations of tlie capacitance Cdifr were studied in a series of 
papers in a mean-field treatment of the interaction [94,147,148,144,149], with 
the result 

^ 1 + \ (324) 



Cdiff 2E,iy{fi;Ty 



Here z/(e; T) is the non-interacting density of states of the dot at temperature 
T, cf. Eq. (113), 



fri^) = 1/(1 + exp(£:/T)) being the Fermi distribution function, and the J\f- 
dependence of fi follows from the self-consistency condition 



C, 



diff 



dfl 

dAf ^ ~Cu{fi;T) 



2Er 



(325) 



The mesoscopic fluctuations of the capacitance (324) are caused by those of 
the (non-interacting) density of states i^ie). We now show that these capac- 
itance fluctuations can be reproduced within the effective action approach. 
The advantage of the effective action theory, however, is that it allows us to 
calculate systematic interaction corrections to the capacitance that go beyond 
mean-field theory. 

Formally, the A/'-dependence of the chemical potential fi in the effective action 
approach arises, because one has to impose a canonical constraint on the 
ensemble, instead of a grand canonical one, as was done in Sec. 4.5.1. Details 
of how this can be done are given in Appendix F. Here we only mention the 
result, which is that for temperatures T ^ Ai'chA or for channel numbers 
Nch ^ 1, where the fluctuations of the density of states u^e; T) are small, one 
finds 

^ = -1^. (326) 



Note that this is the same A/'-dependence as in the self-consistent theory, Eq. 
(325). In the limit Ec A that we consider here, Eq. (326) simplifies to 
/i ^ -AM. 

The difference between Eqs. (326) and (325) is in the fluctuations of the density 
of states, which are small when T ^ N^h^ or N^h ^ 1- Fluctuations of i^(e) can be 
included systematically in the effective action formalism, see App. F 
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By Eq. (273), the thermodynamic potential Q consists of two parts, 

O = - T In Tr e"^««/'^T^e~^«ff , (273) 

where is the thermodynamic potential of the closed dot without interac- 
tions and the second term accounts for the difference between the thermody- 
namic potentials of an open dot with interactions and a closed dot without. 
The thermodynamic potential Qd of the closed dot without interactions can 
be written as an integral over the density of states t'ciosed(£) of the closed dot, 

^D^-T j devao.ed{e) In (l + e"^/^) . (327) 



The density of states i^ciosed(£) is a sum of delta functions, and thus differs 
from the density of states vi^e) of the dot when it is coupled to the leads, 
see Eq. (113), which is a continuous function of energy. It was the latter 
density of states that determined the thermodynamic potential (111) of the 
open quantum dot in the absence of interactions, and that entered into the 
mean-field description, see Eqs. (324)-(325) above. Using the equahty 



inclosed (s) 



d_J_ 
de 2m 



Tr In 



l+5(£-iO+) 



(328) 



together with unitarity of 5", we can write t'ciosed(^) as 



1 d 

Inclosed (£) = y{e) - -^Im ln[l + S{e + iO+)]. 

TT 0£ 



(329) 



Hence, we can represent f)^, in turn, as 



Q^ = Qib + - / rf£/^(e)ImTr ln[l + 5(£ + iO+)], (330) 

TT J 

— oo 

where Qib is the thermodynamic potential of the dot coupled to the leads in 
the absence of interactions, see Eq. (111). (The subscript "lb" has been writ- 
ten in analogy with the non-interacting contribution Gib to the two-terminal 
conductance, see Subsection 4.7 below.) 

The second term in Eq. (273) can be calculated by expansion in the effective 
action 5efr, followed by an average over the fermionic operators using Eqs. 
(306) and (307). As explained below Eq. (307), one can distinguish two types 
of contributions: the "pole contributions", which are already present in the 
absence of interactions, and the "branch-cut contributions" , which represent 
the corrections to Q due to the electron-electron interactions in the dot. Since, 
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in the absence of interactions, the thermodynamic potential of the dot is given 
by f2ib, the pole contributions to the second term in Eq. (273) and the second 
term of Eq. (330) must cancel to all orders in the effective action Ses- Hence, 
the thermodynamic potential is given by 

O = Oib + flee, (331) 



where the interaction correction Qge represents the "branch-cut contributions" 
to — TlnTr e~^°ff/"^Tre^'^°« . (The calculation of fipc proceeds along the same 
lines as the calculation of the interaction correction to the two-terminal con- 
ductance; Details of the latter calculation can be found in Sec. 4.7 and Ap- 
pendix H.) To second order in the action Ses {i.e., up to one branch cut 
contribution), we find 



^ee — 



. TT 7 . . 7 , TrS{h)S^{h) 
2 A^eh ■/ J sinh[7rr(ii + s)] sinh[7rT(i2 + s)] 

tQ 



X 



f sinh[7rT(g - tp)] sinh[7rT(s + ti + + tp)] ] 
\ sinh[7rT(ti + to)] sinh[7rT(t2 + to)] J 



(332) 



In order to find the differential capacitance (274) from Eqs. (331), (113), and 
(332), we need to include the dependence on the chemical /x, which is done 
via the relation 

S{^Ji,t) ^ S{Q,t)e-'i'\ (333) 



Hence, with help of Eq. (326) we find for > A 



C'diff — C, 



diflfjb 



c 



C'hv 



diff,( 



c 



diflf — ^diff,lb + C'diflf.ee 



(334) 



A 



A^ 



oo oo 



dt J dt'Tr[S^{-t')S{t) - {S^{-t')S{t))] 





X 



oo oo oo 



sin-^ [ dt f dt' f dsTrS\-t')S{t) 

Nrh J J J 



X 



X 



'KT{t - t' ) iAAfit-t') 

sinh[7rr(i - 1')] 

A2 
— sm- 

ch 

7r^r^(t - t 

sinh[7rr(t + s)] sinh[7rr(t' + s)] 

f sinh[7rT(s - to)] sinh[7rr(s + t + t' + to)] 1 ^ 
\ sinh[7rT(t + to)] sinh[7rT(t' + to)] J 



to 
2rr2u +r\2 



(335) 



(336) 
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In the first term, CdifT,ib, we have separated the ensemble average and the 
mesoscopic fiuctuations. The interaction correction Cdifr,ee goes beyond mean- 
field theory, and reveals the onset of Coulomb blockade. Its ensemble average 
is zero, because of the presence of the factor (t — t'Y in the integrand. 

For temperatures T ^ A^'chA, the capacitance fiuctuations can be computed 
from Eq. (334). After averaging over S{t) with the help of Eq. (122), one finds 



var (Cdiff/C) 



TC 



123 



A 



where c is a numerical factor that takes the value c ~ 3.42 for Nch oo. 
The first term between brackets {. . .} represents the fiuctuations of Cdifr,ib 
(calculated to leading order order in Nch^/T only), while the second term 
is the leading interaction correction. We conclude that for N^h ^ 1 and/or 
T 3> iVchA the mesoscopic capacitance fiuctuations are dominated by the fiuc- 
tuations of Cdiff.ib, corresponding to the mean-field theory of Ref. [94,147,148], 
the interaction corrections being a factor ~ N^(Nch^/7^'^T)^ smaller. Note 
that at ttT ~ A'chA/Tr, which is where the crossover to the low-temperature 
limit of the theory is expected, the contribution of Ccc to the capacitance 
fiuctuations is still smaller by a factor A^^ch- 

We close this section with the correlators of the (mean-field) capacitance at 
two different gate voltages Afi and A/2 = A/i + s for A^ch ^ 1- From Eq. (335) 
and (122), one finds for T > A^'^hA, 



3var (Cdiff/C) 
sinh^(As/2r) 



A , /A 

— -s coth — -s 
2T \2T 



The fiuctuations of CdifT,ib for A'ch ^ 1 in the low temperature limit ttT <^ 
AT^hA/vr read [147,144] 



' A2 
2E,^ 
4 / 2A ^2 



{6u{AfA)Su({Af + s)A)) 
1 - 4 (2s/A^ch)' 



(338) 



l + 4(2s/Ar,h) 



2 • 



where the correlator of the density of states was taken from Ref. [150] .[^ 



To reproduce this result from Eq. (335) it is not sufficient to use a Gaussian 
average witli the correlator (122), as non-Gaussian contributions are also important 
for times of order t ~ l/N^h^- 
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(The argument of the density of states is A/'A, since fliictTiations in the self- 
consistency relation (325) can be neglected for this correlator for Net ^ 1-) 



4-7 Conductance through a dot with two almost reflectionless junctions 



This subsection is devoted to the effects of interaction on the two-terminal 
conductance of a dot with junctions to the leads that are either completely 
reflectionless, or have only a small reflection coefficient. The statistics of the 
two-terminal conductance in the absence of the interactions is reviewed in Ref. 
[14]. The effect of interaction on this quantity was analyzed in Ref. [151] and 
we follow this reference in this subsection. 

The Kubo formula for the conductance G is given by Eqs. (281), (279), and 
(277). Similarly to the consideration of the capacitance, wc expand the current- 
current correlation function (281) in powers of the effective action Scs and 
calculate the arising fcrmionic correlation functions with the help of Eqs. (306) 
- (309) and (G.3) - (G.4). 

The lowest non-vanishing contribution to the current-current correlation func- 
tion n(r), and hence to the two-terminal conductance, is of the second order 
in the action, 

n(r) = (/(r)/(0)), + l{I{r)mS',^), - ^(/(r)/(0)),(5,^ff), + . . . . (339) 



We substitute into Eq. (339) the explicit form of the current operators (277) 
and the effective action (267). After calculation of the averages of the fcrmionic 
operators with the help of Eqs. (306) - (309) and (G.3) - (G.4), and with the 
approximation (305), we find 



G — Gib 

^2 



G 



lb 



27ih 



Gcc + G. 



OSC) 
OO 



ch 



dti / dt2 



Gpp — 



X 



2'Kh VTT 

ds 



1 

— sm 
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ch 



(ti-t2)7rTTr5t(-ti)A^(t2)A 
sinh[(ti — t2)7TT] 

X) 

J dtij dt2TrS\-ti)AS{t2)A 

2rp2 





OO 



(2s + ^2 + h)7:'^T 



to 



X 



sinh[(s + ti)7rT] sinh[(s + t2)7rT] 

'sinh[(ti -f- ^2 + -s -|- to)7rT] sinh[(s — to)7iT] \ 
^ sinh[(ii + to)7rT] sinh[(t2 + to)7rT] ) 



(340) 
(341) 



(342) 
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The charging energy enters into the expressions only through the time to, see 
Eq. (305), which has the meaning of the characteristic time of the classical 
recharging of the dot (i?C-time). Notice, that the final results are written in 
terms of the scattering matrix S of the dot in the absence of the interactions, 
see Sec. 4.5.2. In this order of perturbation theory, the term Gosc oscillating 
with gate voltage is present only for the case N2 = Ni — 1 (spinless fermions, 
and both contacts carry a single channel each): 



<^osc = -r^ / dti dt2 / ds 1 = 
I i ^sinh[(ti + to)7rr] sinh[(i2 + h)T:T\ 



X 



^sinh[(i2 + + s + to)7rr] sinh[(s - io)7rr] 
X Ree-^'^^-^ [S^i{ti)S22{t2) + Sr2{ti)S2i{t2)] ■ (343) 

The derivation of Eqs. (340) - (343) is relegated to Appendix H. 

Equation (341) coincides with the Landauer formula for the conductance of 
non- interacting electrons (104). Equation (342) is the leading interaction cor- 
rection at temperatures much smaller than the charging energy. The power- 
law form-factor (the last line in this formula) describes the non-Fermi liquid 
behavior of the dot when Nch > 1- Finally, the Coulomb blockade hke contri- 
bution (342) for iVi = A^2 = 1 is somewhat similar to the oscillation of the 
capacitance, Eq. (314). 

To gain intuition about the structure of the interaction correction, we start 
with the simplest example of instantaneous reflection from the contact, when 
S{t) is given by Eq. (322). Neglecting the fluctuating part of S and substituting 
Eq. (322) into Eq. (341), we find 
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The first term in Eq. (344) is nothing but the "classical" conductance of the 
two refiectionless junctions connected in series, and the second term comes 
from the correction to the transmission amplitudes due to the finite backscat- 
tering in the contacts. 



Substitution of Eq. (322) into Eq. (342) yields 
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where cat^j^ is a numerical coefficient that depends only weakly on Ach, ranging 
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from C4 5.32 to Coo = 4. Adding Eqs. (344) and (345), we obtain 
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Comparing Eq. (346) with Eq. (344), we see that the "bare" reflection amph- 
tude entering in the Landauer formula is renormahzed by the interaction. The 
renormahzation is temperature-dependent, manifesting the non-Fermi hquid 
behavior of the open dot on the energy scale Ec^tZ A. For Ni ^ N2 ^ 2 
(electrons with spin, and both contacts carry one channel each) the result 
(346) has been obtained previously in Ref. [40]. 

The above example (322) of reflection at the moment t — corresponds to 
case where the electron is backscattering instantaneously, without entering the 
dot. Now we are going to estimate the interaction-induced renormahzation of 
the amplitude for a backscattering event occurring somewhere within the dot. 
In the language of the electron trajectories, such an event corresponds to a 
trajectory returning to the contact at some time ts > 0, 



It is easy to see that the noninteracting contribution (344) remains intact, 
whereas the interaction contribution depends signiflcantly on ts'- 



(for brevity, we omitted here the numerical coefficients of the order of unity). 
We see that for short trajectories, with return times ts smaller than the RC- 
time for the open dot, see Eq. (305), the renormahzation of the scattering 
amplitude is complete, i.e., the same as for the direct reflection from the con- 
tact we studied before. On the other hand, the effect of interaction on the 
scattering events occurring "deep inside" the dot, ts » l/T" is exponentially 
weak. As long as the temperature is sufficiently high, T » ANch, and the 
scattering from the dot is ergodic, the majority of the scattering events is not 
renormahzed by the interaction. Therefore, the Landauer part of the conduc- 
tance (344) remains intact at such temperatures. The interaction correction 
affects only a small portion of scattering events, and thus only slightly affects 



Sij{t) = rijS{t - ts), ts > 0. 



(347) 
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the statistics of the two terminal conductance. Note that the characteristic 
energy A^ch^ that appears here is nothing but the inverse escape time for a 
quasiparticle from the dot into one of the leads. 

To study the statistics of the two-terminal conductance, one should use the sta- 
tistical properties of the scattering matrix (115), and Eqs. (340) - (342). The 
result substantially depends on whether we are dealing with spinless fermions, 
denoted by the spin- degeneracy factor = 1, or electrons with spin, denoted 
by Qs — 2. In the latter case, the results are written in terms of the number of 
orbital channels N^^ = N^h/gs (and, similarly, A^^" = Ni/qs, = N2/gs)- Be- 
low we concentrate on the case of the refiectionless contacts; all backscattering 
events occur within the dot. For the averaged conductance, we find [151] 



2nh 



(349) 



+ 1- 



1 + 



where P — 1 (2,4) for the orthogonal (unitary, symplectic) ensemble, and cjv^j^ 
is an Ai'ch-dependent numerical constant, ranging from C4 ^ 4.77 to Coo = 7r^/6. 
Equation (349) is valid for an arbitrary number of channels except the case of 
gs = 1, Ni = N2 = I. For the latter case, one finds 
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where c ^ 2.81. 

Similar to Eq. (344), here the first term in brackets is the classical resistance 
of two point contacts connected in series. The second term is the familiar 
weak localization correction for non-interacting electrons [this term replaces 
the correction to the Kirchoff law coming from the backscattering in the junc- 
tions, see Eq. (344)]. Finally the last term is the high-temperature expansion 
of the interaction correction. This correction comes from the renormalization 
of the amplitudes of the backscattering events occurring in the dot within 
"time horizon" ts ~ ^/T, cf. Eq. (348). At high temperature, the interaction 
correction is less significant than the correction from the weak localization. 
Notice, however, that at ttT = ANch/'n', which is the lower limit of the ap- 
plicability of the theory, the interaction correction is smaller than the weak 
localization correction only by parameter 1/A^ch, and may be still important 
for the temperature dependence of magnetoresistance. 

One notices also that for the unitary ensemble, P = 2, both the weak localiza- 
tion and interaction corrections to the average conductance vanish, in contrast 
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with the usual situation in disordered metals [18]. The reason is that the in- 
teraction enhances any electron scattering, see Eq. (348), no matter whether 
it is scattering from one lead to the other (1 2) or backscattering into the 
same lead (1 <-> 1 and 2 <-> 2). For the unitary ensemble, the scattering 1 <-> 1 
and 1 <-> 2 occur with equal probabilities and therefore neither reflection nor 
transmission arc enhanced in average. On the other hand for the orthogonal 
ensemble, backscattering into the same lead is more probable due to the weak 
localization, and the interaction further increases the backscattering probabil- 
ity. 

In a similar manner, the mesoscopic fluctuations are also affected by the 
electron-electron interaction. We find from Eqs. (340) - (343), and Eq. (115), 
in the regime T ^ AA^^ch, 



var G: 



i ( ^ 



6T iV°h gsTT^T^ 
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with c' ~ 6.49 for A^'ch ^ 1- The first term in Eq. (351) represents the high 
temperature asymptote of the mesoscopic fluctuations of non-interacting elec- 
trons [152]. The RMT index /3 = 1 (2) in the presence (absence) of time- 
reversal symmetry; /9 = 4 for the symplectic ensemble, which corresponds to 
the presence of both strong spin-orbit scattering and time-reversal symmetry. 
For the case Ni = N2 — 1, Qs — there is an explicit jV-dependence of the 
conductance fluctuations. 




647r2r2 



{l + cos[27r(AA-AA')]}ln' 



(352) 
T 



where is the charging energy of the dot and C Ri 0.577 is the Euler constant. 
In this order of perturbation theory, such an explicit A/'-dcpcndence, which 
reflects the discreteness of charge, exists only for spinless fermions, but is 
absent for particles with spin. For spin 1/2 electrons, effects of the discreteness 
of charge appear only in the A^ch-th order in perturbation theory in S, similarly 
to Eq. (314). For a large number of reflectionless channels, the characteristic 
amplitude of mesoscopic fluctuations goes down very rapidly as (1/A^ch)^''''- 

The situation here is similar to that of the differential capacitance, see Sec. 
4.6, where we found that for A'ch > 2 the average and variance were domi- 
nated by the mesoscopic fluctuations of the non-interacting density of states, 
while the interaction corrections were smaller by a factor N~^^{Nch^/T^'^T)^ 
for T 3> iVchA. The periodic-in-jV Coulomb blockade effects become domi- 
nant for Nch < 2, thus limiting the range of validity of the non- interacting (or 
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mean-field) theories to larger values of Nch- This limitation is more relevant 
for the capacitance, where the case N^t = 2 can be realized experimentally for 
electrons with spin and a single-mode point contact, than for the conductance, 
were the case Nch = 2 is mainly of academic interest, as it requires spinless 
electrons. 



4-8 Conductance in the strongly asymmetric setup 



We consider here a dot with two leads, one of which is connected to the dot 
by a tunneling contact with a small transmission coefficient, while the other is 
coupled to the dot via an ideal (no reflection) or almost ideal (small r) point 
contact. The first lead serves as a tunneling probe of the "open" quantum dot. 
See the discussion in Sees. 2.4 and 4.5.1. 

Inelastic tunneling. - We begin with the analysis of the inelastic contribution 
to the tunneling conductance. According to Sec. 4.5.2, the order of perturba- 
tion theory is characterized by the number of the branch cuts the correspond- 
ing fermionic correlator includes. In the leading order, we will allow only one 
branch cut in the complex r plane and collect all the pole contributions in 
Eqs. (291) and (292). Physically, this corresponds to neglecting the process of 
escape of the tunneled electron into another lead. With the help of Eqs. (G.6), 
(G.8), we find 
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where a;„ is the fermionic Matsubara frequency, and Green function for the 
open dot Qo, is defined in Eq. (109). This gives with the help of Eq. (108) 
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where the dimensionless function /(r) is defined in Eq. (304), and the Mat- 
subara counterpart of the tunneling density of states i^rir) is related to the 

38 "VVe remind that the Matsubara Green function is expressed in terms of its retarded 
and advanced counterparts as Q{iujn) = d{^^n)G^ii^n) + 0{—uJn)G'^{i'jJn) 



132 



exact tunneling density of state of noninteracting system lyri^), see Eq. (108), 

by 



sgnr + tanh — 



(354) 



at — 1/T < r < 1/T. It is noteworthy, that the inelastic co-tunnehng contribu- 
tion is expressed in terms of the Green functions of the open dot, rather than 
that of the closed dot in accordance with the discussion in the previous sub- 
section. In the absence of the interaction, the prefactor /(r)//(0) equals unity 
and Eqs. (353), (354), and Eq. (286) reproduce the formula for non- interacting 
model (107). 

It follows from Eqs. (304), (353), and (354) that n(T) is an analytic function 
of r for < Rer < 1/T, so that we can use Eq. (286) for the tunnehng 
conductance. For T <^ A^ch-E'c we find 
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where C ~ 0.577 is the Euler constant. Here we used the approximation (305), 
justified at T < N^hEc- 

To find the averaged conductance, we notice that (ut) = 1/(MA). Then, the 
integration in Eq. (355) immediately yields 



(Gin) = G^ 



V^r(i + i/7v^h) 

2r (3/2 + l/7V^h), 



(356) 



where r(,x') is the Gamma function. Equation (356) is valid for refiectionless 
point contacts, provided that tt^T ^ A^chA. 

Note that Eq. (356) predicts a power-law temperature dependence of the tun- 
neling conductance. The power-law factor can be understood in terms of the 
well-known orthogonality catastrophe [136]. An addition of an electron with- 
out a re-distribution of charge between the dot and lead would cost an amount 
of energy ~ Ec- It means that the tunneling spectrum would have a "hard" gap 
with the width of the order of charging energy, and the temperature depen- 
dence of the conductance would display an activation behavior. The power-law 
(356) indicates that a redistribution of charge indeed occurs when an addi- 
tional electron tunnels onto the dot. In fact, each channel of the ideal dot-lead 
junction carries away charge e/Nch- According to the Friedel sum rule, it cor- 
responds to an additional phase shift S — n/Nch for the scattering in each of 
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the channels. Tunnehng of an electron with energy e <^ Ec creates an excited 
state of the electron system consisting of the dot and lead. The overlap of this 
state with the initial state of the system defines the probability of tunneling. 
The smaller the energy e, the smaller such overlap is (orthogonahty catastro- 
phe [136]). This overlap can be expressed in terms of the phase shifts [136,153], 
and is proportional to (e/Ec)^, where x = Z](^/7r)^, and the sum is taken over 
all the modes. In our case, the charge is "shifted" through Nch channels inside 
the dot, and through the same number of channels in the lead, yielding 2A^ch 
equal contributions to the sum; % = 2/ Nch- This exponent indeed coincides 
with the one in the power law (356). 

A particular case of Eq. (356) for one-channel junctions deserves additional 
discussion. For spinless fermions, N^h = 1, one obtains 



We see that in this case the temperature dependence is similar to that for a 
strongly blockaded dot, see Eq. (179). The dependence can be obtained 
from the usual phase space argument for the quasiparticle lifetime, and is a 
manifestation of the Fermi liquid behavior of a one-channel dot, see Sec. 4.4 
for a qualitative discussion. The temperature dependence of the tunneling 
conductance for the electrons with spin, A^ch = 2, is qualitatively different 
[40]: 

(Gin) = Gif^^V ^ 2. (358) 



The temperature exponent here is smaller than the one given by the phase 
space argument. This indicates non- Fermi liquid behavior, which we have al- 
ready discussed in Sec. 4.2. If there is no scattering in the junction (r — 0), 
and, in addition, one neglects backscattering from within the dot (correspond- 
ing to the limit A/Ef, — > 0), then the result (358) is valid at arbitrarily low 
temperature. The existence of a finite backscattering, however, restores the 
Fermi-liquid behavior at sufficiently low temperature. We have already dis- 
cussed in Sec. 4.2 the mechanism that recovers the Fermi liquid behavior in 
the case of r ^ 0. Reflection is a relevant perturbation; at low energies it 
becomes stronger, reaching r^s ~ 1 at £ ~ To(A/'). The characteristic energy 
scale To(A/') is given by Eq. (221). At smaller energy scales the dot effectively 
is weakly coupled to both leads, and the strong Coulomb blockade behavior 
[cf. Eq. (179)] is restored. The Fermi-liquid result 
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matches Eq. (358) at the upper hmit of the allowed temperature interval. We 
refer the reader to the original reference [40] for the rigorous derivation of the 
numerical coefficient in Eq. (359). 

We see that inelastic contribution vanishes at T = in a close resemblance 

of the strong Coulomb blockade, see Sec. 3.2. Pursuing the analogy with the 
strongly blockaded regime further, one may expect, that there should be the 
another contribution analogous to the elastic co-tunneling for the dots weakly 
connected to the leads, and we turn to the analysis of this mechanism now. 

Elastic tunneling. - To obtain the leading elastic contribution we calculate 
correlator Eq. (292) taking into account all terms which include a branch cut 
in the complex r plane and two additional branch cuts: the cuts in the complex 
Ti and T2 planes in the correlation function (G.7), and similarly in higher order 
correlation functions, see Eq. (G.8). After analytic continuation (286), see also 
Appendix I, we find for A^ch > 2 
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(360) 



where Goit)^'"^ are the exact Green functions for the open dot without inter- 
actions, see Eq. (109). 

For the case of a one channel contact (Ach = 1; spinless fermions) an ad- 
ditional contribution appears, which is reminiscent of the Coulomb blockade 
oscillations 
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Similar to the discussion of the case of two-terminal conductance, we apply 
Eq. (360) to the simplest model situation first and only then consider the 
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statistical properties of the conductance. Imagine, there is a short semiclassical 
trajectory connecting the tunneling contact and the lead: 
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where tg is the time for an electron to travel along this trajectory and 7 is 
the constant characterizing the coupling of the tunneling contact with this 
trajectory. Substituting Eq. (362) into Eq. (360), we obtain for T <^ Ec 

ri ri \ 1' teNchEc < 1, roao\ 



where we omitted the non-essential numerical factors. 

Equation (363) indicates that the contribution of the processes where the time 
that the an electron spends in the dot time is smaller than the recharging time 
^/{NchEc) is not affected by the interaction at all. This was to be expected, 
because during such a short time the other electrons do not manage to re- 
distribute themselves. On the other hand, the longer trajectories experience a 
power law suppression. The physical reason for this is the orthogonality catas- 
trophe with the index x = which we already faced in the discussion of 
the inelastic contribution. 

We are ready now to perform the statistical analysis for the tunneling conduc- 
tance of the system. We will limit ourselves to the case of a single mode point 
contact and electrons with spin (i.e., A^"ch = 2; results for spinless fermions can 
be found in Ref. [70]). To find the average conductance we use Eq. (124), and 
the fact that for reflectionless contacts i/TrVFVl^''^ = Tr^iVchAM. Integration in 
Eq. (360) is then easily performed and gives ior A <^ T <^ Ef. 



where C ~ 0.577 is the Euler constant. At very low temperatures, T should be 
substituted by A\n{Ej A). The reason for the logarithmic divergence in Eq. 
(364) can be understood from Eq. (363): the contribution of trajectories with 
characteristic time te scales as l/tg, therefore the sum over all the trajectories 
logarithmically diverges. 

Similar to the case of the strong Coulomb blockade, the elastic co-tunneling 
is a strongly fluctuating quantity, 

i^Gl) _ 2 

W)^ ~ ^^^^^ 
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where (5 — 1(2) for the orthogonal (unitary) ensemble. 

If a finite backscattering in the channel is included, the energy scale (221) 
appears. Below this scale, the Fermi-liquid behavior of the system is restored, 
see also Eq. (359), so that the low temperature hmit of the elastic co-tunnehng 
becomes an oscillating function of the gate voltage jV, 



The relationships (365) for mesoscopic fluctuations still holds. 

For the case of a reflectionless junction, the periodic-with-A/" oscillations ap- 
pear only in the next order perturbation theory in the scattering matrix S. At 
large temperatures T > Aln(£'c/A) they are smaller than the aperiodic part 
of the fluctuations. Calculation of the periodic part of the fluctuations gives 



It is important to emphasize that this part is more sensitive to the magnetic 
field than the aperiodic part: it reduces by a factor of four rather than by a 
factor of two. (Noninter acting models for the quantum dot predict the aperi- 
odic fluctuations only, and hence fail to describe this anomalous magnetic fleld 
dependence.) This prediction was checked experimentally in Ref. [154]. The 
analysis of the Coulomb blockade oscillations at i? = (Orthogonal ensemble) 
and at -B = lOOmT (Unitary ensemble), see Fig. 24, yielded a suppression fac- 
tor 5.3 for the fluctuations of the conductance. This is somewhat larger than 
the predicted factor, for reasons not yet understood. 

Also note that the average conductance in the case of a symmetric weak- 
tunneling setup [panels (c) and (d) of Fig. 24] grows with the application 
of the symmetry-breaking magnetic fleld. This is in agreement with the low 
temperature RMT predictions discussed in Section 3.1, see Eqs. (138)-(140). 
If the two junctions to the dot are of the same average conductance, then 
the crossover from the Orthogonal to the Unitary ensemble should result in 
the increase of (Gpeak) by a factor of 4/3 at low temperatures. More detailed 
measurements of the magnetic fleld dependence of the average peak heights 
were reported recently in Ref. [131]. 
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Fig. 24. Conductance showing Coulomb blockade oscillations as a function of gate 
voltage Vg in the open-channel regime [(a), (b)] and weak-tunneling regime [(c), (d)] 
at i? = OmT, and B = lOOmT (the crossover between the Orthogonal and Unitary 
ensembles occurs at ~ 20mT for this quantum dot), from Ref. [154]. Note that 
the oscillations in the open-channel case are stronger at = OmT compared to 
B = lOOmT, unlike the weak-tunneling Coulomb blockade. 



5 Conclusions 



In the beginning of this Review, see Subsections 2.1-2.3.3, we have derived 
the model Hamiltonian of an isolated dot. One of the main assumptions in 
that derivation was that the electron states are random. Randomness of the 
one-electron states, a large dimensionless conductance inside the dot {g = 
Et/^ S> 1), and relative weakness of the electron-electron interaction (r^^ 1) 
allowed for a relatively simple description of an isolated dot.[33 The main part 
of the Hamiltonian of the dot is universal; it consists of a free-electron part 
described by Random Matrix Theory and an interaction term, which is not 
random, see Eq. (93). Corrections to this universal part of the Hamiltonian 
are small at 5^ 3> 1. In the discussion of the transport phenomena and of other 
effects which result from the contact of the dot with external electron reser- 
voirs (leads), we have considered that universal part of the Hamiltonian only. 



This simple model is not applicable to dots with a small number of electrons 
where the dimensionless conductance is small. It is also not applicable for dots of 
integrable shape [9], where the shell structure of the orbital states is dominant. 
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Despite this extremely simple description of the physics inside the quantum 
dot, a large variety of intricate physical effects is uncovered when the dot is 
coupled to the outside world via point contacts. These effects range from quan- 
tum interference transport phenomena in the regime strong Coulomb blockade 
(Subsections 3.1 and 3.2) to the effects of quantum charge fluctuations and 
non- Fermi liquid behavior of partially open quantum dots (Subsection 4.4). 
Remarkably, the low temperature behavior of transport through a quantum 
dot with a finite spin is always dominated by the Kondo effect, see Subsec- 
tions 3.3 and 4.3. 

In this concluding Section, we will mention three areas of research on quantum 
effects in quantum dots that are not covered in this review. 

The first area is related to the problem of electron relaxation in a constrained 
geometry. It is well-known that in a bulk clean Fermi liquid the relaxation 
rate of an electron with energy e is l/r^ ~ e^/Ep. In a confined system, an 
electron occupying a discrete level with energy e can decay with the excitation 
of an electron-hole pair. It was shown in Ref. [155] that the relaxation rate 
l/r^ ~ lS.e^ IE\ still scales as oc with the electron energy.[^ The quadratic 
dependence of the inelastic relaxation rate on the the electron energy reflects 
the smallness of the phase volume of the final state, and, therefore, is the 
same for the quantum dots and for the clean Fermi liquid. The proportionality 
coefficient is, however, much larger than its value for the clean system, which 
reflects the stronger correlation of one-particle states of the disordered dot 
in comparison with the plane waves. Nevertheless, l/r^ <^ A for e <^ -Et, 
which justifies the one particle approximation as the starting point for isolated 
quantum dot. 

The treatment of Ref. [155] neglected the discreteness of the spectrum of 
the final states. The question of the validity of such an approach was first 
addressed in Ref. [67] where it was shown that the approximation of the one 
particle density of states by the discrete levels broadened to hjT^ is correct 
only at energies larger than E* ~ ^JE^^K. At such energies, the three particle 
state obtained as a result of one-particle decay will decay further into five- 
particle state, etc., thus resulting in a true decay, without the possibility of a 
return to the initial state. If the energies of the initial state is smaller than E* ^ 
the three particle state is not strongly connected to the five particle state, so 
that the quantum dot will return to the original state. In that case, the level 
is broadened much less than it is found in Ref. [155] .[^ 



This result can be obtained by considering the second term of Eq. (93) in second 
order perturbation theory. 

Clearly, in a finite system each "broadened" one-electron level is in fact a bunch 
of many discrete many-body levels enveloped around that one-particle state. By the 
one-particle level width one means the width of that envelope. 
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Inelastic processes within the dot affect the conductance through the dot. 
This is exemphfied in the temperature dependence of the conductance of an 
open quantum dot, Eq. (349). It follows from this formula, that the weak lo- 
calization correction (second term) does not depend on temperature. This is, 
however, a consequence of neglecting the inelastic process within the dot. Only 
due to the those processes, the phase coherence inside the dot is destroyed. 
Quantitatively, it is described approximately [156-158] by the replacement 
+ 1 — >■ N^^ + 1 + 27r/ (t<^A) with equal to the inelastic relaxation time 
for an electron with energy T. Prom the theory, one expects r,^ oc T~^. Exper- 
iments with open quantum dots [159] give the result oc T~^. The reasons 
for such discrepancy have not been identified as of yet. 



The second area of research omitted in this review is related to the spin- 
spin interactions. These effects become significant at ~ 1, and may result 
in a formation of a dot ground state with spin exceeding 1/2. Sufficiently 
strong electron-electron interaction (large r^) should result in a complete spin 
polarization (Stoner instability). However, fluctuations of the density of orbital 
levels may result in an anomalously large spin of a finite-size sample even 
at moderate values of r^. This problem was studied in Refs. [76-78]. The 
simplest manifestation of the anomalous spin polarization is the transition 
from a singlet to triplet state in a dot when the number of electrons is even. 
Such a transition can be achieved by a decrease of the spacing between the 
doubly-occupied level and the next empty orbital state. In very small quantum 
dots, the level separation can be controlled experimentally by application of a 
magnetic field that affects the electron orbital motion [160]. The singlet-triplet 
transition can be seen experimentally in a measurement of the conductance 
through the dot [161] as it results in the increase of Kondo temperature at the 
transition point [162,163]. In larger dots, it can also be seen by the parametric 
dependence of Coulomb blockade peaks on the magnetic field, which shows 
a kink at the single-triplet transition [77], or a change in the peak-to-peak 
correlations from neighboring peaks to next-neighboring peaks [164]. Sofar, the 
effect of the anomalous spin polarization on the conductance was considered 
only in the regime of weak tunneling. The counterpart of this effect for the 
strong tunneling regime has not been studied. 



Finally, we have not discussed the effects of a time-dependent change of the 
shape of the dot by external gates; this effect is described by the first term in 
Eq. (93). We mention here only the suppression of the weak localization cor- 
rection in quantum dots and generation of a dc-current by a time depenedent 
perturbation [1 65-1 69] . 
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Appendices 

A Correlation of the wave-functions in ballistic dots. 

Consider a system confined by a potential U{f) that is smooth on the scale a 
much larger than the Fermi wavelength. The Hamiltonian function Ti.{p, r). 



serves as a classical counterpart of the Hamiltonian of the system. The classical 
diffusion in space, Eq. (11) is replaced by the classical evolution on the energy 
shell 7i{p,r) = Ep. This classical evolution is described by Perron-Frobenius 
equation 



where n is the unit vector in the momentum direction. The last term on the 
left hand side of Eq. (A.l) is necessary for the regularization of the otherwise 
singular eigenfunctions of the Perron-Frobenius operator. The calculation of 
the quantum corrections {e.g., weak localization) to the chaotic dynamics 
requires a finite value of l/r ~ %/{ma'^), [157,170], however, for the leading 




+ U{r), 



(A.l) 




(A.2) 
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in 1/g approximation the limit r — > oo should be taken only in the end of the 
calculation. 

The normahzation condition (48) is replaced with 

j drdnfm{n, r)fm'{-n, r) = ildSmm' , (A.3) 

where fi^ is the surface area of the unit sphere in d-dimensional space, Q2 = Stt, 
Q3 = 47r. 

The product of two Green functions, compare with Eq. (47), averaged over 
energy is given by, see e.g. Ref. [157], 

{g'^ie + u; ft, f^)Q^ie; f[,f^)) (A.4) 



Qd (-ii^ + 7m) 

r^2 = Ri,2 ± 

Repeating all the steps in the derivation of Eq. (50) from Eq. (39) we obtain 
instead of Eq. (50) 



{H^^ffs) = ciA^ {SasSp-y + Sa-yS^s) ; (A.5) 



ci = lim — 

7m 7^0 



Re7i 


r dni 


dn2 f 


7m 


^d 


^d J 



dffm{ni,f - ni6r)fm{n2, r) 



where the dimensionless conductance of the dot g is defined similarly to 
Eq. (15), 

Re7i 

It is important to emphasize that, unlike in the diffusive case, the expression 
for the matrix element involves also the eigenfunctions of the classical evolution 
operator. This can be traced to the difference of the normaliztion conditions 
(48) and (A.3). 

Analogously, Eq. (51) acquires the form. 



Vc2A^f-V, (A.6) 
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a Y p 




a Y P 



(a) (b) 

Fig. B.l. Diagrammatic expansion for the renormalization of the matrix elements 
in the Cooper channel: (a) The first order correction; (b) summation of the leading 
logarithm series for the renormalized coupling constant. Intermediate states 7 are 
only those lying outside the energy strip Et- 




X / dridr2fm{ni, n - ni6r)fm{n2, r2)fmf{n3, r2 - nsSr) fm,{n2, n) 



B Effect of the Interaction in the Cooper channel 



As we already emphasized in Sec. 2.3, the universal description of the inter- 
action by Hamiltonian (35) is valid provided that the coupling constants in 
this Hamiltonian are calculated with taking into account virtual transitions 
to the excited states which are beyond the energy strip in which the Hamilto- 
nian (35) is defined. For a weak short range interaction, these transitions only 
insignificantly change the interaction constants and Js for the charge and 
spin channels, respectively. The situation is different for the Cooper channel, 
where the renormalization is significant even for a weak interaction. 

Consider the lowest order perturbation theory correction to the interaction 
matrix element SM^app shown in Fig. B.la. 

The analytic expression for this correction is 

SHaa/S/S ~ — E — «"7'? ^^^ ^(g^gg), (B.l) 

According to the structure of the matrix elements discussed in Sec. 2.3.1, the 
largest contribution comes from the terms with 7 = 5, and one finds 

where e* is a high energy cut-off above which the approximation of the in- 
stantaneous interaction and equidistant spectrum becomes non-applicable. Be- 
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cause of the large logarithmic factor here, the renormalization of the coupling 
constants in the Cooper channel is significant, even if the interaction is weak. 

To take into account this renormalization, one has to sum the leading logarith- 
mically divergent series shown in Fig. B.l b. This results in the renormalization 
of the interaction constant in comparison with its bare value, 

Jc^l^ r^^TTTT- (B.2) 



One can see, that for repulsive interaction, Jc > 0, the effective interaction in 
the Cooper channel is always weak, and for the majority of effects it can be 
neglected form the very beginning. For the attractive interaction, < 0, the 
interaction constant renormalizes to strong coupling limit and diverges at the 
energy scale A^. — £*e~^/l'^'=l, which is nothing but the superconducting gap in 
BCS theory. If the inequaltiy Ac < Et holds, one can still use Eq. (35) for the 
description of low-energy physics of a dot; in particular, effects of finite level 
spacing on superconductivity can be accounted for in the formalism described 
here. On the other hand, if Ac > Et, the random matrix description for the 
interaction effects is not applicable. 



C Scattering states and derivation of Eq. (99) 



In this appendix we discuss the precise definitions of the scattering states ijjj (k) 
in the lead Hamiltonian Hi of Eq. (96) and the derivation of the relation (99) 
between the scattering matrix S and the matrices W, H, and U parameterizing 
the Hamiltonian of the dot-lead system. Our discussion follows that of Refs. 
[86,88-91]. 

We illustrate a part of the Hamiltonian (96) corresponding to one of the leads 
in Fig C.l. For each particular lead, we take the direction towards the dot as 
positive X, and locate the lead-dot interface near x = 0. The electron operators 
then take the form 



where $j(rj_) characterizes the structure of the wave function in the transverse 
direction. The unitary matrix Uji describes the scattering in the absence of 
any coupling to the quantum dot [86,91]. If the contact between the lead and 
the dot is defined as in Fig. C.l, with Dirichlet boundary conditions at the 
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dot 



lead 

\ ^ X 

x=0 

Fig. C.l. Schematic picture of the dot with one of the leads. 

lead-dot interface, electrons pick up a phase shift tt at the lead-dot interface 
at a; = 0, but otherwise remain in the same transverse mode, hence Uji — iSji 
in that case. 

In the absence of electron-electron interactions in the dot, the wave function 
in the leads at the energy e — vpk acquires the form 

3 

The amplitudes of the out-going waves a^"* are related to the amplitudes of 
in-going waves by the A^ch x -^ch scattering matrix S{e), 

iVch 

a^' = J:S^J{e)af. (C.3) 
i=i 

The matrix S is unitary, as required by particle conservation. In order to ex- 
press the scattering matrix S in terms of the Hamitonian of the closed dot, the 
coupling matrix W and the unitary matrix U, we represent the Hamiltonian 
(95) in terms of a Schrodinger equation for the wavefunction ipjix), which 
corresponds to the Fourier transform of the fermion fields over the entire 
real axis, 

= / ^i•(A;)e-^^ J = 1, . . . , iVeh. (C.4) 

The operator ipj{x) is a mathematical construction that is related, but not 
identical to the true electron operator ■0e(r) defined in Eq. (C.l). Then, de- 



^in ^i{kF+e/vF)x _|_ ^out ^-i{kF+£/vF)x 



(C.2) 
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noting the wavefunction inside the dot as 0^, we find from Eqs. (92) and 
(95)-(97), that the Schrodinger equation for the entire system reads 



'^^ i/=i 



v=l 1=1 



(C.5) 



For a; 7^ 0, Eq. (C.5) corresponds to a left moving particle at velocity vp- To 
find the scattering matrix 5^^, we set, cf. Eq. (C.2), 



A^ch 

1=1 

Nek 



(C.6) 



where k = e/vp- At a; = we use the standard regularization iljj{0) = 
[ipji—O) + 'V'i(+0)]/2- Substitution of Eq. (C.6) into the Schrodinger equation 
(C.5) yields 



V j.l = l 



The wavefunction 0^, of the dot can be eliminated from this equation, resulting 
in the condition (C.3) for the amplitudes a^f and a°"*, with the iVch xiVch matrix 
S given by 



S{e) = U 1- 2mvW^ {e-n + mvWW^) ^ W 
This is Eq. (99) of the main text. 

The coupling matrix W can be represented in the form 



W = J^VOW, (C.7) 

where V is an orthogonal M x M matrix, is a real A^^ch x A^h matrix, and O 
is an M X A^^h projection matrix, Oaj = S^j, a — 1, ■ ■ ■ , M and j — 1, . . . , Ach- 
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The first factor in Eq. (C.7) ensures a convenient normalization of the matrix 
W. Because the distribution function of the random matrix Ti is invariant 
under rotations Ti. — > VTiV^ , the matrix V in Eq. (C.7) can be omitted. The 
matrix W can be related to the {2Nch) x {2Nch) scattering matrix Sc of the 
point contacts between the lead and the dot. Unlike the scattering matrix 
S{e) which connects the incoming and outgoing states for the entire device, 
the matrix Sc characterizes the properties of each point contact connecting 
a lead to the dot separately. The corresponding transmission and reflection 
amplitudes are defined in A^ch x A^ch matrices Vc, r'^, and tc, 



which can be expressed [92] in terms of the coupling matrix W and the matrix 



The matrix Tc expresses the reflection from a point contact to the leads, 

corresponds to the reflection of an electron entering the point contact from the 
dot, and the matrix tc describes the transmission through the point contact 
from a lead to the dot. A contact is called ideal if Vc = r'^ = 0, or, equivalently, 
W^W = 1. 



D Mesoscopic fluctuations of elastic co-tunneling far from the peeiks 

The calculation of the magnetic field dependence of the mesoscopic fluctua- 
tions of elastic co-tunneling is more complicated than in the pure orthogonal 
and unitary ensembles, as in the crossover between the orthogonal and unitary 

ensembles the electron and hole amplitudes become correlated. This is because 
the wavefunctions for different energies are not independent in the crossover 
regime. The corresponding correlation functions are obtained from Eqs. (149) 
and (25), 




(C.8) 



U of Eq. (C.l), 




(C.9) 



Ee + Eh 

A 



A 



Ai 



( 



ATfA 




{FeFh) 



Ee + Eh 



A 



1 



( 



27rEh 
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where the dimensionless functions Ai{x) is given by 



Ai{x, y) = ^(^ij^^2^ |-g- + y — ^ + InxyJ arctanary - a;Li2(l - y) 

- i {L. + . + m (1 + .V) - Re [(. + x)L„(-,x,)]" 

(D.2) 

with Li2(x) being the second polylogarithm function [103]. The asymptotic 
behavior of functions Ai is Ai(a;, y) = (xy/n) \ia{xy) for a; <C 1, and Ai(a;, y) ~ 
|ln^(xy) — In^x /{ttx) for a; S> 1. The hmits of the pure orthogonal (uni- 
tary) ensembles (151) are recovered by putting A^^ = and Nf^ = 0(oo) in 
Eq. (D.l). 

The correlation function of the conductances at different values of the magnetic 
field is somewhat involved, 



[Ee + Ej, \27rE, 

EhEe . (NP^ Eh , 
+ ^^^^Ai — — , 1 + — \ + E^^ Eh 



In this general case, where the gate voltage can be anywhere in the Coulomb 
blockade valley, it is impossible to point out the single scale for the magnetic 
field. However, in the vicinity of the peaks, E^ ^ Eh or E^. <^ Eh, the situation 
simplifies significantly, and one obtains the universal dependence (160). 

The dependence of the crossover parameter A from Eq. (162) on the magnetic 
flux $ in our case is 



\^ . ^ ' - - ' {E^ + Ehf 





(D.4) 

where the functions A, Ai are defined in Eqs. (159) and (D.2), and N^^ is 
defined in Eq. (20) with $i = $2 = ^- As before, there is no universal depen- 
dence on the magnetic field. In the vicinity of the peak, however, we recover 
the universal result (163). 
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E Derivation of the Hamiltonian (238) 



To include the second junction we should modify Hamiltonian given in Eqs. (207), 
(209) and (212): in the charging energy, Eq. (209), the variable Op{0) should be 
replaced by the difference 9p{0)—9p{L)] we should also return from Eq. (212) to 
the generic form (206) of the free-field Hamiltonian, and replace the boundary 
condition Eq. (213) by a barrier described by the Hamiltonian 

Hri = --|ri|Dcos[2v^^p(L)]cos[2v^^,(L)], (E.l) 



which is similar to Eq. (207). Similar to the one-junction case [cf. Eq. (211)], 
in the energy range E <^ the fluctuations involving the total charge of the 
dot are pinned, 

^(MO)- W]). = eAr, (E.2) 



and fluctuations of the fleld Op have the same phase at both junctions. If, on 
the other hand, E ^ A, this is the only constraint on the fluctuations of 
the charge and spin flelds at two different junctions. Therefore, in the energy 
domain Ec ^ E ^ A, there are three independent modes left: the spin modes 
6s and O^i at each of the junctions, and one charge mode 6p corresponding to 
passing of the charge across the entire dot. It is more compact to describe the 
system in terms of the effective action [140], rather than a Hamiltonian, 



+ 4 



-ErD j dT[\r\ cos[7r7V - 2v^^p(r)] cos[2v^^,(r)] 

+ |ri| cos[2v^^p(r)] cos[2v^0,i(r)]} , (E.3) 

here r is the imaginary time, and Un is the Matsubara frequency. The spin 
degrees of freedom at each junction are described hy g = 1/2 Luttinger liq- 
uids, and the mode corresponding to the remaining charge degree of freedom, 
is a Luttinger liquid with (yf = 1/4; its increased stiffness comes from the low- 
energy constraint Eq. (E.2) of the charge neutrality of the dot. It is easy to 
check that because of this increased stiffness [i. e., because of {g < 1/2)], 
the reflection at the point contacts is a relevant perturbation.j^ If the bare 



The quadratic part of the effective action (E.3) represents tfie fixed point of a 
four-channel Kondo problem[40], and the reflection represents relevant perturbation 
near the fixed point. 
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reflection amplitudes |r| and |ri| are smafl, then with the reduction of the en- 
ergy scale E they grow as (E^/EY^^, independently of each other. We intend 
to consider the asymmetric geometry: the bare reflection amplitudes are very 
different: |r| <S |ri|. Then, the crossover to weak tunneling through the x — L 
junction is either occurs upon reaching the energy scale ~ £'c|ri|^, or tunnel- 
ing is weak from the very beginning. At energies E ^ i?c|?"i|^, the variables 
dsi and 9p are fixed most of the time, so that cos(2y^^^p) cos{2y/Ti6 si) = 1. to 
ensure the energy minimum of the junction. The system makes rare hops be- 
tween various configurations satisfying the energy minimum condition. Upon 
further reduction of the band width, the effective backscattering at x = 
follows the renormalization prescribed by Eq. (220), and reaches the limit 
of the weak tunneling at energies ~ Ec\rir cos ttN'\'^ . This energy scale co- 
incides with Eq. (221) if |ri| ~ 1, and replaces it if |ri| is small. Note that 
similar to the condition (225), the spatial quantization of the electron states 
in the dot does not affect the crossover to the regime of weak tunnehng, if 
A -C £'c|rir cosTrAAp. In that regime, the system performs rare hops between 
the minima of the "potential" 



UiOs, Osi, Op} = - Wcos(7rA/' - 2^ep) cos(20F^,) 

-Wicos(2v/^^p)cos(2V7r^,i), (E.4) 

represented by the last two lines in the action (E.3). [Here the renormalization- 
dependent energies U > and Ui > replace the corresponding bare values 
entering in Eq. (E.3).] In the minima of the potential (E.4), the spin of the 
dot {^/^T/2){6s — 9si) is integer or half-integer, depending on the value of A/", 
cf. Eq. (224). There are three types of hops which connect the minima of the 
potential (E.4). Two of the types, ^,(0) ^ ^,(0)±V7r and ^,i(0) ^ ^,i(O)d=0F, 
are already familiar to us from the single- junction case, see Equation (226) 
and the discussion which follows that Equation. These hops correspond to a 
change by 1 of the 2;-projection of the dot's spin via spin exchange with one of 
the leads. The third type of hops involves a simultaneous change of all three 
fields: 

es{o)^e,{o)±^/^/2, ^,i(o) ^^,i(o)± V^/2, ^^(o) ^ ^^(o) ± v^/2. 



In such a hop, one electron (charge e and spin = ±1/2) is transferred across 
the dot; the spin of this electron and the dot's spin may flip in this process. 
The corresponding tunneling Hamiltonian, 

^ oc - cos (^4^ cos {^^ij cos (^0f^ , (E.5) 
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supplements Eq. (226) and a similar equation for the x — L junction, 



Here = 0si(a: = L + 0)— = L — 0) and 0^ = 0p(a; = +Q) — 4>p{x = — 0) 
are the discontinuities of the respective fields at the junction (for the field 
(pp the discontinuity is the same at both junctions). All three Hamiltonians, 
Eqs. (E.5), (E.6) and Eq. (226), have the same scaling exponent, as can be 
easily checked with the help of the action (E.3). These three Hamiltonians 
represent the easy-plane part of the exchange interaction. The SU{2) symme- 
try guarantees the existence of terms like Eq. (227), which restore the isotropy 
of the exchange. The exchange interaction changes the spin of the dot by an 
integer. If we account for the finite level spacing, then the lowest energy states 
would correspond to the smallest possible spin of the dot, like in the one- 
junction geometry considered above. Kondo effect develops only if the spin 
of the dot is 1/2. This doublet state is realized in the dot periodically with 
the gate voltage, when costtA/" < (we assume here |ri| ^ |r|). In Subsection 
4.3) we concentrate on the doublet state only. Returning to the fermionic vari- 
ables at £^ <^ £'c|rir costtA/'P, we find the exchange Hamiltonian (238) which 
generahzes Eq. (228) to the case of two junctions. 



F Canonical versus grand canonical ensembles 

A crucial step in the effective action formalism is the redefinition (254) of the 
charge operator n in terms of the fermion operators in the leads. As was dis- 
cussed below Eq. (254), strictly speaking, such a redefinition requires a canon- 
ical description of the entire system, i.e., that the total number of particles in 
the leads and the quantum dot is kept fixed in the ensemble. A grand-canonical 
approach, as in Sec. 4.5.1, can be used only when all observablcs are consid- 
ered at a single value of the gate voltage A/", or when physical observables are a 
periodic function oiM . For a non-periodic A/'-dependence (i.e., for mesoscopic 
fluctuations of the differential capacitance) , a canonical description has to be 
used. In this section, we discuss implementation of the canonical constraint in 
the grand-canonical theory. In particular, we show, within the effective action 
formalism, that for T ^ N^h^ or A'^ch ^ 1, the canonical constraint can be 
accounted for in the grand-canonical description by the dependence (326) of 
the chemical potential n on the gate voltage M. 

To go between the canonical and grand-canonical descriptions, we use the 
inverse Laplace transform that relates the thermodynamic potential f^ATp in a 
canonical ensemble with Np particles to the thermodynamical potential f2(//.) 
in a grand-canonical ensemble with chemical potential /x. 




(E.6) 
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J 2mT ^ ' 

-TzTi 

Here "Tinj, denotes a trace over all states with a fixed number of particles 
Np, Tr is a trace over all states with all numbers of particles, and Np is the 
operator for the total number of particles. 

Hence, what we need is the dependence of on a complex chemical potential 
/i in a strip ttT < Im /i < ttT. The modification of the effective action theory of 
Sec. 4.5.1 to include a complex chemical potential is straightforward, because 
all intermediate steps are analytic in for ttT < Im ^ < ttT. In the final result, 
the //-dependence appears through the /x-dependence of the scattering matrix, 
S{ijL,t) — > S{0,t)e~'^^* and through the //-dependence of the thermodynamic 
potential of the lead. To find the latter, we include /j, into the charging energy. 



oo 

1 — oo 

-/,AA-//(iV,,^)„ (F.2) 



c 



where M' = M + /i/(2£'c) and Np^L is the number of particles in the leads. 
Hence, by the second fine of Eq. (F.2), wherever we found an exphcit M- 
dependence, H should be replaced by Af+ji/ {2Ec). The term fi{Np^L)q denotes 
the quantum mechanical average of the number of particles in the lead and 
arises due to the normal ordering of the V'-fields. 

Let us now discuss the implication of this scheme for the calculation of the 
mesoscopic fluctuations of the differential capacitance, that were the subject 
of the Subsection "Non-periodic interaction corrections to the capacitance, 
Ki^{s)" . Using Eq. (332) and upon inclusion of the /i-dependent contributions 
to Q from Eq. (F.2) and the /x-dependence of the scattering matrix, we find 
that the /t-dependence of the thermodynamic potential of the entire system is 
given by 



n(/t) = o(o) 



c 

oo 



-T dt dt 



liAf - ii{Np) 

,{t + t')[l - iiijt - t') - e-"'^'-''^TrS^(t')S(t) 
4(i - t') sinh[7rr(t - t')] 



152 



to 



2 N^Y^J J J sinh[7rT(i + s)] smh[7rT(i' + s)] 



X 



f sinh[7rr(5 - to)] sinh[7rr(5 + t + t' + tp)] V 

\ sinh[7rT(t + to)] sinh[7rT(t' + to)] J ' 



where (iVp) is the expectation number of the total number of particles in the 
system (leads and dot) at /x = 0. Mesoscopic fluctuations of Vt are dominated 
by the fluctuations of the second line in Eq. (F.3), which are small as NA/T 
or l/Nch for T ^ NA or N^h S> 1, respectively. Hence, to find the leading 
//-dependence of Jl(//), one can replace Q by its ensemble average, for which 
one finds, 



(0(/.)) = (0(0)) - ^ - /M/- - pc{N,) - ^. (F.4) 

With this, the /i-integration in Eq. (F.l) is straightforward when we choose 
Np = {Np); for T 3> A^chA the integration can be done by the saddle point 
method and amounts to the substitution (326). Interaction and fluctuation 
corrections to Q can be treated systematically by an expansion of exp(— fi/T) 
around the average (F.4). The differential capacitance (334) follows from Eq. 
(F.3) by twofold differentiation to J\f. 



G Correlation functions 



Correlators offermionic operators. - To complete the review of the properties 
of the fermionic correlation function, we give the general formulas allowing 
one to compute any non- vanishing average of the fermionic operators: 



/ Nch 



k=l 



k=l 



k=l 



k=l 



g-i27rm(Ar+Woh/4) jj 



n 



^ --sin7rr(4-r^) 



k>l 



ttT 



In this case, the approximate formula (122) for the average of the scattering 
matrix gives a result that is wrong by a factor 1 + (2 — /?)/ (3N"-^^ , because time scales 
~ Nci^A, for which the approximation (122) does not hold, arc important. Exact 
evaluation makes makes use of the relation {{2Tri)~^Tr S^dS/ds)) = A"^. 
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,t-\sin7rT(r/,-r^) n ^ 



1 sm7rT{TPk-Tjf) 
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1/iVch ' 
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(G.l) 



where n^'^ are arbitrary non-negative integers, m is an arbitrary integer, Nf ~ 
NchTn + Y.ji'nf + nj^), and we use the convention Ylf=i aj — 1, for N — 0, and 
njLi = 0, for N < 0. The function /(r) and the time scale to are defined 
by Eq. (305). One can easily check that Eqs. (306), (307) (309) follow from 
Eqs. (G.l). 

Correlators of fermionic operators and current operator. - For the calculation 
of the two-terminal conductance, we need correlators that involve both the 
fermionic operators and the current operator /, which is linear in the boson 
fields. 



evp d 
27r dx 



(G.2) 



Here the diagonal matrix A is defined in Eq. (105). The necessary bosonic 
correlators are obtained from those of the free boson fields (p by substitution 
T — > T ± ix/vp, where the -|- (— ) sign is for left (right) movers. For the 
correlators of the current operators one thus finds 



\ ^ ^ ' ^ ''1 27rV27ry A^ch ^ sin2[7rr(T ± iO)] ^ ' 

The correlators of the current operators and two or four fermionic operators 
are most conveniently expressed in the correlators of the fermionic operators 
themselves. 



Nl ^_ Nl Nf, ^ Nr 

(r.7(r)/(0) n V^L.n.lr/^) n ^L,nM) li R,n, (t^) li Mr!")) , = 
i=l j=l k=l 1=1 
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N-L ^ Nl Nr ^_ Nr 

(T./(r)/(o)),(T.n^L,n,(T^) nv^^,n,(r/) n v^«,n,(T,^)n^^^^^ 

i=l j=l k=l 1=1 



Nr. 



\ / 2 — 1 j — 1 k — 1 / — 1_ 



X 



Nr. 



Nr; 



i=i 

'Nr. 



k=l 



1=1 



Nl 



Ni, 



Nr 



1=1 i=l k=l 1=1 

9n{r) = An„7rTcot[7rT(r + iO)], 



(G.4) 



where Ni, Ni, Nji, Nr are the arbitrary non-negative integers, and rii — 
1, . . . , Nch labels the channels. The expressions for the average of the fermionic 
operators were given earlier, see Eqs. (306) ~ (309), and (G.l). 

Correlators of fermionic operators and the charging operators F, . - Calcu- 
lation of the tunneling density of states Eqs. (290) - (292) involves not only 
the correlation function of the fermionic operators but also the products of 
the operators F, F^ that increase or decrease the charge of the dot by one 
electron. The corresponding calculation is facilitated by the observation that 



F(r)F(0) = e-^«ff^F^e-^«ff^F = T^exp { j dn [h^s - F^H^q{ti)F 

.0 



where the one dimensional Hamiltonian is given by Eq. (293) and the commu- 
tation relation for F, F"^ are given by Eq. (287). In bosonized representation 
one finds 



^eff ~ F^ HpffF 



eflfJ 



TT 



E IM^) + M^)] +2n{Af-m- 1/2) 



(G.5) 



Because this term is linear in the bosonic fields, further calculation proceeds 
without any difficulty with the help of Eq. (303). The operator m in Eq. (G.5) 
does not have its own dynamics and can be fixed to be any integer, e*^'^'" — 1. 



We obtain 



(r.F(T)F(O)), 



/(O) 



2/iVch 



(G.6) 



where the dimensionless function /(r) is defined in Eq. (304). 
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For the averages involving fermionic operators as well as the charge changing 
operators we find [compare with Eq. (306)], 



(27r^;^)(T.F(r)F(0)V;.,(rO^L,,(r2)),= 



f{r) 



/(O) 



( f{r-n)f{-T,y 

sin[7rT(Ti-T2)] V/(-^i)/(^-^2). 



l/JVch 



(G.7) 



fir) 



/(O) 



2/Weh 



( f{n-r)f{T,) \ 
sin[7rT(ri-r2)] [f{n)f{T2-T)J 



l/iVch 



Notice that the structure of this expression is equivalent to the form factors 
in the problem of the orthogonality catastrophe [153]. 

Averages involving larger number of the fermionic operatos are found in the 
same manner. We will give the formula expressing correlation function involv- 
ing the charging operator and the arbitrary number of the fermionic operators 
in terms of the correlation function involving the fermionic operators only: 



i=l j=l k=l 1=1 



fir) 



/(O) 



2/Afch 



Nr. 



Nf, 



Nr, 



X 



{Tr n ^L,nM') n ^L,nArj') U i^R,nM) U^R^nXr^), 
1=1 j=l k=l 1=1 

1/iVch 



""'fir-rh^^ ^fir,^-r)^j /(rf) 



n 



n 



n 



n 



i=i fi-rt) ■=ifir-rj^)L\ f{r^) t\f{rp-r) 



(G.8) 



where Ni, Ni, Nj^, are the arbitrary non-negative integers, and nj = 

1, . . . , Nch labels the channels. The expressions for the correlation functions of 
fermionic operators were given earlier, see Eqs. (306) - (309), and (G.l). 



H Derivation of Eqs. (340) - (343) 



Here we present the details of the calculation of the current-current correlator 
n(T) and the two-terminal conductance G up to second order in the effective 
action S^s- 

To zeroth order in the interaction, n(r) is given by Eq. (G.3). The correction 
5n(r)'^^) to first order in the action follows from Eqs. (306) and (G.4) as 
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8^ ± / / sin[7rr(Ti - T2)] 
X {cot[7rT(ri - r ± iO)] - cot[7rT(ri - r ± iO)]} 
X {cot[7rT(ri ± zO)] - cot[7rT(r2 ± iO)]} . 



(H.l) 



We change to variables ri = a' + a and r2 = a' — cr and integrate over a'. 
Since the integrand is analytic in the upper (lower) half of the complex plane 
for the + (— ) signs in the denominators, one directly finds 



(H.2) 



The second order correction to 11 (r) reads 



Here we find for Net > 2 



(H.3) 



SU{t 



I (2) 



1^ 



l/T 

^ / dn dT2 drs dT^Tr S{ti - T2)AS{t3 - T^jA 
± i 



X 



X 



X 



sin[(Ti ± T + iO)7rT] sin[(T4 ± r + iO)7rT] 
1 

sin[(T2 - i0)7rr] sin[(r3 - iO)7rT] 



f{r2-Ti)f{T3-n) 
firs - ti)/(t2 - T4) 



1 l/Afci 



(H.4) 



We left out terms that have all poles at the same side of the real axis, because 
they vanish after integration over the times ti, . . . , T4. We now parameterize 
Ti, T2, T3, and T4 according to 



Ti=(7' + (7/2 + (7i, 
T2=(7' + (7/2, 

T3 = 0-' - a/2 + (72, 
T4^a' - a/2, 



(H.5) 



perform the integration over a', and use the Lehmann representation (316) for 
Sir), 
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It 



-iti+il2 



It 



•7t+ll2 

71 -it 2 



(a) 



(b) 



Fig. H.l. Shift of the integration contours for the variable a in Eq. (H.6). Branch 
cuts of the integrand are indicated by thick hues, poles by dots. The integration 
contour for first term between brackets {. . .} in Eq. (H.6) is shifted below the real 
axis (a); the integration contour for the second term is shifted above the real axis 

(b). 



8 



sm[{it2 — a)7rT] 



1 



X 



X 



[sin[(±T + iti)7rT] sin[(±T -a + iO)7rT] 

1 

sin[(±T — it2)'!rT] sin[(±r + a + iti — it2)TrT] 
f{-ih)f{ih) 



(H.6) 



Finally we shift the integration over a to the complex plane, using the ap- 
proximation (305) for /. For the first term between brackets, we choose the 
integration path below the real axis, while for the second term the integration 
path is chosen above the real axis, see Fig. H.l. The result is 



2 2 ^ 

5n(2)(r) = ^^ I dh I dt2TiS{ti)AS^{t2)A 



-oo 



X 



sinh[(ii ± ir)7rT] sinh[(-i2 ± ir)7rT] 

sm{n/Nch) 



oo 

+ T J ds 



to 



sinh[(s + ti-t2± iT)7iT] 



X 



+ 



1 



sinh[(ti ± iT)7rT] sinh[(s + ti)7rT] 
1 

sinh[(— ^2 ± iT)7iT] sinh[(s — t2)7iT] 
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^ / smli[(fi -t2 + ,s + to)nT] smli[(,s - t,)nT] \ '^^^ 1 

smh[(ii + io)7rr] smh[(-i2 + io)7rr] j ' ^ 

The first term is the pole contribution, the second term contains a contribution 
from the branch cut. Corrections to n(r) to higher orders in perturbation 
theory can be classified in the same way as terms with pole contributions, 
one branch cut, two branch cuts, etc. As discussed in Sec. 4.5.2, there are no 
corrections of higher order in the scattering matrix S with zero or one branch 
cut — these have been accounted for here; all higher order corrections have 
more than one branch cut. Now, using the conductance formula (281), one 
directly finds the interaction correction (342) to the conductance. 

For A^ch = 2 there is an additional, periodic in A/", contribution to the current- 
current correlator 11 (r) to second order in the action Scs- This origin of this 
extra A/'-dependent contribution is the fact that for Nch = 2 the product (309) 
of four different Fermion operators has a nonzero average. We thus find 



g2j'2g-27riA/' 

(5n(T)osc = g^^2^2 J dTidT2dT^du 

^ 

X /(0)-2 [f{T, - n)f{n - rs)f{n - T,)f{T, - nf" 

X [Sii{Ti - T2)S22{r3 - T4) + Sui^ - T2)S2i{t3 - T4)] 

sin[(ri - r3)7rr] 



l/T 



X 

± 



^ sin[(Ti ± T + iO)7rT] sm[{Ts ± r + iO)7rT] 

X '''^''r'f^^ , , +c.c. (H.8) 

sin[(T2 - i0)7rr] sin[(r4 - iO)nT] ^ ' 

[Note that the products 5*11 6*22 and 8x2821 appear with the same sign, because 
the latter product has an extra minus sign both in the fermionic correlator 
(309) and in the current-fermion correlators (G.4).] Again, terms that vanish 
after integration have been left out. 

We now use the parameterization (H.5) of the times ri, . . . , r4, integrate over 
cr', and implement the Lehmann representation (316), to find 



e 1 e~ f f 
m(r)osc - 3^ .^2^2 j dtidt2 J da [8u{ti)822{t2) + 8i2{ti) 821(12)] 



"0 



X /(O)-' [fHti)fHt2)f{-a - th)f{a - zt2)f' 
sin[((7 + iti — it2)7iT] 



xE 

± 



sin[((T + iti ± T)TrT] sm[{it2 ± r)7rT] 
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sin[(f7 + iti — it2)nT] 
~ sin[(-(7 + it2 ± T)nT] sm[{iti ± T)nT] 

Next we integrate over a by shifting the integration contour into the upper 
(lower) half of the complex plane for the first (second) term between brackets 
{. . .}. The only contribution is from the branch cut of the integrand, which 
reads, with the approximation (305) for /, 



+ c.c. 



(H.9) 



(5n(T)osc = j dtidt2 J ds [-Sn(ti)5'22(t2) + ^12(^1)^21(^2) 



to 

1 

X 



X 



± 



^sinh[(ii + to)7rT] sinh[(t2 + to)vrT] 
1 

^sinh[(ti +t2 + to + s)7iT] sinh[(s - to)TTT] 
sinh[(s + ti- t2)7TT] 



sin[{is + iti ± r)7rT] sin[(it2 ± t)7iT] 



+ sinh[(. + t,-t2)7rr] 1 

sin[(is + ii2±r)7rT]sin[(iti±r)7rr]J ^ ' ^ 

Finally, using the conductance formula (281), we find the oscillating contribu- 
tion (343) to the conductance. 



I Derivation of Eqs. (360) - (361) 



To obtain the leading elastic contribution we calculate the correlator (292), 
taking into account all the terms which include branch cut in the complex r 
plane and the cuts in the complex ri and T2 planes in the correlation function 
(G.7), and similarly in higher order correlation function, see Eq. (G.8). This 
gives [70] for A^^h > 2 



fir) 



/(O) 



X 



dndT2 \go{ri)iyWW^go{-T2] 



/(-ri)/(-r2) 



11 



sin 7rT{T + Ti — T2) 
/(Ti)/(r2) 



+ 



J{ri + r)f{r2-r) 



/(-Ti-t)/(t-T2). 
- 2 



(I.l) 
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It, 



it^+l/T it^ 
1/T 



it^+l/T 



1/T 2 



(a) 



(b) 



Fig. I.l. Deformation of the integration contours for the variables ri (a) and T2 (b) 
in Eq. (I.l). Branch cuts coming from functions /(r) are shown by thick hues, and 
the integration contour is represented by arrowed vertical thin lines. 

where Qo{t) is the Matsubara Green function for the open dot. One can easily, 
see that the pole contribution at r + Ti — T2 = is excluded from Eq. (I.l), 
because it is already taken into account in the inelastic contribution (353). 
One can now deform the contours of the integration over Ti and T2 in Eq. (I.l) 
using analyticity of function /(r) at Imr < 0, see Fig. I.l. 



fir) 



/(O) 
X Re J 



2/Areh 



dtidt2 



g^{h)i^ww^g^{-t2) 



fHh}fHf2} ' 



11 sm7rT{T + iti - ^^2) j ' 



(1.2) 



where Qoit)^'"^ are the retarded and advanced Green function for the open dot. 
Deriving Eq. (1.2), we use the fact that /(r+l/T) = /(r), ^(r+l/T) = -^?(r), 

and g{u + 0) - g{it - 0) = tig'^it) - g^{t)]. 

The resulting function (1.2) is analytic for < r < 1/T, and we can use the 
analytic continuation Eq. (286) to find the tunneling conductance Eq. (360). 

For the case of a one channel contact A^ch = 1 (spinless fermions) the additional 
contribution appears due to the anomalous average (309). 



GTir) = 



2 

X Re |( 



fir) 



/(O) 



2 1/T 



dT\dT2 



fir + n-T2) /(-ri)/(r2) 



-i27rA/'-m/2 



tofiO) fir2-T)fi-T-n) 

goiTi)iyww^goi-T2)]^^} . (1.3) 



After deformations of the integration contours analogous to derivation of 
Eq. (1.2) and analj^ic continuation we arrive to Eq. (360). 
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